{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "e1fc6530-2f9c-4c03-959b-b5e0a6345b43",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import healpy as hp\n",
    "import matplotlib.pyplot as plt\n",
    "from glob import glob\n",
    "import json"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4ec3d333-f75a-4efc-b066-ac17fc0de47a",
   "metadata": {},
   "source": [
    "## Files for Release"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "91e20f3d-8408-4f40-b0c2-b322a762c333",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Location in /data/ana\n",
    "ana_path = './data_files'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "cfdb8119-ba7e-45af-97bb-0e5dd5626b3d",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Available map files:\n",
      "  ./data_files/IC86_sid.fits\n",
      "  ./data_files/IC86_sid_4-4.25GeV.fits\n",
      "  ./data_files/IC86_sid_4.25-4.5GeV.fits\n",
      "  ./data_files/IC86_sid_4.5-4.75GeV.fits\n",
      "  ./data_files/IC86_sid_4.75-5GeV.fits\n",
      "  ./data_files/IC86_sid_5-5.25GeV.fits\n",
      "  ./data_files/IC86_sid_5.25-5.5GeV.fits\n",
      "  ./data_files/IC86_sid_5.5-6GeV.fits\n",
      "  ./data_files/IC86_sid_6-6.5GeV.fits\n",
      "  ./data_files/IC86_sid_6.5-100GeV.fits\n"
     ]
    }
   ],
   "source": [
    "# Map files include energy binned and total sidereal maps\n",
    "map_files = sorted(glob(f'{ana_path}/*.fits'))\n",
    "print('Available map files:')\n",
    "for f in map_files:\n",
    "    print(f'  {f}')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "1b1982af-d447-48ed-abf0-072d045bb0c6",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Energy distribution file:\n",
      "  ./data_files/energy_distributions.json\n"
     ]
    }
   ],
   "source": [
    "# True energy distribution for each energy bin\n",
    "dist_file = f'{ana_path}/energy_distributions.json'\n",
    "print('Energy distribution file:')\n",
    "print(f'  {dist_file}')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1f8daf7f-b163-474a-bf96-1d38429df616",
   "metadata": {},
   "source": [
    "## Functions for Plotting"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "a3c72cb4-384c-4155-a749-87e859bab7ee",
   "metadata": {},
   "outputs": [],
   "source": [
    "\"\"\" Mask a map based on a minimum and maximum declination \"\"\"\n",
    "def maskMap(unmasked_map, dec_min, dec_max):\n",
    "\n",
    "    n_pix  = len(unmasked_map)\n",
    "    n_side = hp.npix2nside(n_pix)\n",
    "    theta, phi = hp.pix2ang(n_side, range(n_pix))\n",
    "\n",
    "    theta_min = (90 - dec_min) * np.pi/180\n",
    "    theta_max = (90 - dec_max) * np.pi/180\n",
    "    theta_cut = (theta <= theta_min) * (theta >= theta_max)\n",
    "\n",
    "    masked_map = np.copy(unmasked_map)\n",
    "    masked_map[np.logical_not(theta_cut)] = hp.UNSEEN\n",
    "\n",
    "    return masked_map"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2f9be88d-c430-4db7-ac07-816eb90ccd6c",
   "metadata": {},
   "outputs": [],
   "source": [
    "\"\"\" 'Top hat' smoothing over a given angle in degrees \"\"\"\n",
    "def smoothMap(m, smoothing_angle_deg):\n",
    "\n",
    "    npix  = len(m)\n",
    "    nside = hp.npix2nside(npix)\n",
    "    smoothing_angle_rad = smoothing_angle_deg * np.pi/180.\n",
    "    smoothed_map = np.zeros(npix)\n",
    "\n",
    "    vec = np.transpose(hp.pix2vec(nside, np.arange(npix)))\n",
    "    for i in range(npix):\n",
    "        neighbors = hp.query_disc(nside, vec[i], smoothing_angle_rad)\n",
    "        smoothed_map[i] += m[neighbors].sum()\n",
    "\n",
    "    return smoothed_map"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "63bdffb0-1a82-4da2-9dcb-e6175163a3ab",
   "metadata": {},
   "outputs": [],
   "source": [
    "\"\"\" Li & Ma Significance \"\"\"\n",
    "def LMSignificance(n_data, n_bg):\n",
    "\n",
    "    # For every event, we generate 20 time-scrambled events for the reference map.\n",
    "    # Alpha is the ratio of the two sets of events\n",
    "    alpha = 1/20.\n",
    "    \n",
    "    n_on  = n_data\n",
    "    n_off = n_bg / alpha\n",
    "\n",
    "    sigma = np.sign(n_data - n_bg) * sqrt(2*(n_on*np.log(((1+alpha)*n_on) / (alpha*(n_on+n_off)))\n",
    "        + n_off * np.log(((1+alpha)*n_off) / (n_on+n_off))))\n",
    "\n",
    "    return sigma"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "dfc3d0e6-d890-4cba-9bbf-34b489d0da2a",
   "metadata": {},
   "outputs": [],
   "source": [
    "\"\"\" Takes a file consisting of data, bg, and local maps. Returns desired final type (relint, sig, etc) \n",
    "with smoothing, masking, and fitting options applied \"\"\"\n",
    "def getMap(map_file, map_name, smoothing_angle=0, dec_min=-90., dec_max=-25.):\n",
    "\n",
    "    # Read in input file\n",
    "    data, bg = hp.read_map(map_file, range(2))\n",
    "\n",
    "    # Option for top-hat smoothing radius\n",
    "    if smoothing_angle != 0:\n",
    "        data = smoothMap(data, smoothing_angle)\n",
    "        bg = smoothMap(bg, smoothing_angle)\n",
    "\n",
    "    # Return desired map type\n",
    "    if map_name == 'data':\n",
    "        output_map = data\n",
    "    elif map_name == 'bg':\n",
    "        output_map = bg\n",
    "    elif map_name == 'sig':\n",
    "        output_map = LMSignificance(data, bg)\n",
    "    elif map_name == 'relint':\n",
    "        output_map = (data-bg) / bg\n",
    "\n",
    "    # Eliminate NaNs and mask\n",
    "    output_map[np.isnan(output_map)] = 0\n",
    "    output_map = maskMap(output_map, dec_min, dec_max)\n",
    "\n",
    "    return output_map"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "176abefa-935b-4927-b365-9358192de22b",
   "metadata": {},
   "source": [
    "## Sample Plot"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "e94a909d-11eb-4dcf-9749-8e9a2ba78664",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "NSIDE = 64\n",
      "ORDERING = RING in fits file\n",
      "INDXSCHM = IMPLICIT\n"
     ]
    }
   ],
   "source": [
    "m = getMap(map_files[0], 'relint', smoothing_angle=5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "fc09def8-3718-49e6-8b5e-bac61d491547",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoYAAAFzCAYAAACjPchOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzsvXusdd9a3/V9xphzrb33+76/c+BAawHTpJqKtiZVU7HIrbZA21BEg9V6S01sQg1qNSpeoifEVgXphdNTRKsRLGKLbYOEIoKlUG6WKKkmgo2tUlHK5cDh/C773WvNOcbjH2M8Y4w55lz7fX/vbd++n2Rn7TXvc655+c7nKqoKQgghhBBC3E1vACGEEEIIuR1QGBJCCCGEEAAUhoQQQgghJENhSAghhBBCAFAYEkIIIYSQDIUhIYQQQggBQGFICLmDiMgXiMj/+xLzf6OI/LuvcptuKw9pXwkhLw+FISHkRhCRnxKRpyLyroj8rIh8k4g8fg3r+T0i8kPtMFX9ClX991/Dun5KRH7rc077/SLyz7/qbehp9/VlBTUh5P5DYUgIuUl+p6o+BvAbAPw9AP6tG94eQgh50FAYEkJuHFX9WQD/A5JABACIyF5Evk5E/h8R+bnsEj3fml9E/k0R+esi8o6I/ISI/CN5+N8J4BsB/KZsmfzlPPybROQP5P9/UkS+pFnWICIfE5G/N3//B0TkR0Tkl0XkfxWRL3iefTJLZd6Hj4vI/y0ivz2P+4MAPhfAR/N2fTQP/0wR+V4R+SUR+asi8rua5X2TiPxxEfnzeT//soj8bXmciMgfEZGfF5FPiMj/JiK/vt1XEXkE4L8H8Gl5ne+KyKeJyKWIfKhZz98nIr8gIuPz7Cch5H5BYUgIuXFE5DMA/HYAf60Z/DUAfi2SWPzbAXw6gH/vxCL+OpLQ+gCArwbwLSLyq1T1JwF8BYAfVdXHqvrBjXn/GwC/u/n+xQA+pqo/LiKfDuDPA/gDAD4ZwL8G4M+KyKc+5659FoC/CuBTAHwtgP9CRERV/x0APwjgK/N2fWUWbt8L4FsB/Iq8Td8gIr+uWd7vzvv3SUjH6g/m4V8E4POQjtcHAfzjAH6x3RBVfQ/pGP9MXudjVf0ZAN8P4Hc1k/7TAP6Uqk7PuY+EkHsEhSEh5Cb5dhF5B8BPA/h5AB8GkgUMwO8F8K+o6i+p6jsA/gMA/8TWQlT1v1XVn1HVqKp/GsD/CeDvf85t+FYAXyoiF/n7P5mHAUkkfZeqflde9vcC+J8B/I7nXPbfUNU/oaoBwDcD+FUAfuWJab8EwE+p6n+pqrOq/jiAPwvgy5tp/pyq/piqzgD+a1QL6wTgCYDPBCCq+pOq+jefcxu/Oe8nRMQjic8/+ZzzEkLuGRSGhJCb5MtU9QmAL0ASNZ+Sh38qgAsA/0t24f4ygO/Ow1eIyD8rIn+lmfbXN8u6FlX9awB+EsDvzOLwS1GF4a8G8I/ZcvOyPwdJ4D0PP9us5zL/eyrB5lcD+KxuXf8UgL9la3kALm1Zqvp9AD4K4I8D+DkR+c9E5K3n3Mb/DsDfJSK/BsAXAviEqv7Yc85LCLlnDDe9AYQQoqo/ICLfBODrAHwZgI8BeArg16nq/3fdvCLyqwH8CQC/BcllHETkrwAQW/xzbIK5kx2An8hiEUiWzD+pqr/3fe7S89Bv108D+AFV/cIXWpjqRwB8RER+BYBvA/CvA+jL1KyOhapeici3IYnQzwSthYQ8aGgxJITcFv4ogC8Ukd+gqhFJ7P2RLHQgIp8uIl+8Md8jJMHzC3m6fw7JYmj8HIDPEJHdNev+U0hxer8P1VoIAN+CZEn8YhHxInKWS758xgvuY8vPAfg1zffvBPBrReSfEZEx//3GnEBzLXm6z8oJI+8BuAIQTqzzQyLygW74fwXg9yBZS7/lBfaFEHJPoDAkhNwKVPUXkASKWbm+CinB4n8SkbcB/I8A/o6N+X4CwB8C8KNIwufvBvDDzSTfB+B/B/CzIvKxE+v+m3n+zwbwp5vhPw3gHwbwbyMJz59GssS9invn1wP48pyx/JEcR/lFSHGUP4PkNv4aAPvnWNZbSEL64wD+BlLiydf1E6nq/4FkHf2/srv60/LwHwYQAfy4qv7Uy+4YIeTuIqrP42UhhBBynxGR7wPwrar6n9/0thBCbg4KQ0IIeeCIyG9EKpXzt2bLJSHkgUJXMiGEPGBE5JuR3PS/n6KQEEKLISGEEEIIAUCLISGEEEIIyVAYEkIIIYQQAHe7wDV94IQQQggha+TZk2xzl4UhIeSB8NVf/dU3vQl3ng9/+MM3vQmEkDvAXU4+ubMbTghJUPDdPiggCbkXvLDFkMKQEPLaoPC7f1A4EnInoDAkhLx+KPTI+4VCkpAbgcKQEPLqoAAkrxsKRkJeKxSGhJAXgyKQ3BYoFgl5ZVAYEkKWUPCR+woFJCHPhMKQkIcOhSB5qFAoErKCwpCQhwaFICHbUCgSQmFIyL2Dwo+Q1wOFI3kAUBgSctehECTkZqBQJPcQCkNC7goUgITcDSgYyR2GwpCQ2wqFICH3AwpFcoegMCTkNkExSMj9hiKR3HIoDAm5CSgACSEtFIzklkBhSMibgmKQEPI8UCSSG4TCkJDXCcUgIeRloEgkbxgKQ0JeFRSBhJA3AcUieY1QGBLyMlAMEkJuEopE8oqhMCTkRaAgJITcJigQySuCwpCQZ0ERSAi5i1AskheAwpCQLSgGCSH3CYpE8pxQGBICUAgSQh4WFIrkBBSG5GFDQUgIechQIJIOCkPy8KAYJISQNRSJBBSG5CFAIUgIIe8fCsUHCYUhub9QEBJCyMtDgfigoDAk9w8KQkIIefVQID4IKAzJ3YdCkBBC3jwUivcSCkNyd6EgJISQm4cC8V5BYUjuDhSChBBy+6FQvNO8sDB0r3IrCHkWFIWEEHI34P36YUKLIXnt8OZCCCF3H1oQ7xR0JZPbBcUgIYTcXygSbz0UhuR2QEFICCEPBwrEWwtjDMnNQ1FICCEPC9737x+0GJIXhjcEQgghPbQi3gpoMSRvFopCQgghW/D5cLehxZC8L3jBE0IIeV5oPbwxmHxCXh8Ug4QQQl4WisQ3CoUhefVQEBJCCHnVUCC+ESgMyauBYpAQQsibgiLxtcHkE/LyUBQSQgh5k/C5c/ugxZDwwiSEEHLj0Hr4SqErmbx/KAgJIYTcNigQXwl0JZP3B0UhIYSQ2wifTzcLLYYPDF5whBBC7gq0Hr4wdCWT01AMEkIIuetQJL4v6Eom21AUEkIIuQ/wefZmoDC8x/AiIoQQcp/gc+31Q1fyPYMXDSGEkIcC3csnoSuZUBQSQgh5WPC59+qhxfAewAuDEELIQ4fWwwXMSn6IUBASQgghSygQAdCV/PCgKCSEEELW8Pn4ctBieMfgCU8IIYQ8Hw/YekiL4UOAopAQQgh5fvjcfP9QGN4ReHITQggh7x8+P98fdCXfcnhCE0IIIa+GB+RaZlbyfYOCkBBCCHk9PACByBhDQgghhBDyctBieIuglZAQQgh5s9xT6yEthncdikJCCCHkzcPn7xJaDG8YnpCEEELI7eAeWQ9pMbyLUBQSQgghtwc+lykMbwyefIQQQsjt46E/n+lKfsM89BOOEEIIuSvcYdcyXcl3AYpCQggh5O7wEJ/bFIaEEEIIIQQAXcmvnYf4tkEIIYTcR+6Qa5muZEIIIYQQ8nJQGL5GaC0khBBC7g8P4blOV/Jr4CGcOIQQQshD5pa7lelKJoQQQgghLwcthq8QWgoJIYSQh8UttRzSYnjTUBQSQgghD4/79vynMHwF3LeTghBCCCHPz33SARSGL8l9OhkIIYQQ8mLcFz1AYfgS3JeTgBBCCCEvz33QBRSGL8h9+PEJIYQQ8mq56/qAwvAFuOs/OiGEEEJeH3dZJ1AYvk/u8o9NCCGEkDfDXdULFIbvg7v6IxNCCCHkzXMXdQOF4XNyF39cQgghhNwsd00/UBg+B3ftRyWEEELI7eEu6QgKw2dwl35MQgghhNxO7oqeoDC8hrvyIxJCCCHk9nMXdAWFISGEEEIIAQCIqt70Nrwor23D74KiJ4QQQsjd5cMf/vDrXLy86Iy0GBJCCCGEEAAUhitoLSSEEELI6+a26g0Kw4bb+iMRQggh5P5xG3UHhWHmNv44hBBCCLnf3Db9QWGI2/ejEEIIIeThcJt0CIUhIYQQQggBQGF4q1Q6IYQQQh4mt0WPPGhheFt+BEIIIYSQ26BLHrQwJIQQQgghlQcrDG+DKieEEEIIablpffIgheFNH3RCCCGEkFPcpE55kMKQEEIIIYSseXDCkNZCQgghhNx2bkqvPChhSFFICCGEkLvCTeiWByMMKQoJIYQQctd40/rlwQhDQgghhBByPRSGhBBCCCEEwAMRhnQjE0IIIeSu8iZ1zL0XhhSFhBBCCLnrvCk9c++FISGEEEIIeT4oDAkhhBBCCIB7LgzpRiaEEELIfeFN6Jp7LQwJIYQQQsjzc2+FIa2FhBBCCLlvvG59cy+FIUUhIYQQQu4rr1Pn3EthSAghhBBC3j8UhoQQQgghBMA9FIZ0IxNCCCHkvvO69M69E4aEEEIIIeTFuFfCkNZCQgghhDwUXofuuVfCkBBCCCGEvDgUhoQQQgghBMA9EoZ0IxNCCCHkofGq9c+9EYaEEEIIIeTloDAkhBBCCCEAKAwJIYQQQkjmXghDxhcSQggh5KHyKnXQvRCGhBBCCCHk5aEwJIQQQgghAO6BMKQbmRBCCCEPnVelh+68MCSEEEIIIa8GCkNCCCGEEALgjgtDupEJIYQQQhKvQhfdaWFICCGEEEJeHRSGhBBCCCEEAIUhIYQQQgjJUBgSQgghhBAAFIaEEEIIISRDYUgIIYQQQgBQGBJCCCGEkAyFISGEEEIIAUBhSAghhBBCMhSGhBBCCCEEAIUhIYQQQgjJUBgSQgghhBAAFIaEEEIIISRDYUgIIYQQQgBQGBJCCCGEkAyFISGEEEIIAUBhSAghhBBCMhSGhBBCCCEEAIUhIYQQQgjJUBgSQgghhBAAFIaEEEIIISRDYUgIIYQQQgBQGBJCCCGEkAyFISGEEEIIAUBhSAghhBBCMhSGhBBCCCEEAIUhIYQQQgjJUBgSQgghhBAAd1wYfvjDH77pTSCEEEIIuRW8Cl10p4UhIYQQQgh5dVAYEkIIIYQQAPdAGNKdTAghhJCHzqvSQ3deGBJCCCGEkFcDhSEhhBBCCAFwT4Qh3cmEEEIIeai8Sh10L4QhIYQQQgh5eSgMCSGEEEIIAApDQgghhBCSuTfCkHGGhBBCCHlovGr9c2+EISGEEEIIeTkoDAkhhBBCCIB7JgzpTiaEEELIQ+F16J57JQwJIYQQQsiLc++EIa2GhBBCCLnvvC69c++EISGEEEIIeTEoDAkhhBBCCIB7KgzpTiaEEELIfeV16px7KQwBikNCCCGE3D9et765t8KQEEIIIYS8P+61MKTVkBBCCCH3hTeha+61MCSEEEIIIc8PhSEhhBBCCAHwAIQh3cmEEEIIueu8KT1z74UhQHFICCGEkLvLm9QxD0IYEkIIIYSQZ0NhSAghhBBCADwgYUh3MiGEEELuGm9avzwYYQhQHBJCCCHk7nATuuVBCUOA4pAQQgght5+b0isPThgSQgghhJBtHqQwpNWQEEIIIbeVm9QpD1IYAhSHhBBCCLl93LQ+ebDCkBBCCCGELBFVvelteFFeyYZ/6Q9+JY7RAwB2LgAAIqSMj5r+n+O2hlat07437TD6AFWBdxEAsPdzWcaQh03BQ0TLvLY+1+ySzW840ZPbcQgDnCjeOezhXUSIDtPsy3rbA+VFMQUP5yJ8XuZ+nMv4nQ+LZY/5u61/sd55SPt9HHGcB6gKQnAIIe9XSMdV865ocIBTaB6PqdkXn5c/1+OZFpK+S6jD3UHqOFFIEIjW/XTTchEuCBABNwPqAZmbcXlaZ8N0+V3ier8ldt9nW48tM80jmqaVoOkzLsen9WgeZgcpDz+GvN/NtMd5MUwOaePl6pgmmNJ4nfJOHSdgGKCXl5BxgF4dmh3Px1664x3SejV0O9mg0/HkuFeJjLu8wrwtPp1PkrfZtlG8A8YxDXcC2Hyu27fFsse8zHwc2t/ZhnX3Rh3T+V7OiRDS8QsRus/Ls+Oap407mydtqw5uMVx93UYbp912l/XZae/z+S55H6MCTqDSLM8usYu0zPksf56n4ce36jqmJ+lzvkjrCW+l82j3wXS+vPXoCgBwPk64GI94PB4QVXDmZ8zqyn0TAMb8/9MwYo4OHz9cpMOS7x9R1/fWkO9p50M6b+2ec+brhfxkPGCXL8oPjk/TeDfhsU/bduHSOTlKnleOZZr0vbspAHjk0v75fNG57sKOWu9P7+kOHoorTb/z2+EMAHCl6Vz7+PyoLOPnj2/h7TmN/4Wrx2UZ70xp2OWUlnF5SPNO+T4Z871unut65+MAcYp49JAhrp56avdHO6yi6T6bj3d7ecuQ9s+5NM4P6ViNY/5s7v2S5z8fZzjRk8+B82HCHB3O/IRjHPDWmH6PoTknzvPvOObj2z/v2mnn/Cx+GtIxmtRh7wLm/FvYObB3M2b1GCTAicIjYlKPTxov0WPrsXPBzpGrOJZpbJsm9WWaSf3imWzjbRoA8IgInX3NntE23Tv5XAhFS9T5W60BAFfziO/43I+u9uEFOX0DfAYP3mL4HZ/70SLYIgQRgmPwOOaLtb8QbNrBRQwulgvoOmwZ/bJMkDooHBQienJ5ThRRBU4UrhGVrTA1TFTa+ramAdKNOargOHsc57QtUz6pbf8BYJBYLhC7iZkonBqhatvuffdAjSfOT0W90bWCUKX+5cGab2buqhGFNnleXxmS/3FhvV4J/YAs6CSJQRcakYj6oC0P3I2fR4e62UC5Jy+ncc3fkP4kKtQBMsfyUNchi+pdPv5ZJEgWKUUozhHq/XIlvcgDgNnE4gxVhaoiHtIDMT59inh5ifj0CvFpuqHbi6KaQMwiUKdj+Xut2D44X7ZDGyEm7T5uCb8NIW/HACGkv6jQwzEJv9nUfF5WLwqjlmVKFt69YDyJ/WaqgAh09Gn/nIObs6jNv6eOS1EoUctfWVyz7+WS3tpfu0QGgT8q1Av8pIhD8yKUP+3ZGPdZuO0j5ODgHs2Yp/yQPo44H9MMj8d07pz55iLJtKIQSA/AJ+Nh86Vy6N+uADydR3gXi2CMKti5kP+WorDlLL+ZmfgzUbiTUMRd7B51b7l0vpsoHGWGR8QOATsEeKTtO6rHUX2ZbrHe7g3UhOUoAR8a3wMAfGh/mccthXFUwdlugnMR+3GCiML7WMSh3f/Enf597SRoReF62nzfHJei0PkIVcEwbL8Ajj5i9Nvj7FlgQr4V8LMuj/M+nxO9KDSRZMfMnkVPw1jOnx4TmPt8Lgyrmzkw5WdqUIcLd8SFO+LMTavf6jJmUZ6Fnj2D2/1bLVv95vCWqLIQj0/DbrG/7XHoRaGDvkpR+FI8eGFoDK6KH7uIe3E4dILLkPxG9Wg8YudC+m7Wx3IB1GWbtXBwsYjDfnntvPGEsDPLpL2xP9nXm/Dgt8WhnaBxw/I4hfymk8cVK2W+2Ftx2IrP3RDqG2Ue3ovDmB8yetw45UwU2rj+umzEoQ5pvP0E0onO/ki5KR//tcEA3jROXtaJw7wSh5vTmZ6Zq7Uwbd/2MreWbw/5U+JwExM2Zi085J06HE7MAEAV8enTlcAxcRivskhsROEbpYjDZ7/wilnxtkRxy5aAOjXuGuFXrLM2TW9ZzZZB+z1kWj+86sLyp1+K/+uI/sR+5n1wx2yZvOZ4lGe3fV7zcmvWpFM8GdJ5Zve7QxhOTmuWFBOFg9u+OM78hJ2bcTFsXLS2XcUyOC++mxWwZdcJiF4UmqDsLYdbVsYLScv/oL/Mn0kAmuWyxUTpligEgDlsC43+3gkAbpfvr8PSetz/zOLSzdGNEeIUwy6sRKG47rs97+b1b2fPwPIc60ShMbiIi+GIQSIGF+BFV6LQMPFmwujU8w1IwtJJXInCxz5bs4creETs3YzH/oALf0ToxKmdG/YZOutgi4PCS1xZC23aQxzKNrfWwjm61X6YKDxlMW29kQ6Kd+f9yePwpnnwrmTjS37wX8QxeJz5uVrysnWu/8Hb73bR79yMYz5pzFVs4/v5VaWMtxupE8U+v4WH6Ip4bD+33rzNihhVMEWPKXgEFUzBQ7M1UJGEoHMRMS87BFe3sXkr3Q0z5uBxsT/Cixar6JYryKyHQQWHKe2HfaoK5tkhzB46u/QcVVlaCc0NMjmoaBJ5jQiUWaCjlnlcFo4yJ0Gosn6m2TPA3Lv+IMkqF9M4iUjWwcZyom5pSbT7mMRkbWm3WaKmbWyEX3UR53UuXMX5H1W42SxBtq1pWWJi0lzEWsebSzntU4SEmMRGBCSLQLky69dcXcAm6o5po7aE4K2kedLJbleHh5BcyVGT6xhYuMPF++K+NdeqWR0XxGQxxDBUa+xul46NCDBNwH5XBV+2nELzPECznljXN/iyLeqkbJvuzGUtZd90zA9GsxKODggKeKluZE0vDG0oQb026guKZrcy8nmlDoAThJ2D+jQ+DoKwA+a9FAvhlF3JYUzWwrgDogfiRYSOEe58hs+i8JPfusSjXTqfPrDLblw/lwde6xp8GsbidZnVF8/CqXAdu38B6QG5cwGPxgOG/MK7dzMiBE+GKzhRXLhjmf5Xjm8DAC7cAQEOH/Lv4krHhaDzEvFBd4n3dIcnRRBGTDrASbIQtowSECDYIeJKB1y4Cb8ULnAmE96O5j4e4fNF/ItzchW/E5OP/jLu8LOHDwBIYuljh8dwovj48RxRBe9MZ4gqePvqrJw6hoXjmNVwOmYREgXONaFHQeC8Is7ZolgsiM2OCDCMARoFwxgQo6Soh5CWJQIMQ3UPh+iwH5vjlp8J9rs/Go8LI0f7LDKLoRNtXL1ZjLnQGFxiEUv2IjFIwCEORTRNWdR5UQSV8jma4MwiMMBhlFC26bE/IEIwSigCsMdlwTfFKoB9WW8VgiYCW+FognBSX5bvRMvLTkAVhu0+mlXThs3qcAgDZnU4hgGDC7gKI5wovvNz/9jmdr8EdCW/LPajXHVvvNe9zQDp7dfegO2iGLs4vfZCus6tG6Ir8Ta9q/jUfO02mlCz2EETfsXrtLHsfrn21uibm7WtYwq+bF/PfpyxH2cMPhZrpbkpLLalKh4AswOiQHKc4cryV0Sj1HkazLWsUsXZlihMy8rz2HXei8nrrHq9McnV7UlC75p5kR+2eb2xeyEvFsihWiRbUbi5OUUUdSOy6BDvkutU3FIU3jXEQafGXWnxhf7ELWvsDu6WhTDmE2RYTqvHzgp4OGEhnec0zXFaisLFwswFHZNgneZiETyJuZP72MLmHFAnUJGFKASaFyOb18b75TLMOu4mIOwF7pBEIQBkrxriRb5e8zUZJo/HF4ciDnqe5gdaG2pigs9ilfdDjld0AYOLCxdaf2810WiWRLMOlWU3FrcPDJdZuCXL0Yf8uwDWVr4PumTZM1FojPkmcYQvYsIe+LvsRr7o4hMfuQMeuUMRE1tMccCHxvdW3qVPyoL6yVi3o48jP3WPd5072WWLYrEeti/bmSGLevssRvi8LBOFW9uyJQoX6zevlIQS47dltACS0WBSt4rfNKuhCa5DHHCIA7xoefbYp4lGc++bla4/f+w82RKHtv6puwm3ruRT2DS9hbEVhWl7HKI6jC4U4dtjhqBj/rzKF+FrEIUvBS2GHb/tB/7l8v/OBUQIHLR89jEUxUXcuSGOjcm5nfYQBowuFAtfK7T2XdxOa3VsidkNbXEZZlG8mmtsxtW8tOBNndvCApztTdQMJhbnuEhIGZY37FMxk1N0RVheTQPm2WOePEITTF1cySYIzY0cl086d5SUKBKToJJQH3ZmmStistedJhBt+CzrJJNM9kjUhBTbjO6B3wrCYsHs3MZ2RpaH8KzL4ZaUYprCA/6Y93fSlKSito/5WB/tM9/gLeHkkF2VV4f0wx2O0DlAzQ389Ckg+UXgTbuC3y/21MrbK71AGga09ykxkbgbawKKib3yO5l5zQ5oPnk0pmGuE5oidVqL3TSLpCXz2Dps2dOUprXpvEv/q1YBnwVrSU7J9FZDIIt+l9zA5lbuxeJiGd249jmpg0McBHEUhF3+2ydRCAA5VwLHt/IL4KMsvnLiyfkH0nm0H2ec7SY83h0xuoC3dlcYJFltLGnAHuAWG2af707JNfZ0Xu57ayUcJIlFE1sWu2ixjJbM8Km7d9L25Av1U8Z3sJO5xhU2Is5chY8sAaUTi3aPNgFoHpdRIgKkiTuMuLIEgqyeLfnkvZj27RfmtwBUcXCII650wCGO+Ph0gUk9fnlK1sS3j+nzE8dkLXzvmJZp92h7FkyTJVkpYnBwOfbQufS9P5Zb+PJyXoWhar3fexfrfTzf18/yPd+WaHGlrcg58xMuhgnvTnuc5cSTZN01l/7yGWHi3iyBx0aYOTResmwddNAiFs/zjXTfJJy0WMIRkCx/JgbNqtda90rCSQ4Ivwz7xfCWy7jDpL4ktNj//XPPzvMax5qOXBv7OKsvyaGH6HEVRhzDUKado8N3f/7Xr7bhFUGL4aviuz//6zdjX0oGVbYQ9sHT/RuRTbP1JmVvS6NLsXm9hdE3CS7t/Fs3gdYS2SavlGVYrKGL2ZWcLUsbp4w9F5MLeigi7+kxXQCheWNP+7GMSWyxm8wwhu2kGmlEIaqYc8eNfZyX01QTaHYPa+dSNm2wkXzSW+3Kd5divCQq3FzdvQsX8mpZS+uNWRFt+vIO0bkAbb7FssblstRJWWcRhU+nIoKfJ/4OGm+/KHQ+CUIThd6X/9En15xA+phE++59FXMbx+uk9RFILuQpWwbb5Zxia1lSt0eO9aGmXorFsWSjA+Vu/DyxhlALNYj15aGxMG6dY0ZOqEXYJ6t63KVY3HgRSryvJZ4AwOMTFkOgCoEp+nJPNBezibv2hfc6D0yN5Q64CgP2bi4WlsvsivMS8SljEolbMYBluxAW01y4qbEAhoUobPdjCxM15r5sXckfGt4tovSQ/fSHuE6eMPHQGhZP6/bMAAAgAElEQVQe7Y5QFeyGUKwc9tMXo8BGAsgqISUnqpjlzuJCW1EIrC2PadhSFO7y86h/Bl4Mx/yXj+FG/OfWMTRLILAUhQAWorDnkcWudtZc+yxZ6G5Of40oPIUJdvt97HdrXcetZXBrWVFlsU8lxCrf5O18tWnq9+14xtcoCl8KCsMTmNvDQcub8Zb4M+KGudxo4/PshmIBsnajvE4ctpiFcG7czr1lsReHqxjJRhzOc4o1tGGxc+keO8tj61oGluIwqmA3zCczq9Ush/bZa8VGHLpZ1hnEAPzTHOdiorKzEtaF2fC8XxuJbjZPEWPSLfOEONwSiUb0TXA/cFIc2vBrLUKSspRX902zhJrL0+Lh5s4cCmy/AdwBLL5Q9unNvlgJO7FYxN01Iq8sq7jhTzw8bNlh2w20wKyIrltvFnzXnSOGhRKcmtbOx+c579xhe5vLs3gVinF6u8aLHOg/Ls8niyXr7ydb2aFl/TlEZms40CS2lXvY+h76KCcVXLjjZpKHuZPNPdyLwnadF25aWARX68rLeOJMaCzdhUZQhwhX7uNbSS+Ps7jZEoXA2ouzy2733s0LtIkjG6IQVeCZkNzt8rI2hKV3Eftx3hSFQH3m2LNqd8Itmp6NKTHklLC2fX7vRGLF1AgmM760onCx3WIu/ioKgZo0FJqTvBd1Zik0tkQhUF3bc5OAadbClev6GlEI1N/dzqHLeVdcyNe5rm8Dz3gNfpiYv/9Lf/ArESGYo8eQ3b9OIq7CuLpYjtFj59YXfxsXU2MxtsVha6o2N7HLCSBAWx/JwWkSgCE6jD4kkz4iZriU/IKUUSaiGHwoGXDmYgjBXMnZrayCGAGNDmhuTFpu4MjT1xvAlqvb9kVEyw1u9lbPMEKbG6wEqaVmgqQ6gzlmEK66k9NK8jwAMEmymMjSxWzJJWkZadkSpQi7IsZsmZLd0g5AEYkCgabgf5emkYjyXUIaJ9CyzP6+2Qs5dWl+c4mLoiTOpG0RuFkRR7e0IAGQkmUKODMnzCHFrjmXhYhmi9YhuV0txnA63mphKI0VTqMm0acRcn6eStOYUPMeUE0CL4QkFjWWZBSU31NqLKJz0BiLNVGyi1eGARoixLv0ackhJga93/6/tUK2ruUYa+JJa+kTgYQIHXytT2jJJCKLYWX5loAiApkjMEhKTMrnD4D63XAoMYoW3hB3DhIUIdcuDLu0DnMjG3FMSSfhLKZj6BWyi5gnj/NHB4z5+t37ubgK+9pz6eW5EjTdL4/Bl/tdQHUZA2thaa7Ki+FYwmQGiavQGrtfekSMMmdLkrkwbftiEWmjhCIa+nU+cROu1OOsER8eilFyAoqmfbmQGQf1OJOpuJFtO1zOSvWIeDueY+8mvBvOVqLk0XDAJ6ZzXAxHXM477P2Md49LoXQ2zpiDw36cMQWPOb/49q7gEFJCn+REQj+EUg7M+QjvY04uUQxeoaj3aXuO2O96Ps4psaMzQjgozoapuPUvhuNKuLt8PHduLqFWQDWWtM/BFIsaF54yIIlC11iZ5+hx7muSi7mOx5LdHErWcD0XNMUISsQUh0X8Z+tKNithEvSppExQ14zP7ny4lCXcJJ7aZ81A9pjUYZSISV2qDRx9YwCqLwLlemgFsER8++d8A24ztBg+A7sRzp0p+Jh/7Lb+0VbpmTahpBeA1wUxPwu7EPuL3i7MVf1FH9LfdfUXraTN7HOh6vxgyZ+l1mH+PjelDNobr3cR3sXiJvdDKMVUeyQI3KERi7ZZJXM3/bXhJSWWcKrTLJZZ/DLmQ8lffV6OYjtpxKU/qzfYsrKw5GWsDmOeLhTXcBWmZTstfsysiI1QXQ7Py/AO7jAnMTF3O5sFzyJrFk1c4S2OIdbGMlcsgue5ArNvxFYINdnDhpu72Um1FjaiEGhczNU/t1h/EYWGracRmJu0ruUtUQgAMdY4w4ZiJZxjM8y2P4/rC2vbdtss/bnrJWU2CxB3+Zq1gtblM8WxmhvZilmXMK3mtBr3Mx6dLd3Hx5jipI7RY+9TBYYaW7hbTHvVBdYbbTkwC5NpRSFQrVaPhgOCCh7lIOB9V6twaixAOwk4kwlnMq0sdwH1hdzuiWfFLbk8kPmSRZ/3tS+lcdbu014ETuoXsW8mOD6wUX8RwMnatbvdvHIF233Yys34zhpo383r0y519AHjEIooLMlB3tahpWwKsF2nEkBJDhlcWBSmLs9JdSu36WlP2jL0oMXODye6cCXbM3MssYtWf3edVAKkItZXcSzP6b52IJB+owBXhN+p8kL2nJt6y2+3v7O6Iox73eAkLgqn31Zu/xbeIFZs0i4CyzraysI6VZPLLvo+46wXci2tdbEVXX2iizGfcClXgXZ620rciW1nHqeda9kymkuB6y5g+nnwY0wxTJZRXGKquu0Ky88y/7FaBZ+lqcs0/T23F19u+b2fXwJKsskzrf+W/ZyXGbsuFKfcxxJSTJjMuso6NczCZCLIunBYvbxi2bIi0KfcpbeFPuHExNHUPSS8B8YxWwfjdlyg9CrerjmtywCK21eG7tg0sYDtdNdutx1vE+rm5rV5T3WWaZYjQZPos3nt/JGleGxFZFo2ltNb7Oq43m4/KcIunWPhrEnCOpoVMS9jZzFt6fMqX99P9kloDRIxR18Ew2aB6jBuTlMshSdi4fuktvPcRWNsypv0PPFPceaORShU9/CyvqEN7wtU25Hy0GQphK4EoXHoxJ/VQTRx2iaeAFjV0bOizakczwxVWWX6GtobDzqxByS3cl+D8BRWcqzFwpbss+3S1X7f+Rk7P5/MON7CRFItNG2iKh0TswL2ySqjBJz7Y05q2r65967lWidyOX216C5dve34oA6T+iLc6zKrCNyKNyxlZ7LV0Lahz0IOKggqJQu/fY7fliLW10FX8jP4js/9KP7RH/59+S03pebH5k24jZEZXCxmaMvWAlBqEWpTeM/EnHcxFco8IbDatxf7PyCd9JZwYuu1dXkXgehqyRmVIg4nINc01FzbSmsCJxSiyUmqKnARQI4/TM/YxrIHYA7JVH42zuWG1sY3Dj7iahowDBHHw5DWtcuu0aMDvEKRXb4hu4cbUeiCIELhZkEctMYV2nWuQByam5a5aKOsaxNml3PzEyTxJWm3zLAnIQ9r1gEkd3HMFsc4VGG6ej7acK3LKO5poNRUVJfcx5rXDQU0i8OaXa3QnYMcI+JugJ8O0MGlxYacXTt4wFrdTcdkQYuaWuBNqCVabhHStqyLWsXhblxYB83VrGr7Y51HkhvZ2tppiKctfDbc3L52sg8DSio+sC3g2um3xGIIaXjI01mdwsGneMDnFYkASu3FiOzWza5l654y53WUF5l0TtWYVRN62TsxNG5qJFEIAPNF/n7eXjfpBPTnNWnhYtd0s4gulaFqXICDj3g7t3gzjjkBZdZ073s3t36zl8k2jjqqFCEySCj3t10nHCze7MxN2LsJo8w4a+IELfO4LUfjG1Fo1r7dhpAdochRK/U4afoJAgQ7iXgnppI80CREdgh4O9dKnDRZmpLASPt64Y74WK5t+In5HHs3Lzp5WIMD+/9inEqLvFYUKpA7obhV3DdgMYdZIA4pqdD7iBAcvF+3O40quBgnTNHh8XhMv1UXS2iVLnYuJFHYxdZHdSVruLV82ZkSs6VsyOsLKjjGAYPEUs9wKG5hW2dYWQTNVWzDtzxrMQu7vZsW1sJWFE7qF/MedEBQtxSDUr15Jv5WMflINTmn6HGIHqNIPgYp0/jxcMTb0xlGF3BoahMCTVH3fGz/3D/4n6z25TZCi+FzYD9mcUU0BT2NtoXcVn9lczOfEoA+u3dPlag5VSLGbibtW9YgyVJ46s3LZTN9fSba22f+3u1HdSnX720vZrMutBxnX290Lq6y60wglh7JJyyA1tbO9T2UbfwsqbuJYrMFHoAq7sz13L0At1nF7fTrbbF1nl6HZUivXNGdKLAM0rIsE5GWTVosh+t90ussWs28t0oUmhVs3wWhO1laATcEnglEOdvnmMLtBBTpjov0mcQ2vmQqbxxHGxZPnJArl3Y+1l0rPStPUyy9ParpJWWjJ3WpX9jXP+zPKcuK70Rhu5wwSmnZ2GfkA0A4y9f+WUCcHPadC7lYUDqrS9+2rNwLrul88qyasCZEelH4Ab90w5qlzkRhT4oVjNdmGhv9rzP1bcq6uLGyjq47yhO/dhWfSlYwHu+2uxP1t5/eZdxj8dxyzbQXpaVhOma9KKyNGsLis6WN+fP9TRTVSGL7bVnIYxGYZtWtojB972oO5jjCU4Wqr+KISX363DqpsXbxH3SZGFK2OaaYwlPZzEnsejwNY3GR19qL6fvbuWj5oQuhaEUhgDsjCgFaDJ+bP/PZ34gv/5GvKBZBo3Xzmjhs3cYlKLdpf2fOkd4N0LpUzOLYxkW0GcyW7df2Fi0WRRUMucr74CKC1kLWVrw0RkXMZq1kPUzbYAKuuJKDhypKr9He06cqGHPgdLtPqgK0RVOtbM4YoXlwHGPqeuIADArM1sWhO/iaBFRxPVuiiG+/J2tKKw6L5S3Pa0kf7ThItSaWtnlSb85FKwoWSSh9QkubwVytgnUXJNThZaGNELXyNm5WRC/5u8CFWEXifsxlT6olKR3UZYydvndZizbfJpqyMrLbQY/HptRMjpW0OEGRZCks9QB3SehawklOspFhSG7ZLskkF21brr9Z/spSaBZEG+ZcEnt2/OdQk33aRJRuOWq1DPN+llVbcgkAHd06/tDV33zrmFniSiqCrk14QhWFKYbQOp1IOrc1XeP2PIweKeFkh2ypB+KTAEwO4+MjpuOARxcH7HxoXIo50URlkVRQXMWdALI6bU500T7Tpjcvh7VPA2pMm/VCTpmuocTr7fPF2QqS93SHD7n3Futuk01s2i1rYVpW/X/SmtVqn0d1GBFLRw0TDyamHrlDqWdorkn7HwA+abzELx4f44PjU7wzV+tq6wUyS957xx3241xCdAzrOhVFsduF8pJu99PiqUGqN2v34dJkoKlMcT5M5beQ7plhzxonsfzeNl97Sh5iqsObklZy4kaTaTu4UJI17fjXuPelIOzrEi66pGSLX0uE1OSROC6tfE33nTLdNS8ilkxSytzY8y5XF7FuLO91vY5tO+27xRDOWnVBafnoY3Ej/5nP/saT23IbocXwfWA/7qkTri8X01Oyi/ON47puJu16LHB2q35iyGViTm1TmyXcIt229jWuaneU/D1KEosq6wA41Dcv+wNqXKJhNzNxSL2TR3tQ5k9zC7dxf83ulri/zgJYxs3Lvy1WP02/K64OL8KyjWnsp++tj277c7HK7Eo+WarGjFaWwDIkl2sscYZ5uMUZitTC1leHtfi5DajC7fcQkVo+ZrdL8YPA5rYWUajNSbBVT7AXya1gWyxQchFqSS74waf/u/qHuh+rBbB0FGl+yOt6KTcWQG0sm31CiYRYh3XuYADpvLcYwq5Xd/qSYlhVqqXQ3MUtlpEMAHPO62lyI0q3E3+RiwifTbWjQ3Q4tyLGq/guwTF6zPlBXa1F6xM+tQarohBYelPa6QyzSF7GHZ7kEjWWWNImgfxiTJW6J/WLciXAOha7Dn9xap/dbA2TGZ+cu660430nhgDgyZh62T8eD6uXeWPwEfshlCSPrfi+/TiX8bthxm6YMTRizkThfpgX8fDSZNm2664GiSoij9GffKaYKLI4uus6n9hybdmlqwlqbCGwjtErsaMn4lRbS3ZrzT4Vx2p1Bbc6mLwX9iV/wKYFsEoy6fe/F4X9c9i25a6JQoDC8H1jaeZ9xpLR1zZa1fxqxGHEdo2vren6dTiLJUQVcX1z8r424ouIw1PaIgafeiEHwZS7nLRs3fjTOkwcLsWgXaulu0m2/J2K5WvFocQT7t12vVteicYsqILSus7ct8/Q7UUk2rataxba8nLM14nEkjStXJ9QY8s84Zq0eDvJYlGGjcKNN4wlgxTXr7mITSja91M1CS2h5lSxaROEGxa7Qhv7t7WNZ2lb1BJU+mnbLidA8+Zk57Gp+hMPytIG8ZptBGAZ8pvb2FyUOrjVeWNJT5YZ34pDAJg+kO9R1gJvn2MLc2HkR/ukHM9zEeOTdeyuEXmGiK6SMYDTyXqGZaqaKPSdG+Esq9va3WRebFOZzgTatWu7fvykLmc35zjtknSSjutlthpaPcVPHpIV0/Z7zje1vrWfYXHhxZrX7YPVHRx9WNWWrNNomQ6orQgN+y1L0uNGiZrF8BPFqk3Abb0AmJVwOOFqtvU8SxTGTlzb90/kN5t3c01EE3BXOb6zTXxpy8sYbSs96238XliGtrwXdkXw9eES7X4Xd/GJknM2/Nt+03+6eRxuO3QlvwDf/jnfgC/7oX9h1fZOm7cY+zS3sL1RWMKIuVlGHzCHFKAbmzevNivZlpWCtFNm2zEOKRbHp8KZqqkYdysOo6bWR2HDouid5ph5zXfFmoySviHvU3Uri9PiykjDfOow1jwEYxQMQ0jWwiEdk0XtQ6dwiCnuaXbQXGsrFV4UaADgAZ2XFrvVp7m+LSl0qDf3Nmawnafcy5rlmgHUhGX0yQNubclKgootsr3fJS9d+uhdy6jfW2NLHCTFF4bqPi7b4ps5HeCmiOgdfC5e7HL3DB1cEoDHKSU6HHNrtlzgWsYBepwgw3hrOp9IF9dXRWEWtLlmYElAMcHUiEU1t25zvhWx2YvFtuNJkxiyVe5Hx6GODwFWNE7HATKH9GkJJtOc2txZ4orqyuJY4ik1JVcBWcyZxa+zEJZPq3PYHzsb5qoobDub2P9hVx+EJgbnbEUMZ7XtHQCEJ/k4ZFHoh7AoUdMX2G/rqQI1Ge5qHjG4kFqk+WnhUpujW4nCRSx07p88SCxCcJTadcNEYYDgQgKCOjxxVytRaNt3hhrL5qFFFB7V4UxSMIwXYERK4LMteSd6XEjACMUEKe3ZUr07h2MnMEaZcaUjJh1wFUccsxXqSkd4KK7iuBCFDopH/oinYVxZlex+ez5OmILHfpgxR1eT+FprX5NgaPdUCw1qjQDmqhfRIuqtBV5bmy/1uc7PFF+FdXIprws6mws1Ntbh1N51aRndEoXFrezWXVWAKqInrfc/w7reTOrxi8fHZfjgIq7MpRx9SgTK4VwpMWTAKKEkitSi6h4DwqojyxR96rYTh5JxPEdpKpKYZbxaCtNxkJIAGjXVLbwLtQqvgxbDF8R+9K2TfMu127t0S+eQsD7JtkrUGCkbTEqAtlGsftntbBd3nS/XhGpa8O262oKpMOppD6Rasei82GhdRfKntdEy66F99i5lsxbKsExAsZIZduhOWuxOGGOs5qG5m1dWxsYlfE2zhoW1pl/GqW3qh68SWvrpzTrZu559FgKd1bH9UeQQNpNSAKTYu6tDElK3TRTmEi96PNHKrC1obX/tchr3rni3FoNAPU6DTyK0P07Wu9inOL/iju+3oY0ZbJdr05s4LduUH9q9pdCG967kPN2q37GXelfeMKqJarEuq5dNt/Z8lsbP5/m6vFiOD+cRMkkRheO+sy7lRAWzOllMmFkG7dN6s1/OuY9wtqIc5gHH4FdlQ1ouhmO5d/Z17ObocOGOCOpKbGHLVdxhRMAxF57eYpSICYKrTpi2NiCHJApbzE1Yevt2j0grcm1Ww6PVy7Pi18XVnMRiv+8fGNedW7bYDQE+F6o+5enZDdstVfelRuG8sESuXMiWse3nzfE9feu6mlSydBOXfXCpMPop62Fd7lqGTNHXvw33L7B8zvbWwUMcYO3rbFqgWm8PccAx/9n6bHia3hJHwmL+rW0A2msjLecui0KAwvCl+K7P+wiAKrpKsWmr/t601XGiK7dxf7L1nVCMtoj2Fm3R6vYmYr2aRxdLdfse/wyXjkgScqtWTJmapGKfFl8o3Wd92wUacThmd3gWh9HE4bBcX0kaKb2I84hsPVwVuW4EoLTmgcVCk8XRpfo/y/lNnHUWw01Ojc/btuq4Yqs3t2IRh/lYldhCB3+IjUUpT2huyHFphSv1/Myt/Kz+vm8I7VrMyTAkkeg3hBtQLXzWR7kfnQWftqKoXYcTWIHp4hLu1yOSjl/fUq+PKcyfduzLp2Udm9Ds4wTtsy9W3buON4Rdqm3YDTQPta1/SIlalpkcrbvJDnBBMZ9viFQA4SKVhtKLAGyUQdkPqfOGxRYuNqGxNrWFe1tXc9/mDViHtvRJK5bdee6zC7sTiuZytBjD/rOvWTh2B6/EtUkqrWJbe8iH5yIvPyIlqkyWaNfFTV7GfbIIapeRXaavNe/aIteG1TMEakFvw+7hp1zs5nkxS2Ffm1aAxX1+9EtB04u9k3Hynbu7Ly/ji3hM300s9okmhs137o44d8dVwkm7PVvu38tY6wbaOde+LJhoLL2Jmz7GabmuZBbP6ktcZIu5ntdhYXVf2qLVvRhMMZlVZN6FOoXP4nY8Oe4w3/V5H8Hv+Ev/0qKOoOasYEiOEWxOuLYLirmZ7aQcfcCMlI2sQRbunJSlnC/WXPcLQKn9NWvTkievSy3Q2KU6iaOLpbdxiILBp6K1ow/VBG+tlnxMCSdpr9ICTRya1VAAqGR3sjRrriYPVcEwxFzqxi1c6dEB8BHiI+LVkNaVsyZT1pzWrLuQxKG4upo2w7f0WZ6BOKZxcWwEo6QcBheS23nzJdal/UrXf/dgl84dXXzFzVdttqm3NDYuZftfBU2MmS621WoyWicLd6yRpuoFLrfKwzjUItd5e2U3Il4dAHG1Rd5Nogq1Fncxpmzk3S65vnNWcTllLI5Qm7Z1w5AylEUWJWlK5vIwLGMBzcVrq9+N5buEmMScZS3HCIxDEmPOQUJIxzTEZHW0GolOksicA3SfWvMVcXi2T/vokJbnXfred17JGcKwmoUlpGHpKq5JKUBbaK+ITxOD9jmmZdnzsLiPz208oANKv/CY6xX6s2wx2s21/imqxbAvPG1chXFRLaHvnz6rA3SdhWxuSmsvWlq0SSzFjfduLn1sL3wt5/LIHXDhDqXLCQA8yeIrqpQSNU5rNvIIxSjJPWwuZKDW3QuQIlwiksvZlhchi/66l3FfrINBHRwiLuMeOwm40mW5kza+3LJ0+zIpZb/GY8nq3UoU9E4xh+RabsVgecHO351oqQsZIdnF27qMlyLnzE+L38DqF7bbYbGCJhZ9Xk6QdIxDc6xt+aNLmdGP/DKbvJTH6foOm7t3ys+v3v3bxvud+ynF5ef4QIcm1CpnE8/qMUhInXqgRdxbpvQcfc5aFkRtO5S4khnvcjF3J1iFi/WlePoak9/5eX9s83e+a9Bi+Ar4rs/7yDOysvKF1b0N9idbcSs3WYHt9/LZ3oSaLLP+Ri6ydim3zCEVS15ZBroEkb7OYbGilUDEvB9zvrEGQZhdsR7OebjVQVQVhODW22UuZa+Ig64sbGVVDtvZvrlLiZvq32J8FrTXJar0ZWYW677mainu7ywuT02jTZzYeoJrxrlkQVwlQzSicNXiDVjEf94oG5Y/AJsWMwC1bM1Qy9dsLnY8kWTTdTspqysWv92J8bnkj1kC+2NqFkiz1p5yRZ9IcOlP+fX3RhS2n2WCvJ1Wpma3nOBUX2SztMd9vrbH5b2or4faX5vFfZwtXqeyP0/R99rtsXvaliXrkTuU+dyG6X9RPqzbrgmCMznd1SR1p8Aqo9k4wuM93RVR2HOlqeXa2UYZhEn9wp3ZY2L5FH1yYDvcspJHH1bLaOsGWgxny5mfFjGjW8f83E8rF/DQWA0HiZuiEAAe+SP2uSC57SdQE4haUQhUS2vv/u2TQCJSrcD+uVl7EyeL4Fwst0v3rnUrKVnFunyuGjb9sWxH/+KTx3e1i82DeB+gMHxFfPfnf335v72Q7Y3NbrxbpmybrsVM16ddATW9HlhmlQ0u4szPi36X5mY2TTecKIRq7ZZWNyXRE2a2JXHqYnaaOMMYq1hc4TX/dcvziuiTSDxxb06bd8IjXuIO+5hCS2DR5q9Bnf3JdmmZjXkWo4f0dzLG0MTtxngVqW7mrfmdO1042XnA+U2ReNNo1JRAshqh9a1DN8Y/q+yOSOMCbqbtrYhAFXY9G0KsLTejWyWAWlHYDt/4f1GUfFHfMH+eyk4u09XxbYKKWRvDXoqlsMVK1MyPFOHxOrB2t08P77NxxtmwFjeDhM2s41MZtnU+S45Y1sUr45v4wv1GX97H/gpP3NWio8li/RIxqbu2pEhfrBoArvKhm3QdMwcAE1KySfpbnytTsnkV93FPVFf68y6HS2mt+qzyLib6WkYfF/GCPaMPm/GILW3ZmBZ7Zmy5U3duXvVF7sebpffUclvMVWyi2WoRrgpSN27hVUmazj2clld/j0MYcAjDwv3czt8Ob7c1LVvKs7ctvTZHX9Z3DOnPjhmwfP7fB+hKfoV8zxf8UQDAF33/7y8uZWP06QY7nrgAzfIXrIB2vsA1C8koy6zAtk1RyQrbKG1zPkyLmlR9H+UQHUKUYq20WEBruRTgUusta8xusU55+WodSZr9iSUJxdX9sMxmUcS8rhBcGV4KRmdxqDG1wzNLUV98WpptAVATTnKxasQ8T3OfjKjj+7qEbcKLFbJOK7YZkYP+daGPS+xjtjQWi6PaPEvRqi7FgC0srlviMGcoqzq4q5Rs4qLWzFXnoAMgJrS9A6apul+9h7oZ4gbo4RZ0QIkBbf/m4k6OXZYvsBR3VpQ6hOR2tjZ0wNoq1/eRNsy1nEWhWQFL0ejB5fHZMrD3kKkes2I9NCumbZe5hIfkIi9uYsOhurBjrOLQ5pO0bSt9EuqxsP7HbfFrlfqS1HY7kQhMj7JVo0k4Kd1Njg76JCeU7OqDfHAR5+NUwllaL0Ta3PSgPMbtmDULiUmHWjFKWLTrBGqig3235DkTkE+GJAA/MKQC122LNOMqt6I7qsdbYlnMOdYNgpDFhoNigoPP0YRX2dV8qcAubpoAACAASURBVMCVepwhVLcxBFc5SWRCfRFvhcqYM6In1IQaK2ydunDUhIWgrljGLtwRl3FXROEUXHZ5mvW0umkdqisfjRDf+YBZ3EJcC2oM4a4ESS9/k61nxankEtuGNrazZuK6pkRN3qb8252vWuZ5nLuaVGQhAaWMTOlIUKt0tJzqpmPbN0WPcz+tYj/t+D7tQhymRvxZ0Wpg2cbWjpWDlgYUR3PtQzYbTNi8UeXeCUKDFsPXyKypt2T7xn2q5V1PiUW0cheN+Tw25vJWIALNm3ynNva5GXqfoWbf+4bqVtOwtFjqbyYWb2ifNr7frSCIk4Nmd7K25W4UJ0xigA6arG72LM1FsKNHdWWvNyvFAobskutejt10jQUxb/tCFJYZ67JxYpNbkdqXxVnVN+x6DxYLULdccyua9dAKXNv3WsYkr8CSTXIrPCluzWtMrW+YPhFlxamex7ZvpxJqTlkVS8LO9nJ19IvPOp8sP/Ni+uzi1XpNnPYdUbrl9yEDJc7whOW7LWmzhdUs7Mmx+4ijQs82LHdDwFnui9xntwL1vrJFH+LSIqJFzPSJc9ZS1ISICQ0rT2MuyD4R4a1sPdxZt4rOWnSVv1/q8hwxy6CNb0Vhu31borD93ovCqbkHX8bdZnYt0Ga2LsdbAspWYqHLx2/0YWVxHX1I/YxdWFkI+7JCfWHyOnx5ou2bzOTWOthbAu176Srj6nzn7lieKyYKaxmaWkYGaMKl8m+8FoVVvE2NmN4ShUAT89eIwlbI9cWotzx0wPWi0CyIz2pMcR+gxfA1YJbD3/IX/9VyIk7Bl6QLl5M9Sm1DpAutjZ1oTdtOtNQljCrY5bf6uXF1tDeAqOuLIKqkt29Xk1lCdLlyfV5Gc0Ox91gAJftTVRHD8sIQl+sQzq6dcfm/xRbaNLbcnLSyCD8TLZY2DHkZs0BHIELhXIoJEkmxgioovYktxtCscaI1oaNY9jRf9CFf9PbALwcvT2+HVuunRK3irt3NfABFa13Ctj6hmPVQGktnrHUSxermmREpaG2/5xtr65zEoTuGJBCiQHIupkwzgLBo7Wb7pNOcxOEt6J0s3qdkmFx3sfRO7oVS37KupS9g3dY1HHxNMCnTN8tq6wmu2tlJTfaIWIpD1TQ8JLewhJDEZk4yEZE0PMYq+iyj2cSnd+lctQLfUavr3EkVhVLb4y0shfaCMDrEMYU3xAGloPV8XoWgepQeyWYxTMc/wo0Rzil2jdUwqGBEEgwitfZqa2WyewaA3E5zWxRuPXiLYMhi8PFwRFDBW8NTnPupCAkThale66EkmrSi0F6Gr3TEqGYV9MliJ6mGoYfiMg7YSx0fIBiRso4tc9lBcWhES1TBI5lwpUMRkm0W8qTDQjgGuJw4kdy0W717R4kL4WPHwmLUemOB5thBe15Ya9HBRXiJxRpbthmppeCsLtXDzTVv7Tg6iYvf0aax/60VXfubjbm2VyqPliyqezcjIiWcmCD0iLjIonJSX4RjW2GjPV6tILREkKlz67b7MUVXRGjMAeZVwNdzz+oM9rV/l+LQNYmc1Q1syVS2XgALCyKARUWRCMH33lNLoUGL4WvkL/zmP5xdtSawqrizP5fjSfp4l6iyKP1QXA6yLJhttAkpS5FYXT9pXfUmA+R4lo3YFQFKeYRF7IrU2lptzcNSkxBYx9/Zdproyu31FtbDlmJpyzdK644ypphDu0/bPbi4gXsx1+5P9c4XXFDInP+66PRaEzGLrJJZ2hyg9hPr7yYIF7tm22oC0ra567bSGx/UZ0HTWYxKLGJbegUoSSkyDKUbyo1jwjS/CRRReCpBxk6uGGs5mnjCpCZSStMsRKF1OrmmQOfCCme/y+A2p6/t6fK6TPx1n+aW3ooFFdW15RFYlJdpYw7thUSCFhHYTm/laVbLm4E5i8J4rpBsMYxT3aY2ecEeiiURpes8sXhxvMZq0otCs0CZheyxfTZZx1dxrP2Qmwv4SsfiVnaoMXJt/N+VDsWC5LuL/1IHXOqwSi6xfTp0lizjCI9RQq1P2Fi+2tqJp6yLxlaMm8Wo9VjL095ya9ZB31n5qtheuuWvE4XttrQiqCW14qvHuj57Yq1NiLjoRmNZxf3zrfQy1lqOprUEl4YMxRNmiSSubMsgVcRticJ2eCsKW1rrdlTBMXhchXHVwSSFTPhNUWixjn/hN/9h3HcoDF8zveDqaxmaIMxGrk0Xc9s+r49zaOldFX33lOI+7mJO2oKpi+F5OWMe7n2yNDjredytX4Zkkdw8qzRZuMx6uDwIeVy7UstQ7ixwJgrt/ltivE1QFbG4/eCSmGP8tsRjULhZ0/jFPCYObV+W6zzlIn4WpYtFX9/Q1eXIrOvi2Vm06G5b8Mk4rrNlb4G1sGAWsL6MjgnGzuJW6NrPFUH5LOF7SkyiEYX9+RLT71OE3splbNvYLfDUHbWdP+paDMf8O5/I0pawPg8sI9kykW1f5rM8fm8vVum7ZkE47JbnQl970K7rY2/Juab1nd1rereniRS7p5j72Iabhcnq/pmwu+j6IvdZwVFTm7pfjmtFHHLcYF+celKHSx2KqDiViXzVxAy221SXn92OeTob31sLiyi65oZwKkvbLHi9y7d0NPHL42ri+1Sx6nVx67TcfeMKdtl6+LzbaUknARsVJhpqFnK6WR/jUDKEl+OXxab7+Sd1eBrGlSi8CmOtMbhKVJHFeWsu5/IZ/MIIUz7j8vNUi9v7yi0xJdxf/uI/9IfK/5/zvf/GIsEkqkAATPm7oLqOgSQqF+ZxQdNFZF5c7DbeFJS1OgKwmXHmJZYMYBd8bl+VahuGKCUA2tY/Nl1SUrHVgBBSlrEbtMYMDjG5QOE2LIf5c8ou3CBNjCKAubnoTGyNMbmSQ3YXR0EcFBgAf5WK/MYxxw/aaswFHSRZONvtUJSWeuqzO9p31sQI+JAzoU03KJaxg6ZRirCwUbkOXT4cotXdXdR/YxW0ZJiIZbhmSmJJx0myaIjewU1WuFkhc0ziMAISAnQ3pvZ4IZZYPD1OORHldriSNYRUo9IsbvOcOpqooPRDXFjwLMg01u/Whi5qXc44pGlOZv9m1725hUXWLw993cEeq3tYrLS+fg9aRWHJKD8hOo2oJYxikXVs7uNmV0oLPIfsQraXiryuspx+m/Nwr3BNjKG5kVUFYxYTx+hXGcK9OEy7kxLrVKXEYqGJbSuFsBsvB5Dck222bSsKW1eyWaLOZCrWQuMy7nEm00LUJSsUsENYdB4B0nWV6uRVMbhDXMStGe00Jjw9NJfISS3vWlFo9RMn9UUUpjZ6frXM/liasGvr5x1jynQ9zEMyKMQa92nHdJAI52sZsp2bG1Fe/wdQMqDLOl11w48SV5nGYzEY1N/R2rs5ST2OAxxGCSkhB8tnhP1uAa6IPI9YBKGTiGMcUmWOfIwP0RfLJpCEWHLJe4wu4BCGUn8QyEaSKCUr+BgG7HxqI+h8jefvDShtPUJrsRdyDkArCiMEc0g1FC+nXTnX/9Jv/Y/xkKDF8A1z6o3jujeRVtTZhW71nNrh/TLMHdDfoA2zHJbEkyJIl62Y2kKpzmlNSAGq9dApxDdlbkrSSrcztgm2iC4reTGdrOdXp7WUjL3WFItddS8D1aq4OrQmwELzaeVrGiQn/J14iU7T9JafU67jxrpYStW003lsl8VBGi5B1/MMDghaExsswaIkXOSDcSqh4wYQ278QUm9kv74FqfUl7i199r1vT2ecqB1YlrtlHSzp9X2cQb/hJ6yYUgXbYl0nRGFqebdxrQtWGe/9+Rj3TWxi81kshDm+0Dy00Tz1Z53VaQyYZ7/ZZeMQhmvdxKtkh6ZM1mId+d6yy67qJ8NhMX7v5tQhxFdRaJy5FHPYi0IrEXPK0pcshcOme3fKMYY9o8TUPq+7uGxaq5u46nZyIlHlFKG5h/ZWVbtPl/CeLiGwt/LZMU+lZNau4lK/cOP3tVI4zyMK2+8m+uw3aUXhKGEhCg1zJdv+1VqVrojCdh11vrQMK/zdisI5uuKGP4ZUysbEoFn/2udg6ne8EeN4IlGoryn8EBJNtqAwfIP80Bd+bTnx7EYBVEFmhqgyzTNa6PVxFUDbc9nVCxGycknYCd/GHA6uxrdYncNTdbNMIJaEmuaBUVro9cKuWNxkqdaKOOw+GzEVd7qwtJVZh/Snz7J9d27fk5z2OtZey6eWYVbEXgzY9vbzyVoMiibrUGshSp/ZMlS6XuQRrfVqDquCzOJdKg1zS+IMe/dx6ZtcYghDni4ksbaVxfyMmn8rV/MzBONqvgZpy+i00/RWwFMuZxOz5XO9Dm2TYxbLtI2wZeWPnIEcRkHYY1WwPZwB84Vio+Yy/BAwz27xcgc0Lc26osd2Xff3orS52718h5wgYm7NUprG3J35e5/x6iTFrfU1AB+5A7xEBMhmgestcWbCrq9FaG7fAIuBq/uVEk7qdWJlaC41KWwTLH2GLbB2Ydp4E5491uZ0LaibMjTNOto4P6C6jMv6tWbN1mVly6ufsHOhuInbQtXL7V2KP2PswgGAWozaMPFtBaq3ik8DWBWXvo5ST1CXNX3nnGhkSTat4GzXVTKXISWuM5V/S8Mt6zhCcIy+bFMbsjVH9+CshQBdyW+cH/mirwEAfPb3fNUi4aRmDm88OHIWWi8Mk9l9nVBi2VUp06qKw7aUQZ92b/XLAJRkmLZVHpAeBHbCzNE1z9uIGM3l7cuz040ht9Vz6aHWC8I2Vq8VhzENFJVU17Bxv5rYVAfEMbdJizlOzyfrnt1q3VStRDI37t123bbYxtW8EHImBNttaP/v96MRhy67gM0NLSFv+9YwJynZpsmKjkiucjdpdtmhJM1ahrI4SVnK8KVdnO7GlJk7TcB+DxcV8XCA7PfQw9JycxNsteiT3ZgSZsYBmI7AuINOE2Qck7u5dSlHLMrPaG5BlxbUC7Pk4hVVKKSM114s2oE91SIDWfA3LeuKe9q57Kb2kBjLOiTG9GOHcLLMTF04tl+ibFjXR7stTxPObFy1kqfhTWzCoOXcHIaaZCLNfcfquZVNaERhhCBGX8b3YtAS4hbbjCQCW2FkiQtmaQIsm7VaCFsLovVDBpKo6xMwehET8lVi7uVqyRKM2e2L3ILNXAZBfRGEU5OoclSPmBvpJddyykAO2Q1tLlIb1h+3/n5t8ZWHMODcT3gv7GosoYtASPNakk6Mktqd5nZuVZjV+/3CmyNJbO79XFyxLrvvvYRVhxLDBLq5hwNqLcZk4K7PHhPF9hsCwCEfu1Raxqx2btHxpSQNNcNaS625k2NO/rDj1ya0AFg8z9rllQod2T18zC31jp3lr81UtvhBcxdfhbpeL4of+sKvxUOFFsMbwgRiixOtntbmptLehLeshnZTui44vF2H3Yh6a6CTWuew9Bv1Ad4pfA6As3Z5FgRtCSnj2BRgteevGUJO1EJcCK3sypUgkCZLWYKkWMj+eW06Mmcsl3ptbZ3D3mKZ53uWd6A8e9q4xN5w1cdP9vtkX3X5uUle1pZ7GahiwOLJYucSXWW+Oreq2+dyBnDJBL5FxKtDKqkD1BZ/U3o46pQeptq7lecAxFhKvxS24gPnXNex1IvsxWP3abi19benlKLpBKd029tup8zLcatwhA4rch1zm7vpPFs29rI6ZyQAYb8sUdNiotBeEPdjFenr1mDJAtVbwrZq521x5pcu4n0WhYt5YfGHS1HYZv4C6wSQ1Bu5S6LJAsXcma0oTOO7QsuZAMllaE7bSWyZdq9tReFiumzZau/TgwuL4tH2f++G749jb2ltacVgO6/FiJr4a9vWjY3V0Jbfivb+uBlborDHRPOcYwb7NoB9t5E+E7oKxyrO7LMXhf12tVZVsy4CKBZCi4VtRWFrAVdNZXHaZW49nx8SFIY3yI/mk694ivJJ2YrDttwNuula2gLYQL3AVhdgE29oQtDcGmUaixPMf0UsduKw1NJy6Y1rHAP8kP4MydnFteWdmdTy+JiEoFi2crbQSZRF3J9ESUWuN55BcddYKqSJLbRYPknu2dYSqB7ZktSJthPCDkBxzbXzaGNBMoFn5UjiNVai4jUv5WqW39t5VQThzG26sFvBoee75bbv8vdWDEaFjMvpboq1W/kI1Y22edmdvBCHZpWbrVBld2KopljLXjjafE3JmPZzq5TMKjTAhHjnzu+nLyVszGrdCUBJcSPVvXyNQKxdTwTDldZQg4x5X+0FyR8txCBbCwE4v1y+cxHHeTu5pMWsdFFlNW6xDboOWfGdYDFaobR3U9dpZE4WPpkX2chnpdvJsmD0Mtlj3V6tiLpFwspS7BRhkr+/F/do2961Yugy7tGWc7naEIlOFE/DrnTkMAG269znlrxjbuU2ZnDnqlu5FY3DwgIbynTtMbVlOqTM5pJ9nIeZJbKN7+sFre1HW2KmPQ6HOOIQx0W8oe3TYh+fI05v0b7umun7YuSb4VQnClGX+MK+93E3/C9/8X/0zO2974ieysC7/dzZDd/is7/nq0p8oW8sh4Zv3gytQCeARSeTVtz5xrKX3BXVlD9ILMWv25tBW69paszq9jZldRhDtujFEsgriNHV0C4bHqS0v4tzql1YytVo/t/EVBGG+dNiFPNwCSivMTIDbk6txOw+6yYplj2J6dMf88HQJj4wAm6uQtLui6U4Nho3b3ZRt3GFKsiZoSjuP82uYMnFjlus9I3do8oLe2z+b4Smy+VpbL9KrboZ8IdYBIKb8u97GeBy+zZ3NReBJFcT5PIKcAK9vIJ4h/jOu4Bz0KsDoHHTpXtj5O4sMg6Q3S6VrGljCYehJNHIMFSr3+Cr5dQ71KzlThSKFEFXrHpbYjJPm5bXvKzY/0Gho8u1BaXOE7ROY113yotIM7xNSmm+l5hRl4ab4I+5JE0cHeZzV6zHx0eCOOZ41H2Ksw1n6cUinKVEpXAek0V9FyFeMZxlYTUGeB+LN8CSIUYfsPdzcTHbg1VEi2ty8ZOtLFzW9i6kNntNhxO7Nz3KSSijBFy4I0Y3Y5RQeiI/8aklnolALzFbB+diNTSXY59M4qFZUC47o7TJEe10AQ47pFqFxY2pqRfyUT0mHeAl4qge74RzjPnN8DJn9XwiXCxE4Zmb8E44wygB74Y9Rgl4GnbYuwnvhj32bsZ78x57P+Pt6Swn4GTrVhOXNzS1B9M4X+7v/TEfXUhhP1a0unnxjyrlN4iQhSjdEl9RXRHpezcXq18bq9fG9KXffM7PlFQFo+9gch1TTkJp6xdubdcieaSZrhyvJjwKWLqNAZTmEjbcimYHTVVBzOX8Y7/tP3zubb8jPFuRn4AWw1uCma77GoJG6N6Mtt5+ZnWbsT8puLaa9vuLsC2PA6S3sjaofOwshGY5bNflXCwB7a5xUTvLrBtS1nKxGHbxSOp1MRhRILMsXbGNBTEOWkrBFP7/9t425rauOwu6xpxr7fs+5zzvWyyN2ii2TWgkrSatNAUSY0wULfhRo2iqMTRCYjRE4w+gfVu0EflIoYmGWJsQbfgwAn7FvgFbQPxBjFV4iRU0sVKhlUZCf2DwfZ/nue+91pzDH3OOMceca659n+fzPOeccSXnrK+5Pva+917r2tcY4xqE1kaPWpXmWLSiaot5qwq5g64XwhgS1O5Gu5jAzFdSdzDClpchapGQ5lGgUeVxuLZxWOz74gq2dxbku0q+7/vwDT8v5ar0onq9CSG6r2Hls/Zyrxh8vRbSOipoohwKoR3sZUiUxkmxCRPdzB98aag9DE3XH6qUzyxrZj/Iab5e1O48tL0j/S7UQw4G75JmQYNSaEnhCMm7GtWWM1I4qlsS9rSVszOslLBKK7xKCgXS8k6USiGAMh39Wm1ByRnGMLHsI9XGckw5t6iFVyWXR1JoISHwlRIe86Lk9C5sWEJR7HYudkCPacFd2LtWeQtl3McN93E7kL/xvS7H3bVjiVYUo4WX11BIvoRQ5e8xhnTHQg0h05YUbtzSCSwpTFyeKwHc/U3U829SkGPXt2KRc1Jo95ECFDuOh6Ib+dzKMcfq4s2IGfY9fQNJ4UeCK4afUfxD/91v0Q+xQL4Q61A5eLA5uHETz0y4hNR9YSVxvHgZDmGBQU2015RywJ56gtq6vNRrzgEpEVi2S8/kayjqoPwVVRmkQgqFCAoIWqjRhZ1RyFunHtYwmqpuqiqWabzWB6ZRC6W9HmDIqTxgwzFHMNdQtEWnItZtoiCGnVWFtK9Nrt0SznJNhfhKGlVI3OU5xmvWYyzv1wfpQw0XvV9eKF1r95P3i0rD7xY1hh8eQETI771Xlj9DyqGQVbpcgBib2bVUVYcIWmIhf5YALrGobERtvSiGo2fhGOKfWcqMECXSjNe/m1H+mCpBHdvtydgxVzI0RbC7hrq830fdni+EXfIMX7Rp+dEAbJ9vH9L9RQZfGFgy4vNaBXzXQo33lw1rlArXphgC7X6xKnHpH6LjvHbcqEqhkBBRq+7iXvPZEp6rjA+8Ex9wTztW2vHVy1fKtZgP+Uo7VvQk8WBhwz2Bs8Up0URKZFwJU4eOyFw5IlLGxgsCslYhyzFFFRRS+OVcfnBtecF92PBevvTh7Eq8rJ8fAA3d2nWtJ3DWFnGAFQLa9nL9x7zwfrnPEZfKW9v3116nrQy3hSCiRFoIoZzlPAJ9YUl3TRL+BakFjSyX89E0P94qlaMgcmYjMxLCMT3L4iff/DxCVwzfRKiv4PDln7XPs9Mzjya5iVxz7/fUDFFbvqG0Z9JrqcUq4oMlaCpi7lrrBZOPGCODZL3NOZRwcXeRAMKxy4MYQWvRh3lLshhU16pn9Tk0OYblIMDMIgYoqovmDg53kVlzgiDehzDTSgpt6Joyg3LJCSum1ea4g2LZFo7nA3AsTqnq4v6sPoiqYpieX0o4Ugoj7lpOIcWA8KKoHeF5r3p8FsD7DlDQYhT94ZqShpsBFILFPM8tnK0fx9yC/D1HtW/s5yyK5tjOcO1l4VMhZAxZy3jJJ6yf0+NnlZpN06BA55WRVy6kUNZd47RAZEvh4JknbgYjKRznBaWQLR3IguQUCinUc1ZycUc77usvts/Fh9PCj83YzYjKJz54QsaEwKnP4dgZI/e9jqfnmayP4I4UqhpHqSOFrWUfd6RQK3tRbHiEFK6G5KnHX608BlpRyDoUjMwcKywkHJ65NCtI3NTC0adRSKG2BjQK4RqSvo59mB7zB0OnQo7PJjn24xCxsupfM6UWBTEeCk3KNcxDxXqetGBL0XQvCd1+0lDiQ7OltwhODD+j+LP/6O/F//Crf8+UHBI+PDkcfcpODbfr12f00xp7LK/Vy9CGpWJNoA6Bq4BTQ9GGHJItRpmRQ+LDQ6/kD9Ix7FyhPZUHY+u8AqCWoD+GlA9G2JXc9Qc/rpOcxXIwM85MiRtB1H3kNWJC9k7uWtrhQgtr6rKQw+f1xvxsCCfXtnn87K6FlA3COy8+M8UoAt77ilQlh+JxeFDdCGA+r+6d9Sb+oOQQ82Ulgyd/N10/hp4HUijXbklhdzkX6v0tu23l3+EzGxi0cN/H3GAdvErHwpFZIcKtZa10NqSwO5+EV2cGixVn9y1tOzfcFN7l8tkdSaHtbTzDWGwiGHsjt/GMC+1KBG3oGCiksFcrGykUG5fSMaQnhQC6sLIUhcg59XhDDvm4HTiqi/LMmJFCAFNSuHPAVoWDGSkMAznXY1YyOob5rW/hY1paT2I1pe7zGIUQWpXRjr9FCvvXWIo2Z51Q/sd/7AffBrXwI8FDya8J/sE//Vvn1ciwHUvajWn0JAOKureErMUrYygaQOf7dNYfcgwr2zxHNsUp+qXUYpSTsHIiYDOhZSFVuzw0qTPrJdMtpSseoX6/sFelsYaZtcBjA6CKXn0ocysGsdXHM79DCQuKUqnkoIacJaTMYp1Tj6nWaeZnqxTLKOQa9578Sqhbilk0xFxD1GErr2F5kEKUBBAQ370CRKDHDfRQw3hEwPsPpRDl3XcBCi2s/BnwOdTK4VpVTUTAupbpZe2LUIAaVqbWHk+seoTExdD7Fj4VStbClmDWTX6IzMigtOqbmXCbjiVK/GSai8IrRD/VUHKqanBeCHklXN+RHwI1lPxO2f36+fJDKtf+yPl5/YCsGctddQuIWY3rn11qYcjactpsARuZdTMEan18i5l1I0hLSFoM8bml5A8Kgfr88oDn8VHVwufhsfMtbJXHexdWPrTB46Bq4mhlIybYmqPIQUPFz40voqiBYz6hqIfv5gsulJQsvltDyWJVk5RYrdpmzRZtaFs++bFe28lJDuJWjZqBEp6VUKwU+tiewVZV3HM0f7MWKrfzojIqMRvIpq0gtoqiEFct3BA/xIGoWWFidj0CKarRaJWGeG9rUhmE9/e1/N1y1M+aLYy03pm2cQQAbTsriCG/jWTQQ8lvOs7MNkU9TEydjH6GWS6HTd4d2zUB6PyyptdAffeDWTcEgYSVy7/aIQXoQ8vyL6Apa/YHYSVsSsIwhJllGEH7IOdFI82aH9jCzfUmGcu/TkE8hKEbGbT5gaIMMkGrmCm1f7LcXedMyDEExKpGEonJkfqmMdIJpRYm7Pexy6dLL8oDi+9W8P1RGaQXLwDUsPITfnqfGmooWL0NR2iXFPMGSg7gWXh5tqyGmwNObIb0B4BsH8mfqH2TNn9TGFKox0+MLBXJa+hUwrHIJK/t3/J+uxZe6+tb50rhZUlqjTK7X1g15iyiMKqLQCOFwFGlAgppeh6PPzyeais3kkKLM1JYjrsY5dB2NGnnG0khUAjlSAof8qohYaAnhSsl3IV9SgqBdv8cSeEsDG9JobY0nZBCWS/K3UgKH9NSFDoT5t3YdBFBH2a26uJWVUO5FiFbdv2oGLawcOz+CeQZcxYm1nHGh1DeD7t+/LzKs2/8Da9FnPVZ9BaSwo8EJ4avEaRy+ezxLetFQr8OIQPufuW1tj9jYm+oX6anikp3ngAAIABJREFUQhdyNDnuaHOhCdCad9jnHwopDGsupND+MyeQEHFH0mTW3hFGIcioecRme6gkcC1Er+Rqkao14nPI0ZBDo0529jJG8aN8/DfesSgN13z2x3wqRHn7WYr0fJ5LlWs4mWdh5RfPmgH2E8bOnxZ420soWUigFMqM5PCM/J21zzPdS+xUoaH7nuydksOB5B1OZ47PgY7jZVlMzIdiFCGFIqINXeOQak7hSArlO5ZTaYN3Wfoffodw5PAlGlNPzjDeKwRClMb2d0Kc0hOPoFkLvHKdVVE9IYz55Lg2Z3A2RoigJYVyDYEY7+U7VT9Hk+1IbRzQXrMW81DS7iJCFMtxslq+CBGTjiXluI1wCym0sKFjqzSq6lfJnhxPiV6O2DjgmpdOAbxqaJk0PLwZIribqS2YGadjLrvmtA/qYztvrYhOC97fj7Y3Vil8mWYOEkFzX8IPDg8lv6b4VX/qe3TeSvUCm5MYJw3Vx9ZSI8YeqeOXWAxE5eYxViVnhrrO25tOSi1XJGdCToWgcg0j50fDeFKR+EqlMoBMSqxoJ4w+hFatC1vvc6j3b6PuyXpL4KRqOGzcFMBkwsyTsO/4bOqURqALNYdUSGnYjXpYjxtSCzFqvpw5h5JNtFAzBzLeh5VoX7lULVfEhwSmVrVMDzvooSk39P4jQAT+yleAGJG/8i74WsLOFOMrr1qmuzsgc2mbFwIoBoACEEMJJ8dapSzkTaqWBbVaWcPJAUp8dV2kZmcTW8EHx9ruLlB/x5HfEHtu4WOBDQ/XayJj2q2KsDUxjwG8EHKdcijVxtvz+n1aCOlSi04uKPN3wPVzQggLOeSVq6l1Bq0ZYWGEkLGsCcuSTCejjLUWi0lYWL7jsfroiVMBUMKxAazjxd2gvNyybrSseRY3rJTVsuVz8QEBrKHjd2IJM19oxz1tJSRNSYmY7X5yoaQdSgQJ1FnVRDCuHLsQNNCqjqXCeSSE4llo7XACtZ7ND3zR8PFKuxpcP+ZVpw+6PWk+oeQjZiasIeGhqouiGj7mBUvIeD+t9b2OGlK2VcSJSX33BFYljMRd2FnyBAPlzre2mVpzp8rJ9nJc0s+CVYVb95HmfzuLItmcvjHMu9d2eYdwtpnOUpiEENrPI9DffsdbMeOtqDp+Ch5KftswfuhnrfJmbvCyPnGYJnrLA2HcX6qV7a81Mr9obVWy/FKz2wU2eqeVy7KdAFrLA63sWB96Uskc6oPPfvsNVBVklD7Ks58ONpwshwlmW1UJ04X6MPU4to7X67CqoFUXZZkbuQwmtGzJHlAIISXu73SWZJrCmdYN4/gy06XPj7Pg+wV8b7qhSE7fOyVhLbzzAuHZM1B8Qpb8lCB5j3yt7fGEZOk0tWWijoTZcWObugOsQbU9v1Yl4/D3AHAegg/lWg7XMxxXzayF2EpqwzJO68u5VMJwAeIjgaWoStTCWmwSqtK+rL0dTetkNC86ke+jJYUADg/scfw123Bs89eTH6VfTsXm5WGQO4UANnuScCB3QB82HotGgKYeXjl2/wQbIh54Pd1uDbQtKQSqdQ7taoZdQsEZWyVxogYKKZRprJXF5RhtmlDyvB9rT9/yPrU8TftDfKxGHkPHlhQCLfyqfYS14CN2y4LRR9D+nfuWdMdCju44h44iNaRcw8e2a8ns3whJjTpLZZp9FQEnhR8HnBi+xvjJWl211l//sy+KqHlnbfXGm8qTIaMnKhMFUvF4lmuopDBmhNjyDSX0peRwyTX/EO3TGgrxA2reIIbQaiVqHFplc6fOmUsSImgLRXSbeUPTpZGy7i268XYdQswJfbHJrX2FXI7JMzDXocR1IBDV604MscUAOw0G2Kg2NvysksTBpy+8U/MP7+/PL/TTAtc3TUPHZaok8SnSp6HlujwaUQvO2hi+TGi9hoh59E0cL4VsSFpUTlEw0U2Dho4lxFyXpTe4LEtF/tp/f6Q95bIIKaxhzTB/vwIYyaSZzLbJveQ63DsaoSzr3019TqsohRKteK++CCnqaMSpV6hl/FXDz/21SchX1MSzMPLGy+k2e1xLCi+040LzMLiFEENbsWyLakZ/WDtvw9Fjm7qyb4u62Ly+sYXc+2nVziNiByPbpQHCIXybov6zP/rVENvkJVqMPx6uRtUr+/Xh7JJfeExfsrBqoU1rAkqh456PCQvj7dErjj8eeCj5DcKtymVZuxhFDziGoa25rZqD2pvZ8GsSaBVgY/s8GSeVyhJWZv3XQsqyT0qhhJWBUq0s17CFZoK9DdNEGiouYeY6L9XLxK2LilHoVKg0ljMHA+qtqniTPEJ9c+VYk9iG2pBEatWodptEHSO1Vn3m2grBRdemT9bLOQuJ5MO8VCpLWDk8mnwoa4Z93ZpqKGbYNZSMx0cgRvC7pWo5S9XyK7hv0LKA7u6akjkzvZYpGf+/uoxA4BhbKJmoyyHUNnc5N8KmJtmYk0PpupJbO8TWgcWEmMd5NGWQ1wCmUmTCsfzLa5mmSwshpwuQKunfXxRT63zHyAsjX1gNrQGoqXWMGXd3JZS7xPbdvzOm1rMex2MBioQbZZxtublQUpcD29dXcuy++vIunsdrF0YGgK+K5TP1uRpSvqcN9+GqIV9RDbXzialE3ngpCl6+aNu6CyWtQE4cVOG7UOpCzAEZGQEB+UAqy3xs5+OIhICxHzNQDK4FNjVHSOAjL7W4I2p1rVjXyPsaiPGYF0RkXb9znPoBZj4Sdpv7lyfEStC12EuLjpMWdwuVEPZ93MrfO2Rc04JL3PU67bSco70nl7jrccdx7fr7z9TYpcSGinXbEDq2YWMLJ4NTeCjZMa9ctqRQYNXDMWw8kj0pRFH7Cg23jOGN+ut3CCO3dnstrGyPR8TaSk/myRafaByWNaFeVBGdGkWRCVrBzHIcIZei0gzKn/4YDmVfVRBlXaBOyZGCFb1GOZb8k+e/zRlMjUC0f2Ws5BcKLAktx+mXBS2kbt4vmFDkKiHH+uC4TCpQn1/Az859DOlZKVKhF58RM+yUwGkIM+ZJlXK3XVTF+j4puR5CxqFWfE9a6ul4LVpp+5Tl/lumhNIWmgznknEtpDz+3YYQ8l2vGurLEzPrwa8wxoxlSa0YbKISaigPJ+3Lhof3iDF8eBZBELI0hpG76z3Je1byZkhhOVb53I4dSvRcdb8zr8Iv574Ay5KzxAEPvB6KYzaOiCgt4GbXa/0TJSQcSfrVH6Mtew5dH+dyjPnzfEYKr3npKo2BIyncc8Q1LfrPhoQl7++aF1zCruewJG/2fNgHtbj5Ec4VwVG9lPdJVGcpmBTRYPY6nBR+enDF8A2GLVABjl+smf8hAP31L+qAkLix6hhoIYvxy8xAZy2g21P7stsb4LbVG4uohZm0KKVTEAGg+h+qgvgYaphWSFjZrIWQXNRDVdg6T8Q6rQSNqe1vj9Ha63HfM9l6Hppzixp4eHbYPwC1c8qyjAk7z1uzsdlPxCwhnYMCGrZKzkUx3Aw53XPh29cEMLf2eQ+1nd5jVQyvGxAC+KGoOlyLVqzf4adWnBIiKJBWTpMohuulFKQsS1H1pNjE2saIx6GqjOg7lISx6jiAUlUOqRE5AIdq9/KZMeR8kDasgmgJIVBCyEIMJfTP1aJmq16F18+V6fZOI4bpHuCFke65VNVfSsoF3Vcl8K78TS6XXQtOiLj5DQ7FaQA6b9OzLidAyUkkk/clOYTPlysuYced5NTV6Vev75bzVIJ0HzaslLQNXvMuTLgP19LKrqp5I/kSUihq3tjyToicFGyU4/aeiFJsAjQjaiGLQCGYRclr+YQ2hPuQV6xh73L9umswOZlynWfvqdjWSAjYdggB0FX9bibMPLafOxK4eV6oeNhaD8BZL2agqImA7adcngnXHLtiRvuZsUqhKI+W4Ml7s9eUhLGQpeW1tmPJ39FenZPBl4Irho4jxi/PqB5K7sosBxGAqghSqDKSQgAdeRzzCVfjiSjb7biucjpmxJgRQu2YEqpaafMOx44pujNXexku6mFV7kTdK/O9qjaCyYRsJzUXkhfGoYb8qFcJG6GoUwkZ23/j+U2YmBid5yGAvgiF+/0AE9K2/Xhz21cLVUQxXCeKUFUR+VJIFT+79ErapTywqOYZ0jOZHu1uPg3Q5dKUu8yqAmquoVz7yQ9eyrnkEY5h4cl4JZbDpvH5rkVPI8E3x+EYDqTQIq/FLomNL2XcjGo45BiKhZPyA+kGVL8v8n1K6XiuM1XPRgvOetGO4wtZSCUUWX8IPho16UX1LRSCZ0PJD3lFQAkDWwuYEaMXoS2as356lhSeHWMkhWMBipBCLS6Z5RPeIIVAf9/s90tdH/uRFO4cEYn1NWWjAgopbNvmBRxS6CHX0BWPTEhhufZ58cc1L2pQPZJCe8yZomi3y/ufOBzO1XLdQ+lWwjafsieFspeTwk8eTgzfcMzI4XgTsG7xQHOUP47rb34Cm1dkW+YBRzXShpBlXIne5Y5Y2orlzggbaNKYVC0vhSjyZSCIQyJ/jqWimddz8qd9l/Vc/VTNhRcCCGqcPSOckhfYVS9jUBGt4kRVzTP5koCxodmPD82xDVwLeQ+hTc1rq8vL8auvvZWfXcB3l67HMmDIYe2xrGHmT6OlXu2TfAgjAy2ULNsk70/IouQkPmU8zZNKdqMMlhnql88whId1aloZ8kLFxHog7EoCxbtQUxNkez9lLdQ6XpSkaQDlu3jWSlOI4FMG+cB5WHn87tsqVAtbnQtgWoHcjKX7ApXWnaT3QbxFClvouQ9jC0nMKKTlykuXAyhTS1bFFmcsEGnnKiHkWy35JEXHkkK1BzOkaOdClNSX0ISArZ2QkLNDOzlDziRMO/P/E2VvRth280y45qiFSaOXYOuGFWuueXvdY+FNZuqECPuagfZZtO+Jk8JPFx5Kfssg4eUzy4kZobN2ATYRPZgHwPiwsG3ytuEGUvwP0akSqVac2RyTbZP2fABLRRoTcipqHGcC5F9dBwB0rcnMWwsva1g5Q4tSNMdvbEsHGOWtvgepX27dTLgpd7mFpDGQOi1OsSFGzSVkXZYCE1USTchYwY0scIAqjyHV6V4KIfTccvydNeTZrouBVHIewzXV928vr6Xaw5DtQFLnuYaa+eFBw8sUA/LDQyFxeULgPgqIAAog0/pOWuLR5QKy9jtV5eR1KeHgtYaYtdgjYiw+EfAaQYmRhTjXPsxWYeRAx77MZsxUbKs+iVyNq4tPIWkPbABIzyq5WYuHYbqQhpF5KcUm6VJyCvMFyCsj32X9gRSe74g1z/DuvvztpOgEAJZY8g2XkLuP4hoythw6BV/uD0IY1piQclAfQ+t1uFDGpfZFfrFclZh8bi3pB89D+azchV2rdgNlfC48IFLWUPJ92HDlqMtK6HjVkK8QQykoEX+/SFx9Aveu8tiaSMuxVkoHkijEsuQQZj33jNiNVdGiHkoxiiWWKyXjR9iTxcykoWOZSujYehcKIZR8vEB82lpuzAsMxHovPnOIAGoL1bpdjm0LeWyrxMStuwqArl3dGFbu3ydS79ty/DkhBMrzR9a7QfWHxocOJTsxfEvxq/7U9yDlcLhZnBHEWx0QRsNSwdhTWcbscuPJ5mE73ixk30oOswmJacXzTuAUyichVYKIRgbpWkmhTDVPsJ1XydneL+tr24ZxQgJZ5o3KOZhp237G3TGmYWHWFnsWSgBP1tttZI5lO7VYYtiulRH2mgawmR8BQg7FDHskhwMxLOuuAAXw++8Xg+z3CyH4WMmhKHWWFMZYlmME3V0aeTPEsLyo+tmRvD7JLTRFJqKg8hJ6M+tReaWn1cLub0VklOteNcxrKCp1DRfv9+169mdlnRDDvTgHYXuH63L9zklP5IURLgkUWxj5bi1/KyGGo3/hiBkxHPPPJD1kVTPrrK3w7uPe9dUdieE7NawsbfF+Ua1KtkbWVomzlcIjMZSq4nHsU8b9AA4hXiGFGtb8AKRQz5/b+R/yqoRqFlIeyaBMx7zBclzxJOzzFMfiDwCH0DHQFxqeWodN7t99uJ468mfHSn/t0eNwlq5gI1J2OkKu0wnhR8aHJoYeSn5LIZL8WS7R+CtOQgczs9Gn2hNF4hZeIsYSUzXC5kNoWc6nVc2ac1j6K5eNwoAAirl8/KXXcuDiXxhYQ8r5UqdLzTtcquR24wmvm2u18zxvsQ/Zaq/lOka8ESXsrOvpeB7iFsam8fKUhB6J60F9NHmJo9I4Cy+zKYhQSCFFAPh+LcpaJVuac3g36btcQ8tB8hAlxBw+BqNsqQaemUkPIeRmaM39+ieQL7G9F5I6MEz1cqj/B7T3E5HO+y0bUgiUkDLlPocwr9Wa5tKWZZudaqvIOqV4fJ1NwWnrJBw4KzTRvOMh18tCiN+iPxqPPwBknSUK7+UL1rAr4fpyKp+XqwkLjz14LSwpBAqB07xD7Xm8HKZXXvBgupMkhO6fXKddHpFAp6QQKITwKVJ4lssnWCh3tkFnpHAGSwrHrlWS122bGoz3fRsKnpFCrXw25zgjhfojv+aui1F1N0beN/NjJZp0IieFrxauGDrw7T/xhdNfk0LoJGlZQkrAXHlgJg0DBGJ1rh/Jo5xtG9TD3eQ4acihqoXaTi/36wGAa6UyixpYl1U9fGzhY0rUDKezGTNRDelAvPrpWAhiW9mFBAQjrAW1rMF8OrOsMT19iblNp+FMc93mWW1D1XotVSVU1VBUzUqqwoMJH6uyuHdTSCeSbTNjd/DDI6gqhxQD8nUrIWDgo1cwS0g5xlKJvPZEleqymncv9fN0aaFDVQcl97Aud+FjGXv2LBZ1cjAFL/v3Q/NgYm1zCsVKSIpMpBJ5fyGKYTuOqITpWZ2+qH+XZ1UZrJXIy1KKuCTsG0JGoP77On7fJWxpc4uFYMjVigfqQlnbbAYwni/tQ34nIeVY1n1+eR/AUS0EmpIoBtKiAArxm/VJtkTQ2ssIhNjZ7iN6Hxn+MJaA2nZ4dnmGCO7GCwnUY01CxoX0FlI1VjkDRTUU26AANipinBR2HEPLwPwH+ugVOCp/dj8h+qO5tP08jPmItmCpvPayPHZrsdsEM6LrhPBjhYeSHR8dv+JPfq/OizUN0L7k1gKi+RNmzT2S/YBaUVa3CcYbgfwBRcGwuYc25xDoSaAliGppY82wB3IodjXYSW1tmqWLTCUMLePL+lNiKM8sS7asPQ6z6etctyuplPy+c3I4EsOXgr1WhhpqWwLZXndT1Shzzb1kIHM5p8k5VGLIDOzV3qZONbQs5HCf5CNWs2ze9xJy3q6FSH3Aew8tCzgl0FLJoJC9ID2TQwkrh1CIIFEjhpJniI9IDOW6I7W//ZCDaPkHG6PtkRhm0xM5RypK4TonhnmtuYULwCsjr83QmtZcQshLwlJzDLXTiSGGmUtHIvlujuq8hfzAA8rTRUhhIC7h4ppbCECJoWy7hF17JN+FEl6+q18CIYbSJSQia/hXjKfb29iIYrGvuSgZlCKRgwegKoBBq3hX2rscQiGNluDJPmO/5BG2NzNQQshxqF62xFBURKsWig2OVB1Lrp7g0bSyk8ISuSbBy5JCm/cnfy/b73jMJR1/FGwp6mdFXCls/h/Q55MzyudlNMQeP2Ej2XRC+InAQ8mOj47xyymFILYCbQwv6S/FieXNuHxmk2G3yTPZ9lsmYg0pA63SUm1s1Iqmxl/j8A9Qu5pmJ9P/u2U30r0nZv9ywSfjpJMFmRCgtuerRMGcR6uFtXq1HeelMBJYVQn797wUqzDyYgijOL0MZs1qiL2ERlDHac3jo3UFETXFbgItFpEQ88u+NnuMWf9mMbCuaidr27xRUeVG4l723DLW7jO20zs51thjOddik1x9CoH2ebJh4tNQ8iJjJIRcvwf1M74sRyITNDe47jtJ9D9c97BNlCQbbgSEYGQ8pNbv91J/HW31A6x5dFwI1GNeDwUhgrFdXUbQddKzePQELF1JSkj5gddOjctK0BZsHA/VxpYUlg4nRwJmMZJC4JjTaD0MI2W8ly6TQr9etdP+9NVixragG0nhGNIdq4P7968//uy4OtZ8LpipC//y8Jmx9/4xLYGG483OM5JPJ4WfPbhi6DjFt//EF25utya3AA5q4s2cGPPrEjha5kgbPSLWQhTtslBDyVnHluKUrr2ehJbFFHszSdhVHQxjccqgHpaT1HWi9o2G1va5MIwNW/ES7FRGky9YxrbQsxpnV7VvrHA+hI4HdOFjbse1UzlG6ApR5uFkDW3vGVTz+GgTFdF4VO6pkCMtUKkhZlEPZUqhGGMTgfcdFCNYVMXdhKNRCKCog7xvugygeBmiJ4mqHqr5dZ8TqYUodl6Uw1qMMlUMl0HlNmRP2uCdqbpCtLXzTM0rlGKTsq0qT896pXB7B3Xajq0h5Of1wyNm1vfVzPquqbVCCKP5LrZ8Lkbm1o1IQ8aTQrRILSQNNNXp2VL+Xpe4q3q4hoSFsrbCe1bDyV9Vw8libp1A+Fx4QELAPW1ICId+xEDLE3zU3MCnw70zH8QR4k8Y0axpbF6dKpjmy21J4Wh2PUKOORaXNKJ6JHJiWyMFKA9763ksOCNbNlx81pHK7vNUXrnFzLZstg+jfNakul3SEs7GOxn8VOChZMcnh2/78e/T+TFkPOYZznJLzjALMXRhiQk5BOYEcSSHQgglH1HDygCgxLDm0lRiKM+m+Ni+T2MYGBjCzXQMLytBo9vkEDCVzMxTcmjPbXMQLdTeZui6cvBNZEsGe0VRuqMcSCHQQsmWHKasx1NiCPTkkKiEmCm0yuVKDO318LaXXMTHxwMRlLAxuL4ISwZFgYx9HqF0RsHSxqonIxGwxELwTCVy1zdZC1BqGDj0OZ5lmxy4Ls4MsmsxSgsZlx7IWpAEYJc2dwuqBQ1VS5qyTjqbAC2nEAD4vv5N7mp18EU6ndTOJjafsE4/CDG07gSaq1iLDdaYOouai+l0IsTwLmwIxHger1ipVSqvoRFAyTFc7S8aNPJgq5I1dGu8BW04V0ynJUxcXsOc2MnyjBSW19wTwxkpPPMp3GqLOTttr6upiTo+96bgQgzf24/FXWP7QeBICs/yA8fxM8wszM5szcYexi97rj/3Hb/79PyOjx1ODB2fDixJHH8RflAbjFkbJOAYsrDrpKXeIZTd5SDW/aUSztjT8LXeqOs6VQ+NShiUONZlIY6TnsVqZzNus4TxpGVdt5/6JYqi1+9r9x+Nrrvnqpxjsl8zfhZSKK/BkIi9L0qRMK2sBzcvRCWB2e4/KIqVTLKuN8qQkD3tXiJjX65ARZXDRRRA8TVsuYd6qrtRQayfjTuxtDGfkdm8hv/7e60tDhr3aZXG9TN6JwpiGyu2NACwvRDVUJZFIazTZ+aDU4lhvJeCk0oQazh5RgwX04mocd9KAIfv7mqW75ZyDu1eUgnhvaiGhuy9U/MNP7cUq5r7sCkplP3vQmuBd4YxnHtG4ixsWFfG3epnfGv/M9jrkPEjSbRt7cYWd/K6RCG06qGQQbWkqcebGUmfhYQtXqZ7zVkY2pJC2/5um5BTe67RXQJwMvgK4cTQ8elCCKIlhza0bBOPuzHmGIzyALL9lUdXfEHKpA+zkRwymmoINJIo3oc5B3BGLVIJxfNQK5HrQW218k5d72NKdV0+VwdV6esUv3b9nbegIY1NcRpCzWyMs4UQDlXLovYJKVHCyJWw2AIaIYUjcZSwtZDFE2KoxxKiOIaUcy5qo+TxDduACTkMVIgjhTINVNXDvsuJhp1TMkriosSPYijHEGI4ViJnBtalWslEIIbW5eVS9rFG1xxDUzNPiOGMEMqHO5tuJ6IUAkdimFcqSmHlDhJGTndVLbxj8FpyCkvRSW6t72oYOay5FJcELm2ioyj6ZZ7r949RPAhFeRdvwxkxjOZ7LNXHZX3WXstCBpdapRwol4KTuCsRlFCyEMDnUYyu58RQjaI1j7CvMAZuE8PuWKBuP0GYPDZGhbCrVq4EUGxoLMaeyE3pPBJEqxbaYpONA/YctZgEaB6F1oJmJITtuP2P5xkxOxs/W7YYu7HY8fZcZ+d2QvjK8aGJoRefOD4UvvRrfheA441p5o9mb2yjZcGW262aUCsbjbLRQtYl/JUZ9UHY36Tl4Qi0h2OQqRSyBEZYU3m4XjKwZPDS2uc1/0Hu/elMwYl6yAWzjaAdSFpBAo5t92Qdt/2tXx5TG8NUChU6P8RIna/e1IPQHLOFO2/8htLXWEOmsYQ6uz/TLSFFSJioc111bl2nptRSzLI00iZqn4xZe4Wl7NeHjCnGcuycm1p4OS96kV8UmqNpwuCn7w0N723i5u/YjeunM7WxtGmsRDEOaqIUlGjFcg0rrwAPRSeY+BSGIbUjpaD1NWx+OAHo0jFeBmcRgMUQOsk5XLWy+PyzJvl4W14QwXjMpVjkMa8m/+5YXDQWjeQbzzvbp9gWnMyuTbwHgSMpFM+/jeOxS4rxA3yKFFpfxvG6M6gzqi7FJFELUHItAnmKFNpuUXZ+HDv6KH4QUmj9C0dSCByVQieFrzdcMXR8bPi2H/++0tv4JM/wzB5jbOg+GrRmc7Mj4q5jClAeeDZfcbS4YSYNLwMmxFy7pfBOfeeUmnuoPodVQdRwr1jb7ENYWYQwU0fRKX5G0QMwKS5Bvz039VCup99vCAmnYk8TtiFkPISy65vaG2izsa1JuVMRrV0NE5VcQ6sgCtESNVF6FJsClaYi7maVSrxF8csJXTu9zSiLmYuiKDmQMYBTLmHjzB0xJCGnS7WoWUQdNF6G61KIcG2PZ3tHs/gOzp6bg3F1jkF9JXWfUFvdhUIOlQRGaEWycI39vm6reYVprfN3XMzSq3k1L1zm1wwElB84AJaLCQ+HVsFvv2NqIo9jPiERa66hEEEbRhbFUFrg3cfWag/1eKIYPqvb7kxoWeZlKnmFoujZtmvW9kW26fs8kKpU8wkfeekInx03HkPWybk59rneAAAgAElEQVRnGI2ptbJ3eNyM1dHv5ZYTWEjUsdBESN9mvAkfjHq4c8ReK31HomVtZ4Cj48NoQTQjgPYVjJEbMvPj/sAxZ3Ekn04EP5NwxdDx6iEq4virVZbPbipnlgZqVVOn4qM1imSqCNZ/cr5mdZMRIyNKT9mqSAYJyy0MrAzpnCI2N3wphSP5UhTE0S6El4kqCKsQoSl+aKoQqJA0pkIWOqsctHlRMItSSO2Yuh/Vf6iqZa9GSSgzr6G0X7MKo3kTaQibKjEShdJ08tAii9DCykr6Qmj7E1VfwUFJXI6KoBaUBDMNsZA9zRes54+hhI6B1n1l/EBYa5nZD19RLidFObo82++smwlRpzDnsYgFUKsiVQZrCFkIIVBIIVBI4eEcS3vfAYAiHzqdNDJxXDdeOd34AWdbWFoCaCFdkMq2Mkby5sa8OqCpc4+8HHoHj9Nx/mzbYy1OEZXQksIAPhDLSLkjhaL82X8bRyWHT5FCAIduKTNSuOeIPUc85gVbDRtb6xigkULgqL5Ze5qZ6jfeV22kRn4w26jM4bNQp2ek0EaCZoqkk8I3D64YOj4RfPtPfOG0ys2qiocWeMSHX8NCCGfE0jq45CExulnfmLAaUxW36s0zxZJ/mEKrYL6284eHerM2Bti2PR0l6gpIVFmU57jJJezHsSp9TGacKnH9tFsn+7IoiuVFipLIBMQrKyHVyuqdNYRNWuQiqiM3pTIVYhT2XJRFa5AtuYXmvkGPe1snqmGu89UsW2FJl+QQynJVBXWajKpTzbHBuUyBRjqZK5mkYm4tCIaghtBsauo0X6yFTd0vFBJenqDUq4YTYqhdTSQMD6jFzalaWEPFeSUNHQPA/k47rhSe5LsaKl7rD5aFS0rEJTX1LzSVsCmB8ha0PMPy8ppSKN+rJWb9TtpORwC6tA4xtpb8wktM2idZQsnS+QRoIWdREu+H3MIA7siXGEILLCGTbVpxXLuDSM/mLUesIWHLEXdh1w4icnybb2jPLxDlcexgIgqj3SaVxNLNRCqQgRo6RmsnuOfYEUQAXT6hVCFbI+vEQc83EsFb3rBTaxozlsyy/SSP6T0W8gN7dl4nhJ95fGjF0Imh4xPHL/9vv1/nZyrFzC5jFhIZMc+l6bclk8OYJxXMQhDTXsONtUgFgFYwa2hZ2uoZMmiiZodqZht+HsfbzinWy1BwKHLpwr5yrD7EbNEVqlgbG6Ba6FB/jBoujtIR5Woe0NsxJEzqfWjCxkrwDDkEWkh4xEAQteDj1j0pzUOAmqM45DWqOsisoWSxrpEClDKuhY9tWBkwyp8xtpaiE1VpL+YzpdXIpEqhFJ5I6Hh7xyiHElZ+0V73XiuRc60+xlLVbAB0qX+rRaapu7zy8k1V8kAK7fdNwshSsTx6kQKtKhmAEjEJKUuF8rO4KXG6hLYO6EPL9ybPQr6ZlrwBTREUUga0/MSOJBmlbjTLttcrRG6d9HTWY71EtbO9hjGUbFVSIYBSWAK0/MBrDR3bELKGrW/c88aOUkCfagP0f7OzPtd23K1766wAEGhRIcdrASeGjtcDQhJFrRgfWmOujGDMMwSOyqE8v8X7cCSHo2oIzCuYAajFja1ipkQlB2+j1l85tW9feByIYa5KILflzuDakDxKPM03tHmHUtHcbU88rYIG+pxDDo0E2jxEoCeH5ZoNqagEMWyptNdLAyFM3FrllTewXm+vNAIoBFEUwTqmHIyKIhhqdbKQOqlmzrkRPiGP+17GpQRaV+12QutajrEsZeyYX8h8sKyxnoZiLdNhUA77fEQhgaEVBi2kyuz+rB53gba9K/6EErY3HU4uxrOwmllzzSkEUAqmANCSS4pEDSULMZS3UvIMBTbVQiBEMJAQxKoaynbT73xUDUslclMsRV18FjckJiWEMgUaORyJ4Vm4VjCSMmmDJ6TwjBBKL2XZ53m8KkE8M16eOSoIsWyt9nqFUCxl5Nhle1MIRfUrimDQNndADdGa3MGzSMls3agkji3tZmFh6z0reYyzQhK5Dnk/XBl8bfGhiaHnGDo+VfyFX/s7AUjeSssH3FLAVslZMmEVoN2ozyrh2rijYgiUGyEBXa7irIK5hOSq/cfCCJLXFblUMEcGr+Vfvs9lPpr8sbtSIJAvLb/wkENILXdQ8gVL/iJ1+Yha9TxWOduEIbQ8wEOLP7RcN9ubV+cvx3uGEB4JjwIAa9eOSqJq/iCHAA43bh+zNnESyg0kjKSFg4fqZLuOJuchO772R9ZuJ3KMWa7h5LqsOfVZd5nuozZ5aR1x1Ip10uXSCq+FjSVf1CqGk9Q8U+VurnF8kIvaPXnIzzBaTAEm7QI9KRxxlnMItBzDPUdE4q7iVrDZKtwbpFAUuIe8dqQQwE1SKDl/e25VySsVcic9lm0F8eiiMKp2lvD1/wJ2jkhMnUn1VquQr7nkUl5TyR3cOapqqDnXEzsai9nfc4ykLLU46BYpHHO8bXcpWbb/bOqBk8K3E64YOl45vvVPfP/hl7GEuF7mQQc0ommX7XQMtxzMtHNTJFs1MwB5EEjnFFvFDIBqqz3ai5oYrsYDMaN5IKam6IWEojyq6teuO+wA2Ob+tW3TfMMEoyL2IekuDC3h58HcGoyDesgExMesx9WWeUY9VIK155LzOGmVJ/mEojJyoJ542dAzUalubrJvUftSLuRRVMxA/TohlnbdaKEjYWRRIJ9dQIm7UHI3bz0Nxy4oaKoiMSNdAigx8p2Eogn50kLIohimFeCleBgK+U93VTG85zpW/ghcClAIYOl/TCjhZBmytgIqopJ3KG9njGV+WbL6GQqi/XGkL61veyenA1qOoeQOMpOGliWHULwNA7hTCSW8fJlUKYvqZ8ldIFZFUJaB8t21ip1sm1Xe2mPa17NS6kLPcg2yz3iccV56MEv+YDbhX6k0buFiCR/HQ+6g7jPJlxaMZF3uVWOjgEDc5QcSjvc1UQatD6x0vWljRElk/eHueCPgoWTH649v+eO/rVseH14jZonXcpMr/OF48x0Jog3rjOFlAMe+zJlaH+ZHIR6kxFAMs6WtnjXFllBztHY2Jv1uzEO0hG0MLU+7p9jvspk97lvJITdiGa6Sr9i2SUcUNbuGIYd2nclBpC0VS5tJV5RmrzOEku18X1JbptWzUEPNFrKthpXtfnxZtRNL1ytZu50YZXTt1a0ZOeQ1ABlId0Hfy7G7CdC6mYhCmyO08rhtM+u0s4n5PN+bULJAiCExguYX2g9CTcsY1HBmwrIkMJMSvWhCzS2sfCSGwZBG8SwUYpiZ8Hy5Ys+hyzUU3GnLvCMx3HNoxtkm5HvopmS2jZjlG9rxcq5Z/2DNZxwIp52355SQsSWFQgSFFErIGAAeqnqaclBldFocMiGGdhwDnc3QrLBE1sw6pJy1V7dpNwDwU//k75gPdLzOcGLoePNgiaKQxFt5NzMn/qYYtnXW6NfefAEcuqgoIeRKEpm6dbybauadgFzzDytHiQ9BD6yET0gjF2LYeiFX8mbVRAZoN11Q5DXnFr7U/fMxFNoVsZgXSdbH0OzbvBorMSEgXrOeK24ZHAjx/b0VsFxr8YMhi6IkdgriaG9jr7P8gRqpq+qiWOdQym3dYE9DKYOrSqjzQhYlJ/GygnLulUEhhgGNhJIQt/o3lnxCZuRLREgZ6U5yB6pSmIEsPY8jFVVQhMyqEKYLQVr/7s/r7ktVEas9jVgjAdBqZASgs1GS11yVwrhk/TyGmNqHQvL/hEAa1VAImXQ/yQysdV77KZvcX2t0vVDJb7Rt0i61T/Hz5aok6IV0OTFVymVcI2VCdiTca9fZPD9ZBnpFUEjdnqNa5shxhBjaSmQ5pl33mJdDh5SdY62SNmS/5gsCfXi5EcJWcZxyb0VzZicjr+9gD4P+fgSzDMyriMfwcDt3n9ZbzgtXBt8OODF0vLkQgmgVRMHocCLjgD5HZzTF1v1POrIAo3IoFc11XCWIQg4BgMXmJhFIzbLbcaXnshSnAEA0fZhHg+zO9qYStUMxy2Ra5g2JmPR47rbLObQIxYwTWxvTT1n7OT8ac+UJOZyGl4UcTtQ/DTkvjdB1hM2qj3IMWde1pqPDOlUMl3A4Ntd1HIO+F0oML4X4AU0VzJeg758ohvu9IRFiWl1b3gHA/vyoFKZ7OReXTieGGPKdea2iChpiKIohtdP23YAGYiiqYXkLhBCSqoVCDGVbYsKdqovtuFLZuwxkUSD9k4UYBuK2j8mZGBU7WfeYF1X6Zopgr+y1fQNlSEs9Oa/gFpG0x2r+iqTjdLt5xkpLOyGFQOtzbEnhGSHsz3tUAC1mhLC/r80JYX+ONv+//BNOCN8iODF0vPn4lj/+20wVc1s/+1UM9Hk0s/xDO//UDTrZvstGPbTbOvWwhpkpN/UwKAlsdjexI451ukPJSBQSafIONTLHZ6RQpuehaFENAfTK4aAYjscJGkq261KbSi/hSgptmFlJ3J4aMRyqjWlPqv51mBlWz8hgzqXIZU+tIllUvwxjSyOqLet2USKJGwlUf8MaQu6Oh2ZDk5fW7m573q5rF/J3MSSx5hLyWkhhOY6MY/2DcCWJWFq+p5hak/RBjm18FzK1xCgKSTREXohjTPodEAXxEpMSQ/leSPhYupQAjSRejA3MO+ujzktY2drEWHI4hniFGMr8PHzc5wkKIZT58bhWnRuLYey+W6fwlfnHSf9i4Gg9Iyqh3WZdFG4RwqceYrOws/UVvGXZZeGE8K2EE0PH24dv/RPff3P7rAWfTcAe82xmuYiMI5EsxHAIMaOtA1NptweAH9oDRQtVJiqiVQwpUylk2Qyhq9ttzqJY1YzFK0wtvFyOJ/l7dp3ZR024+aAeMpkWe5YgCjHM3AirURVpS+XaKjEseYeVLNpcQCkqsSreqAzKeju1kNAxAI6xhIxDaB1d1qJA8hJUAeUlAInBq1X/IogZeWnr0rMFlBn7fTTrgr4nQhb3u3Zd24v6t68h5HI+dN1OyvZiT8MBxcAa0F7dCGjdToDSCo9RuvTUtysshhzKOinasopeVQ1nFjazvMPMwP1aPoyp5gIyWpEJM+EuFhNpqyBewo5LaDY2gfLB+FoMqAF0Cp2GriHh5dDyEU1oV443mk7LcVbK2Djo8baqDuZJaNcezxLDq/Y5bmFeSwyvJpQ8ehCO82O1cToZN7tXjcRyRgLFmmuEE0EHnBg63nbMSOKoKp7Bqo15uJEDx0q/Mq6RwVY3YRREmTddVPBY8+B2Qwwf5/OiJIYmwGiIWQgiYIpGDEEUgmkJYl8NLOvMsY0COFMhm5n2E+qhDTlfJZR8UqgiBNBWI7c/RFcc0g46udfJuJzBsVeFiFv1sSh6lCYqYWZjyWP+zjWfMFkDa0MC0yyU/Bxm/qge2iKTJCRxbV1plCRaYiih46VVH5O2eGzrLPkTcmjXSa5hjM0rVIjcZUn69q+2krkLJddcPuNt+KyGj4UYAsBifploN5XRaBM4hHTHfSyRsuRxDPsKIQTQKYh5Qtw2Pn6ubimCdl0ZK+knR7J5RgpvGU7LOPs+231vRT0sKXQy6BjgxNDhAHrrGxt2Hm+qZ9tk+9Rg22wHGjmUddbyBmgEUdXDREAmYGsFKZRIiSIlqHF2CyubLiuSezgWrAzhX8roOrJYM+wyg0PxClMleBJitQRQKpWtCjlRCcNezK5pEmYGoARRyWEoOYkcwlFFtH+Y0ZpGilTGHs1d6w+qimAl45YgivhoyCCljFyrlZVA7hnpeVknyiDtjP150PdlUwPr3oZGPSTv0LqfaFi5kUBRD7Mhhoist3SWvEJbfRy5/JmC8fYkgAJ3TwKrGo6VynadDb2OOYe5KooaUq7zEloOYNwtuyp9tgBF1tkCFD0PZSVbF61Ujjr2MS1YQlLCdhd3PKYFd3E3RtW9EiiQ49pzzHKIZzmEohLaYpOZv+EtUjgjhDS5BustKB6re4pdXrR1V7CFc631ITshdJzBiaHDMcPLhpsBTEMys/CNtYXoKpwnIem0h8N8NutsX2ayPZpFMdxOFEWrJFYSGa92XX3oV89EoFcAoyGRen5L/IRvTAiiHdt5IM7sbawSKaFkm/toK5mle4k5j5K90K/X7ZP+xZ1qKG+ptaaxVjSVbFq7mU4xrOulUwnQq4ei+O1mu4SSASA9Q5sXxfDO5MCtZl5yCs3lsyk4uVWEAqDlIU5CxmBCqGRv1jrPfva7whIpSjEK4Z2Zl2ITu4+0y7PrrIm1DTXrdZjzz3z7muIWTlXGEbPwMYCDfQzQ7GXs+i3Ni0/2qZp3Ekq+URQCzItHZqktFv/rP/Xv3Tymw1HhxNDheAqiJlobm1lOj61+FpsbgeyrYeMuvCwFKkclEQA4F4WAU9tfVURA7W6ARhJLIUoliUYdjA+GME6UxL7yWcKzMDmEbbsNBTcj7LZqJIdMNISPe3Iob4n2Xa6m2eW6kpKxMNjblHMZ8+yZrU3iRtxOOq8omZy0rgNa+Dhf5qbWQgbDlrE/q/lkplOM5BOGBGziS1hb3QGlIIUHlRAAdmlztzTVUskgWgi5HLyuG5TCcTuZ7fazrGHmYKRhm2srPZKrCXbZvx1eQsVLzCrcXpYjObxfNlXPrDJ4HzfsOeJ+2ZRIXULrfWzVvBlptDgLCQussmhhlUvJC1xC1uuRfWxbOqARQr0+4ikZtNczK157qnAEaKqg7dg0u984GXR8CDgxdDg+DP7+L/47APr8ng8CnqgDs1zEMl/VJ6siGkscNutZC1WMomjzEU0BiyiFfY4izNg63RvJ64mhmRclkXGTJHbqYbd/Po61qmGnOlplMR/XTRRCiy50bDmi8ULUsZYkiu2Maf0ntjPAQBIrIZyphEAjiVZR3F+07VJtXOZrcUcTpzrF0BJDVQrt6+qIIfdT9OqhtbHRsOMktGznZ16HQFMNgfY9sdufrc3UOtLxGEu3rnyYrnnREPKMDNr1ew6nY0Z0ar5VDHEkc3Zs16pvEhK2mCmCZ+pen3N4HDvbz0mg42OCE0OH46NACCIwf1DNlMUzJGtVMSGGdj3nk7GbeSg9mspmqWLeqeUjmk4qcRJutnmINgStSiL6MLIqkJmnCuPMI7HsZ8iGeiK2wpKONHbh5Dp2Szq2zPdWNod5QUbJKURPBm2hCpv8RCGGabXE0BDHxYaNK0m4O64r68vUhpJt8Um2xFA8Ck1P5E4xXCbEkGBk3rZf5zFU14e1Wc8ESyK7qvxK1ie5huO8/bEk81Y1vF/ah8J+P6xyaMmhkDMxugZ68mZDxHPlsFfuZshddfHtXs2zTiExZFUM7RWcqX+3LGnKceff/xmcEDo+ZjgxdDg+LliSKLBdGGLINxUC24/ZhobYhKVtCz5m6ghity2X0LP0ZsYeWs7gboy0t74fc7G3IVXzgt22N6JXSCMbyxr0ljj5JLcwmfDwDr0FhWsrnAhbKaho1je1i4oQxT2r6mdD0GEvY4RIMjUiyYFKdxMpGtlSm99zV0zCa1QFMi/FvFq6loSUsd/LWEbq5oNe8/6sWNxsLyT8zZpDGDZge6cSg9iUwLxCQ8m8ACw2MAuAgFZ9TINKuLDe1bqK5MD1zZONZh1Q1guPrLY2ABAWBucWVgZTZ1MTra/hE+RQpkTcVSzfL7sSrLXLQ9w1xPxs2bR/8PPlqnYwC2Vcc9ScxJ0jLmHX7Zewa+7fuN4qjte8YKGMh1pFHIiVAK4h4f29HGOppG80hwZsUQnpMZIp9JB5WS/vx5YCYuCu3/Ce2n5jlyWLQOxk0PFJwomhw/FJYkYWgac/hFMD2kFFGAmizutBCCw1BMZYm7dGEsEEqoUqfZGKCT/P+jenFgqmZHIUrfpn8xWlA8ukOwowCFlbI3PxMeu5dexENeTQchM5kFri5DVohTOvASRj1qD75ks/LzmMaW1ELxs1MFsLmqoC7ndBifD2zBDEe7tfew3S+7gPGbf5bAtNRKwKjSyCLImcKIZ1zAHUCDgYrY9ypkIOUT4rJCSOzTyaosiYF6BYcijG2MyES1UJ17Ow80RlPAsDS9h5y1H9ELccNfS8cwsh2zzEWVj4cV/03Fu2Poahs77RNnwp6vUlE6o+Cx+fdU56SjG0+Ev/9G9/cozD8THCiaHD8Wnhm3/sB6brNYwXTP/akDWPMARWUmjngXnhSrb5hzZv0fgg8tWEy5IhgQ8mRFcLVfq8w5N8xIlvYhnD0/VdEcpJuHmWwziOaWqmIS82x1BUsGvSUDBtWee7fEI7b26NNs/Q5hTa0HCehJGBvogk2X2fydi2fRY+Luc312KIXzbrbR7hKTk8y7eTgpQ16w8Miq34xBJDsn2XZXebS3gyf5Z7aKuW15N83a5N3Uvk9Nqq4BnOuhWdh3VPipW6YpIyDYRO+bOh5k49fCLH8H//zn/35mtwOD5BODF0OF4FvvnHfqDLOcrqOYYpCRwrovt9LREs0zzp11zWm4eiFK0k6hVE8Uc0Xoi0UUtZswTvSs0027TkC1vjIbRDu4nEh7ZvvLKSsTIvx+nX23CzkKRwbVXGYWNVz8JjViIXH5MWf5wWpEyIYdiy7he2jCxehKmEiOUaRQ2k3BTAsLXWdmEHtnfaeyNjaAf2uj5HEyJeh/BxvTSOhqiGForvOr+ZsHAXTia0O54NI9MJebTEz7TRs+vYfD4FIznUdnlLy19ca8s8ALgzptiXJXXkS8jfXUzY6vr7ZVfCt8aEx32px9l1voSCa3g65Ol8NPmGMWRcd8uuj5gpfjbsG0xuYTSEz77thEYWY8ydp6D9HjsZdHxG4MTQ4fgs4Zf91+Xh0FcgonNikQfyaG8jIT07343Jofvwsxho5xaKRjZ5iZlaDpopWgFDyWNnmG1Ntfe+wEQLWmweoslRnKmJOZKGoDkA8bHOx9ZuT5S6eG3kMD5yW/+YkS+k/owHTG6B6c62riMTFp6QwtSIYE15w/a8vS+dqmjUQe1xbCqNO5XQViB3oeLz6+eJLQ1gws92/BkxHMaS6Y5ifQ7J5idOchEBIC6Sl9h3TtHt9RoetwUv7ov8/LhH3C3Saq+NfdwWrHX947bgbt3xuFVyWNvwPW6LFrlc96hq5Jai5jFuKep599q2DziqiLJkc3/3FPT6Uwrd97C13LZqPun4nIN+J/+Pf3YeOXA4PiNwYuhwfFbx9/5XLbeoL7Cdk8CzsFS37+CrpjmI3Eggp9B/S0ylMxnjbNuiT3ISaYMWttiezWVMm7dEUApZxvVdKNl6ItqGGHLN8bxamro3oM2mC+k+6WJJqFUR23gbLt7vG2HcnjdCnJ6R5ivuz6m71i6HsApVnQXNLK9QrkPIIZlrYigBZEJP6p4ihmzGM/X7zuxt0BNDm3toLW4sMbQqonUI0iKOFHB3Zy1r5rfmMa/vqXDy2Q3ekrmUguZAJkP4gCFEPFEMcw76Osdjdp6QZtef/ufmucYOx2cQTgwdjtcFlijOHli3cIscWkxDzR0ZHPKtbP9mm39oQ8a2kMUSv84ux6zf+6/oSCJn4zoBLM/np5Y1KITSKnv2iW4Jmoar976opMsTvMzXl1CyGVeJ4ZkFzYEY2irkTiU0g2YEECfEUMbbg51Y1ViS2auGx+4pQB9iHrsC0cm2eEIiu31PiOMtG6ixH/nsOOc/qF7ih9bJcf/PX+dE0PHawomhw/G64xv/i2JdwUNoGTh/mI5jRq+0jiBqCNmEmRmt8wqMeqi2Ny3/zvZ0bnY31BO4wTi7s8jZmrK2PHBXACIhXyGiMs7uEx8bEYvXfn+LsHGnHApZDHsjeXa+nM+QSKP+UTaFJ9mEkbkfy7ERNyWJ8hZb1c+Ebu0YOYbOW7I3UxD59rEOyqH8jWcqomyznHpptjcUGRRYWznGJXefK6s0LuuOtNdcvSVpjmyIWX8ExXoNavQ+VEPvsn/Muj+FrD+ExGanO3aaH6vL+TXHAsoPK5n/mX/h34bD8YbBiaHD8abhl/6x36E5UzZyCJQE/bPKZ2DIlzIPcSEBzDDk0JBBHohhbkUsYJtX2MhW8zykTglEroTSFLVY4phXYBEiuADL+0O+4SN3BFGVvmvbXyqqheTFayNrtLewrx0npLYRRlYDa8p9qFiIobyGdOnVy7yW90U7mFCvHgJ9yLjLH5TxWpBi9jkjeXUfJZ3EIKmMDRPSZ/cRWKKoYWh0pJOGHsyiOoc1Ie8B8dK250Sahzjm4wFAiKmRxdg+q43gRSV08uMlxL6C337ugT5nUH8YGfJIIZfjAocUDSLG//Vdt3uoOxxvAJwYOhxvG77hP/tdAIrCqL50Q9oWZzIZ+GQewIPCkqi/jYxkUniPrXxGrxhS6olhM9ceyCaP2+vm2Fvn2GPlpW3Lax+yzktPFoX45cUcW8iiWbb+jFatS/ft2oQkymvMa7v+vDZVr/MgNMStI356rpYTeG5Jc9yvI3FyvRmN0OX+OihRW7bj7LHkeLJNX9xAHGV5lXxE83pTUCLJiTorHGTqjzESVnMNcsyxewtn9L+MhBSLGij7ZdL5v/ovfd/8PA7H2wMnhg7H246v/09/t853Lfy6imbqcr94JISWWGaaz6MvWMEYTs5DCLprm2euhQbvwzrPy5G0yTIvLVdxpiDa/SzyeoMkLv01dvmCC1phjKw3xFB51OiWMqRwdvmEZwTNvCY9p+HvpzmGlnCptNbG8UmeYVe4cnY8oA8/7wSIOpoIsMpiItAyvy1b0qYvaPZacPzsWuJ39qj72X/5C/MNDsfbCyeGDoejx9f/4UYUS1+5OVkcSd/hm3U2dlAni5po8gb3fltfQNKTMUpn+5ldwrDNhnRjv80qjEBPzCxJlLHjeXQ+9vv1YeSW43ewn/mwxHD0YLb72RzAPBtryZw96A1SefbosKRRTn5W1JKHbcPyGSEsnVms8jiQx+56erLqRNDheBJODB0Ox4fD1/2BH2wL47dqXE4nTAU1dGmrVX8mhycAAAkSSURBVO3YdLxH2ernLjcRlT9YPiHbQ7/f7BoteeMFoG1YNsrkqBSO19GRxIGIWcWw20aTt/Gs6rhe/4HkWdwgcnwrl7Aeu513WD+S1zPyCJQfAWdh75EY8jD2KZjz/tx3f+/L7+dwOG7BiaHD4fj48HU/+oPnGxP1pCIdCWKHs/C0ju/37QqwGae2NRx6CxyOfYibIyPY5cDHa7OniqzejTq+LueFe8UTQ/jYEkhgStA6cjiEUm+22rU5i+badHMechXHO+OtMLHNV8QQdhbcIq23iCnj+D6Y4//cb/ieGwd2OBwfEU4MHQ7HJ4+v+49/T5k5SGKTe9AYas6DlMb9fiPxooxjbuMtkmhCyYUk3tgeegJbSGG/P9BIK0d0136WC3iKYF6qVt/0Qw5v4UiynjgHj6Hf7ljDH2xQDy2xpKEohhJNqqkn57LnNyriz/3G33r7wh0OxycBJ4YOh+Ozh6/7/b8XwEDCKqmw6lyrnB4OIBxKxo7EEkfO03z5JmHtsdBmyNs7YCzmkA4zUhQiw26Fe0dCaHFy676pIFrcCgefneOQ/1j/Hjz5G03W/dy/9lte8uIcDscrhBNDh8Px+uHrf+SHjmSQMP92H3IJ57mGB4I3CWmO+ZCybgTllyRpJ2reSxE1nJxjli/4MkRwtm/APMQ8vlcB+Nl//TfPj+lwOF4nODF0OBxvLr7+h3+oXzEljjRRD+fHG8cxnSmG83vrbP+XGgicE8NboWTgdtj6Jc7zs7/JCZ/D8RbBiaHD4XA4HA6HA8BHIIajLevrhA/9oh0Oh8PhcDgcRzxVS+dwOBwOh8PheEvgxNDhcDgcDofDAcCJocPhcDgcDoejwomhw+FwOBwOhwOAE0OHw+FwOBwOR4UTQ4fD4XA4HA4HACeGDofD4XA4HI4KJ4YOh8PhcDgcDgBODB0Oh8PhcDgcFU4MHQ6Hw+FwOBwAnBg6HA6Hw+FwOCqcGDocDofD4XA4ADgxdDgcDofD4XBUODF0OBwOh8PhcABwYuhwOBwOh8PhqHBi6HA4HA6Hw+EA4MTQ4XA4HA6Hw1HhxNDhcDgcDofDAQBYXvUFOByvEl9DX8tXPJYFIpBsIP2vgczydPuwjoaNN5efOm4Bn+335PkAnhx6vi8AojL+1pjDtd0Ydzb2JY5r8YFew7D+pV7PrXN/0PNOt/H5Pk8enw/rp7uQTPrxNLkW+3mnYdu4jz2e/ZjKert/2T5bPx6Xb24/HruNp+74/bh+v+Mxx3F27Pz4w7Zu3l6/GUeAuaMcxh+3lTV/4S8+/klm/g44HK8ITgwdbzWueMSvCL8aFAigoFMEKkQxVFGdCAihTglkxhy2y79uuxk/2c5ERb8321n2NevtOG5PurKs62U8umNzXSfblegQgCD7zbd36wjlXGfb6rScu1/XXYNe+2zcyXa0Y8/2e5ljj9fWzQ/7t/V8vr/djnFfNmOgZGncTnasvkbZ1k+BNp6I60fDbuP6kWjbQ53X7WZdQFuW7bpNt7NuD7K+brPLul3XZ11u+2ZEHVO2R91el4kRINuyHjtSRgCXKWXEegy7fzT7yfi2bzl3GZPrmLpvHR/RzhFRrkP2j/V1Rdj5enwwAtUpgEhABGlYLhIhgOo6qsuhnpsQKSAgIH7tX/4aOByvEB5KdjgcDofD4XAAcGLocDgcDofD4ahwYuhwOBwOh8PhAODE0OFwOBwOh8NR4cTQ4XA4HA6HwwHAiaHD4XA4HA6Ho8KJocPhcDgcDocDgBNDh8PhcDgcDkeFE0OHw+FwOBwOBwAnhg6Hw+FwOByOCmI+dBB1ON4aENH/BuDhVV+Hw+FwVNwz89/3qi/C8fbCeyU73nY8MPO3veqLcDgcDgAgoi+96mtwvN3wULLD4XA4HA6HA4ATQ4fD4XA4HA5HhRNDx9uO3/+qL8DhcDgM/J7keKXw4hOHw+FwOBwOBwBXDB0Oh8PhcDgcFU4MHQ6Hw+FwOBwAnBg63iBQwe8jop8hor9IRP/AybhfTkR/qY77fUREdf1XE9GfJqK/XKd/W13/y4joJ4nokYh+86f5mhwOx+sNIvoOIvrper/53sn2OyL6Y3X7/0xEX2+2faGu/2ki+sfN+h8lol+oPqwOx8cKJ4aONwm/BsA31n//KoAfORn3I3W7jP2Ouv57AfwZZv5GAH+mLgPA3wTwbwL4oU/msh0Ox5sIIooAfhjl3vRNAP5FIvqmYdhvBPD/MvMvBfDvA/jBuu83AfguAN+Mco/6j+rxAOAPoN23HI6PFU4MHW8SvhPAH+KC/wnALyKir7UD6vLnmfknuVRe/SEA/4zZ/w/W+T8o65n5F5j5zwPYPo0X4XA43hh8O4CfYea/wsxXAH8U5T5jYe87/yWAf6RGMb4TwB9l5kdm/qsAfqYeD8z8Z1F+sDocHzucGDreJPxdAP6aWf75um4c8/MnY/4OZv7rAFCnf/sndJ0Oh+PtwMvek/4aADDzDuBvAfjFL7mvw/Gxw4mh400CTdaNfkwvM8bhcDg+DnyUe5LfqxyvBE4MHa81iOg3EdFPEdFPAfh/APwSs/nvrussfr6un435GxJ6rtNf+GSu2uFwvCX4ebzcPemXAAARLQC+CiVM/DL7OhwfO5wYOl5rMPMPM/O3MPO3APhvAPz6Wp38KwH8LQkNm/F/HcCXiehX1jyeXw/gx+rmLwL47jr/3Wa9w+FwfBj8eQDfSETfQEQXlGKSLw5j7H3n1wH472v+8xcBfFetWv4GlEK5P/cpXbfjLYZ3PnG8MahE7z9EqdZ7D8C/wsxfqtt+qpJHENG3oVT1PQPw4wD+DWZmIvrFAP5zAH8PgP8bwD/PzH+TiP5OAF8C8HkAGcBXAHwTM/9/n+brczgcrx+I6NcC+A8ARAA/ysy/k4h+O4AvMfMXiegewB8G8K0oSuF3MfNfqft+P4DfAGAH8G8x84/X9X8EwD8M4GsA/A0AP8DM/8mn+8ocbyqcGDocDofD4XA4AHgo2eFwOBwOh8NR4cTQ4XA4HA6HwwHAiaHD4XA4HA6Ho8KJocPhcDgcDocDgBNDh8PhcDgcDkeFE0OHw+FwOBwOBwAnhg6Hw+FwOByOiv8fCl9h341pzAUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 612x388.8 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "hp.mollview(m, title='Relative Intensity', rot=180, min=-1e-3, max=1e-3)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c19bb56b-33c3-4c8d-a6c4-d720d6e7ea46",
   "metadata": {},
   "source": [
    "## Energy Distributions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "12cc8a28-414c-436f-b846-93a04a158123",
   "metadata": {},
   "outputs": [],
   "source": [
    "with open(dist_file, 'r') as f:\n",
    "    e_dists = json.load(f)\n",
    "\n",
    "# Loaded dictionary contains bins and energy histograms, separated by energy bins used in the analysis\n",
    "bins = np.array(e_dists.pop('bins'))\n",
    "ebins = sorted(e_dists.keys())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "ffde0f59-54f1-4fe3-8430-e880b0b906f0",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfsAAAGHCAYAAAC3a/toAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzs3Xl8VPW5+PHPQ9gJYSeBIJsrURQhivbWFq1aXFj0akFt6+WWilVvrbfV0vbqRVur9rbaqq3autsiVPtrQUW7qKi1vtCgUTSIC5sJhB0hRAjL9/fHmRnOnMw+55w5Z+Z5v155mZk5M/NNiHny/X6f7/OIMQallFJKFa8OhR6AUkoppbylwV4ppZQqchrslVJKqSKnwV4ppZQqchrslVJKqSKnwV4ppZQqchrslQoJEZkgIu+5+Hp/E5FLIp/PFJHFLr72pSLyrFuvl+F71ojI2yKyU0Su8PO9w0RERopIS6HHofylwV65SkRabB8HROQz2+1LfHj/34tIm2McS71+33yJyE9EZG8kUO0UkRUicqeIVEWvMcYsNsYcneFrPZzuOmPMmcaYP+Q5dETkMBGJK9hhjHnEGHNWvq+dpe8DfzPG9DTG/CbfF7P9m9h/lja7ME5PRf5w228b80oRuSz6uDFmpTGmvJBjVP7TYK9cZYwpj34Aa4FJtvvaBRYR6ejBMH5qH4cxZpzbb+DRuP9gjOkJ9AP+HTgEqBORSjffREQ6iEgx/r8/DMhp5SPFv+cfHD9L/XMfXtbvnY9XbP8ffgW4XURGe/A+KiSK8X94FWCR2dJ8EXlcRHYCX43MxufYrjldRFbbbg8RkT+LyCYRWSUiV+b43oeJiBGRr4tIY+T1Ztse7yAiPxSRj0Vks4jME5E+jufOEJG1wN8i988QkbWR638Yed0JIlItIq0i0tv2+uNFpDndL3djTJsx5l3gQmA7cE2S78sPRWSdiOwQkfcj73sucB1wiX1VQ0T+KSI/FpHXgF3A0Mh9/2F76w4i8hsR+VRElovIqbb3ahSRCbbb9tWDlyP3RWeSJzi3BUTk8yJSF3nt10VkvO2xf4rIjSLyr8iqxnMi0jfyWHcRmSsiW0Rke+S57QKuiLwMnALcGxnDSBHpHfnZ2iQiq0XkByIiketnisjLkdWTrcD/pPo3SfB+HSM/D7NE5CMR2SYidzqumRn5d9kmIs+KyCGO514hIh8B70fuP0tEPoh8j+4SkVdF5D9EpGvkax9le+1BkZ+vfunGaoypAz4ARkWeG7cS48b3XwWfBntVCOcBc4FewPxUF4pIGfA08AZQDZwBXCsiX8rj/T8HHAZ8GbhRRA6P3P/fwDnAF4AhWEHxTsdzvwAcBZwj1kzpTmB6ZGwDgCoAY0wT8E+sgB31VeBxY8y+TAYZuW4hVhCLIyJHA7OAscaYCuAsYK0x5mngZxyckdpXNb4G/CdQATQmeMvPYQWe/sCPgT/b/1hJ4QuR8UZnv284xtofeAb4BdaqxZ3AIon8IRVxMXApUAn0wPq3AJgBdMf69+gHXAHsdg7AGPMF4DXg8sgYVgK/iTx3JHAa8A3g646vdznWv9ttGXydiZwNjAOOx/rD9fTI13wBcC0wJfL6S7B+5u0mAycAo0VkIPDHyHP6A6uAEyNf2+7IY1+1Pfdi4K/GmC3pBigiJwGHAqm2s/L6/qvg02CvCuGfxpinjDEHjDGfpbn2JKDCGPPTyIz3I+ABrACbzOzILCT68YDj8TnGmN3GmDexln2Pi9w/C/ihMaYp8gt2DvAViV/y/l9jTGtk3BcCfzHG/MsYs4f2s8NHiPyCjszmpwGPpfl6ndYBfRPcvw/oChwtIh2NMasiAS6VB40xy40xe5P8wbEeuCvy+FxgJdYfEfmaBLxnjHncGLPPGPP7yGufY7vmAWPMh8aYVuAJYEzk/r1Ywe8wY8x+Y0ydMSZtcpmIdMJavp5tjNkZ+d7cgfUHT9RaY8w9kddN9nN4seNn6e+Ox28xxnxqjFkNLLaNexbWdtKKyPf6J8CJIlJte+5PjTHbIu99LlBvjFlgjNkbGas9P+CRyFgkcvtrpP5Z+nxkvC1YfwQ9iPU9T8bV778KHg32qhA+yeLaYVhLzrFfuFjL1FUpnnOrMaa37eMb9geNMc22m61ANFlpKPCU7X2WAQYYmGTsg+23jTG7gG22x/8MHCciQ4GJwKbIHxjZqAa2Ou80xqwAvgvcBGwUa1sk1ffEOfZEGk18Z6w1WF9jvgZHXstuDdbXFpXs3+Rh4B/AH0WkSURulcz2uAcCZY73db5nJj+Hcx0/S2c4Hk827mHAr20/S5uBA1gz5ETv7/xZMthWX4wxr2L9gfd5ETkG62f1mRTj/mdkvOXAIGAs1s9KMm5//1XAaLBXheBstbgLa6kwyh60PgE+dPzC7WmMmeTBuBqBMxzv1dX+x4EjGK7H9stbRHoAfWzXtgJ/Ai4h/UysncgWxiTglUSPG2N+b4z5N2AEVmC7JfpQkpdM1+JyiOP2UKyVBUj9b5TudddhBT/nazeleV40f2GOMWYU8HmsLaBMTnVsBPY73tf5nl62/PwE+IbjZ6mbMWZJkvd3/iwJ8X+YADyKtVL0NeCPkdWktCI/v/8P62cpK3l8/1XAaLBXQVCPtQfeR0QGAd+2PfYa0CYi340kKpWJyGgRcT3DHrgX+GlkJo6IDBSRySmufwKYKiIniUhnEs+cHsXaJz8H+H0mgxCRTiJSA8zDWsL/ZYJrRonIqSLSBfgs8rE/8vAGYLhtyTdTg0TkqkgC2XSsfd7nIo/VA9Mjj50InG973kbAiMjIJK/7NNZ2w7TI8y/GyplYlG5AInKaiBwT2UrZgbWsvD/N04gshT+J9e9ZLiIjsBIdM/o3cMG9wI+iSXViJQtekOL6p4GxIjIpMnO+Gmuv3+4x4AKs/fVHMx1IJGdiKjmcVMj1+6+CR4O9CoKHsRKl1mAFl3nRByL7nWdjJSutxloOvQ8rySyZH0r82ejmFNfa3R55/+fFOinwL6wEqoSMMe9gBZAnsGavWyIf9hnXy1iz7iXGmERJcXaXRN53G7AAK2jXOrYdorpgJeJtxlqC7cPBnIH5QGdgq4i8nuY97f4FHI21bTAH+HdjTHRb4kdYiYnbgeuxJZsZY3ZirSosiSxb19pf1BizCSsZ7ftY359rgHONMe22JxIYjDUr3YEVrP4BPJ7h13MF0IaV7PYS1r53xkEy4hLHz1JLhhnwT2D9PD0hIjuAd7ASQpNdvwErp+N2rO/RocBb2H6WInkBy4A2Y8y/0gzhlOh4gQaslYOr0407gXy+/ypAJH5VUimVKxGpwAqGw4wxn9jufxkrOe7hQo1NhUtkC2cdcIEx5hXb/Y8CK40xcwo1NhVOOrNXKg8iMjlyFrkc62jZm45AfxJwDNbsX6mkRGSiiPSKbM1cj5WQ97rt8ZFYR/keLNAQVYhpsFcqP+dhzcAageHARdEHROQPWNsCV0cy9ZVK5fNYx+M2Y53emBpNwhORW4C3sY7rrS3cEFVY6TK+UkopVeR0Zq+UUkoVOQ32SimlVJErmkpI/fv3N8OHDy/0MJRSSinfLF26dLMxxlmToZ2iCfbDhw+nrq6u0MNQSimlfCMizlLUCekyvlJKKVXkNNgrpZRSRU6DvVJKKVXkimbPXimlVDjs3buXxsZGdu/eXeihhEbXrl0ZMmQInTp1yun5GuyVUkr5qrGxkZ49ezJ8+HCyb85YeowxbNmyhcbGRkaMGJHTa+gyvlJKKV/t3r2bfv36aaDPkIjQr1+/vFZCNNgrpZTynQb67OT7/dJgr5RSqiTt37+f448/nnPPPTfpNTt27KC6upqrrroq4ePXXnstRx11FMceeyznnXce27dvB2D16tV069aNMWPGMGbMGC6//HJPvoZMabBXSilVkn71q18xatSolNdcf/31fPGLX0z6+BlnnMG7777LO++8wxFHHMEtt9wSe+zQQw+lvr6e+vp67r33XtfGnQsN9koppUpOY2MjzzzzDDNnzkx6zdKlS9mwYQNnnnlm0mvOPPNMOna0ct1POukkGhsbXR+rGzQbXymlVMHc+NR7NKzb4epr1gyu4H8nHZ3ymu985zv87Gc/Y+fOnQkfP3DgAN/97nd57LHHeP755zN63wcffJBp06bFbq9atYrjjz+eiooKfvKTn3DKKadk/kW4TIO9UkoVmaamx2ne8FTcfVWVk6iuvqhAIwqWp59+moEDBzJu3DgWL16c8Jrf/OY3nH322RxyyCEZvebNN99Mx44dueSSSwAYNGgQa9eupV+/fixdupSpU6fy3nvvUVFR4daXkRUN9kopVWSaNzxFS0sD5eU1AGzfvoTt25fE/QEQlOCfbgbuhVdffZWFCxeyaNEidu/ezY4dOzjrrLNYv349ADfddBOvvfYar7zyCr/5zW9oaWmhra2N8vJybr311nav98gjj/D000/z/PPPx7Lmu3TpQpcuXQAYN24chx56KB988AG1tbX+faE2GuyVUqoIlZfXMG7sXKD9TL+lpYFmCESwL4Rbbrkllki3ePFifv7zn/P000/HXTN58uTY5w8//DB1dXUJA/1zzz3HbbfdxksvvUT37t1j92/atIm+fftSVlbGypUr+fDDDxk5cqRHX1F6GuyVUirkEgXz6KwerKBuD+xL37zY1/EVm5kzZ3L55ZdTW1vLVVddxZ49ezjjjDMAK0nv3nvv5eWXX+aGG26gY8eOlJWVce+999K3b9+CjVmMMQV7czfV1tYa7WevlCpFS9+8uF2AT7VMHw320Zm/35YvX572yJtqL9H3TUSWGmPS7g3ozF4ppYqAfdleKSc9Z6+UUkoVOZ3ZK6VUCWppaYjbuw9Kdr7yhs7slVKqxFRVTorb329paWh3Ll8VF53ZK6VUidHs/NKjM3ullFKqyGmwV0qpkGlqepylb14c+2hpaSj0kEIpVYvb22+/nZqaGo499li+9KUvsWbNmthjZWVlsda19uI7iaRrkTtnzhyqq6tjr7do0aL8vqgkdBlfKaVCxlkOt7y8hqrKSQUeVfhEW9zu2NG+Ec/xxx9PXV0d3bt355577uG6665j/vz5AHTr1o36+vqM3iNdi1yAa665hu9973vZfwFZ0Jm9UkqFUPRcffRDM+mzk67F7amnnhorf5tr69pMWuT6RWf2SimlCufZ2dC8zN3XrBoNZ7WvY2+XrsWt3QMPPMBZZ50Vu717925qa2vp2LEjs2fPZurUqe2ek02L3LvvvptHH32U2tpafvGLX9CnT5+0Y8qWzuyVUkrFzt1HP5qaHi/0kDxjb3Gbzu9//3vq6uq49tprY/etXbuWuro65s6dy3e+8x0+/vjjds/LtEXut771LT7++GPq6+sZNGgQ3/3ud7P/gjKgM3ullCpxVZWTaLbd9rUrXpoZuBcStbj96le/yu9///u46/7xj39w880389JLL8Xa1QIMHjwYgJEjRzJhwgTeeustNm/ezKxZs4DsWuRWVlbGPv/mN7+ZMFnQFcYYzz6AicAK4CNgdoLHuwDzI48vAYY7Hh8KtADfS/de48aNM0opVQrqll5k6pZeFNrXb2ho8Oy1s/Xiiy+ac845p939b775phk5cqT54IMP4u7funWr2b17tzHGmE2bNpnDDjvMvPfeeynf46GHHjJXXnllwsfWrVsX+/z2228306ZNS/o6ib5vQJ3JIB57towvImXAr4GzgBrgIhGpcVz2DWCbMeYw4A7gNsfjdwDPejVGpZRSKuqGG25g4cKFAFx77bW0tLRw4YUXxh2xW758ObW1tRx33HGceuqpzJ49m5oaZ2hLbebMmUS7tF533XWMHj2aY489lhdffJE77rjD3S8qwrMWtyJyMjDHGPPlyO0fABhjbrFd89fINa+JSEegGRhgjDEiMhX4N2AX0GKM+Xmq99MWt0qpUuF1i1qvX19b3OYmnxa3XiboVQOf2G43Ru5LeI0xZh/wKdBPRHoA3wduTPUGInKZiNSJSN2mTZtcG7hSSilVTLwM9pLgPucyQrJrbgTuMMa0pHoDY8xvjTG1xpjaAQMG5DhMpZRSqrh5mY3fCNjPHAwB1iW5pjGyjN8L2AqMBy4QkZ8BvYEDIrLbGHO3h+NVSimlipKXwf4N4HARGQE0AdMBZ2ulhcClwGvABcALkezCU6IXiMgcrD17DfRKKaVUDjwL9saYfSJyFfBXoAx40BjznojchHVUYCHwAPCYiHyENaOf7tV4lFJKqVLlaVEdY8wiYJHjvhtsn+8GLkzzGnM8GZxSSoVEU9PjNG94Knbb3gRHqUxouVyllAq4aJe7KO1y545ULW4ffvhhBgwYEGs9e//996d8rSeffBIRIdER8BUrVsReZ8yYMVRUVPDLX/4S0Ba3SimlbKJd7pR7UrW4BZg2bRp3350+XWznzp3ceeedjB8/PuHjRx55ZKwl7v79+6murua8886LPa4tbpVSSikPpGtxm43rr7+e6667jq5du6a99vnnn+fQQw9l2LBheb9vNnRmr5RSqmBue/023t/6vquveVTfo/j+id9PeU0mLW7/9Kc/8fLLL3PEEUdwxx13JOxg99Zbb/HJJ59w7rnn8vOfpyz0CsC8efO46KL4BkPa4lYppZRyWSYtbidNmsTq1at55513OP3007n00kvbXXPgwAGuueYafvGLX2T0vm1tbSxcuJALLzyYl64tbpVSShVMtL89WC1wvWp3m24G7oVMWtz269cv9vk3v/lNvv99a5w/+tGPeOaZZwB46aWXePfdd5kwYQIAzc3NTJ48mYULF1Jb275c/bPPPsvYsWPj2tr61eJWZ/ZKKaXiVFVOih3ta2lpiDv2VwxuueUWGhsbWb16NfPmzeO0005r18t+/fr1sc8XLlwYa0Bz8803U19fT319Pb169WLz5s2sXr2a1atXc9JJJyUN9ACPP/54uyV8+/v8+c9/5phjjnHry4yjM3ullFJxqqsvis3ko7P7UnDDDTdQW1vL5MmTufPOO1m4cCEdO3akb9++PPzww1m91rp165g5c2bsKF1rayt///vfue++++Kuu+6666ivr0dEGD58eLvH3eJZi1u/aYtbpVSx8rrlrN/vrS1ucxPUFrdKKaWUCgAN9koppVSR02CvlFJKFTkN9koppVSR02CvlFJKFTkN9koppVSR03P2SimVg7q6OpYtWxZ33+jRo5MWVMmG9q/3x/79+6mtraW6upqnn3467rFrrrmGF198EbDOyG/cuJHt27cDUFZWxujRowEYOnQoCxcuTPj6mV7nBw32SimVg2XLltHc3ExVVRUAa9asYc2aNXF/AOQa/KP966MBXvvXeyNVi9s77rgj9vldd93FW2+9FbvdrVu3WMvaVDK9zg8a7JVSKkdVVVXMmDEDaD/Tb25uBsh5pq/9670VbXH7ox/9iNtvvz3ltY8//jg33nijTyPzhgZ7pZTKQKJgHp3VgxXU7YH9oYce8nV8YdX805+yZ7m7LW67jDqKqh/+MOU1mbS4BWvFZtWqVZx22mmx+3bv3k1tbS0dO3Zk9uzZTJ06NeFzM73ODxrslVIqA85l+6qqqth+rAoXe4vbxYsXp7x23rx5XHDBBZSVlcXuW7t2LYMHD2blypWcdtppjB49mkMPPbTdczO9zg8a7JVSKkP2ZftMNDc3x83w3UrgKybpZuBeyKTFbdS8efP49a9/HXff4MGDARg5ciQTJkzgrbfeYvPmzcyaNQuAm266icmTJye8rlDBXo/eKaWUB0aPHh23zN/c3Nwuez8sor3tox9NTY8Xekh5yaTFLcCKFSvYtm0bJ598cuy+bdu2sWfPHgA2b97Mq6++Sk1NDePHj4+1vp08eXLS6wpFZ/YqMOYuWcuC+qa4+6aMqebi8UMLNCKlclcse/hVlZNott1uaWmgGWItcIuJvcUtWIl506dPR0Ri1yxfvpxZs2bRoUMHDhw4wOzZsxMG8Uyv84sGe1UwzuC+ZNVWAMaP6AtAw3rrOIwGe6UKx97bHoqvv/2ECROYMGECYC2/282ZM6fd9Z/73OcyWqHJ9Dq/aLBXBbOgvomG9TuoGVQBWEHePpOfdt9rNKzfwbT7Xos9R2f6SimVPQ32yjfOmXw00M+fdXLC66eMqY67rTN9pZTKjQZ75RvnTL5mUEW7gG538fihcYHdPsNXymvpztXnQrPzVaFosFe+SjWTz4Qu6yu/uH2u3vncfCvsKZUNDfYqNHRZX/kt23P1qRRLdr4KJw32ylP2fXr7En4udFlfKaVyo8Feecq+T59uj14ppfw0fPhwevbsSVlZGR07dqSuri6na/7whz9w2223AVBeXs4999zDcccdl/Hz/aDBXnku3336VHQPXymVjxdffJH+/fvndc2IESN46aWX6NOnD88++yyXXXYZS5Ysyeo9vKbBXoWW7uErpYLgc5/7XOzzk046icbGxgKOJjEN9iq0dA9fqfB75Y8fsPmTFldfs/8h5ZzylSPSXicinHnmmYgIs2bN4rLLLsvpGrsHHniAs846K+fne0WDvVJKqZL06quvMnjwYDZu3MgZZ5zBUUcdxRe+8IWsr4l68cUXeeCBB/jnP/+Z0/O9pMFeuSpZlTylgs6LIjrp2IvslGqBnUxm4F6JtqAdOHAg5513Hi+99BLf/va3Abj88su5/PLL213z+uuvJwzW77zzDjNnzuTZZ5+lX79+Sd8j2fO9pi1ulaui2fdRmoGvwiJaRCcq3yI66dhb4O7fv5gNG6+NtZBtaWnw7H2VZdeuXezcuTP2+d/+9jdOOOGEWJvayy+/POE1xxxzTLvXWrt2Leeffz6PPfYYRxxxRMr3SPR8P+jMXrnOy+z7dDQ7X+XDzSI66diL7Dyz6FE6d94ADAKgvLyGqspJvoyjVG3YsIHzzjsPgH379nHxxRczceLEjK+59957AWsF4KabbmLLli1cccUVALEjdpm8h1802Kuiodn5Ksza2ioZN3ZuoYdRMkaOHMnbb7+d8zWXX3557PP777+f+++/P6f38IsGe1U0Li57nos7Pxm7/V7nT2EL8FCvgxeNvgBq/Zm5KVWsWloa4vraV1VOiut5r4JHg70qHsuehOZlUJVkn3XNP62PZQf/INDgr1R2qion0Wy73dLSQDNosA84DfYqvOoeig/c0UA/4xkAbors3c+fcXLy60GDvVJZqK6+KC6w22f4Krg02Ku8FPSonXMmXzXamqknUzsjPrA/dI71/IfOsW7rLF8pVaQ02Ku82BvdQAGO2tlm8lmz/2Ggs3ylVBHTYK/y5ttRu2TL9imkPIpnn+lHZ/dKKVWENNir8Mhy2V6P4imlUknXftaN1rWZXLd48WKmTJnCiBEjADj//PO54YYb3PxSNdirkMli2T7rRjn2/XvQPXylSkCq9rNuta7N5LpTTjmFp59+OrvBZ0GDvVLQfoVA9/CLXiFq4afS1tYWq5MPpVsrP0jC0Lo2UxrsVXDlsEefs0SZ+qqoRWvhRwO817XwU+nRoxw42OY1WqO/FIL9iw//lo1rVrr6mgOHjeTU/0jfSjab9rO5tq7N9LrXXnuN4447jsGDB/Pzn/+co48+Ou34s6HBXgVXtkfrlMqSn7XwU+nZsyc9e/bknLOtsdhn+Mo7mbafzad1bSbXjR07ljVr1lBeXs6iRYuYOnUqH374oatfqwZ7FWz5HK3Ll+7hK+W5TGbgXsmk/Ww2rWtHjBjBpElWA6NsWuRWVBysTXL22WdzxRVXsHnz5ozyATKlwV5lpWT61esevlJFbdeuXRw4cICePXvG2s86M+BTta5N9NxDDjmE+vr6rN4DrG2byspKRITXX3+dAwcOxP1h4QYN9iorBS+ik6eMW+DqHr5SRS1Z+1k3W9dm2iL3ySef5J577qFjx45069aNefPmISKufr0a7FXWCtmvPh967l4pFZWs/aybrWszbZF71VVXcdVVV2Uy7JxpsFfB4XH2fdbn7pVSqkh0KPQAlIqJZt9Hafa9Ukq5Qmf2KlgKmX2fjmbnK6VCSoO9UpnQ7HylkmppaYjra19VOSmu570qPA32SmVCs/OVSqiqchLNttstLQ00gwb7gNFgr5RSKmfV1RfFBXb7DF8Fh6cJeiIyUURWiMhHIjI7weNdRGR+5PElIjI8cv+JIlIf+XhbRM7zcpxKKeW3pqbHWfrmxSx982JaWhoKPZySNHz4cEaPHs2YMWMS9iFYvHgxvXr1YsyYMYwZM4abbrop6WstXryYMWPGcPTRR/PFL37Ry2HnxLOZvYiUAb8GzgAagTdEZKExxv5T/Q1gmzHmMBGZDtwGTAPeBWqNMftEZBDwtog8ZYzZ59V4VWKeVszzs9FNEvYiO0kL7KiiELQud80bnqKlpYHy8hrKy2uoqpxUsLGUsnTtZzNpPbt9+3auuOIKnnvuOYYOHcrGjRvdHmbevFzGPxH4yBizEkBE5gFTAHuwnwLMiXz+JHC3iIgxptV2TVfAeDhOlYKnFfMK3OjG/nXkVGBHs/NDJUhd7qLKy2sYN3Zuwseam5u15W1IzJ07l/PPP5+hQ63fHwMHDizwiNrzMthXA5/YbjcC45NdE5nFfwr0AzaLyHjgQWAY8LVEs3oRuQy4DIh9k5X7PK2YV8CjdvYiO1kX2NHs/FAKSpe7dJx/hBRzy9vtT31M27pdrr5m58E96D3p0LTXZdJ+NpPWsx988AF79+5lwoQJ7Ny5k6uvvpqvf/3rrnwtbvEy2Ccq7OucoSe9xhizBDhaREYBj4jIs8aY3XEXGvNb4LcAtbW1OvtX/tHsfOWh2trauMCuLW+9ka79bKatZ/ft28fSpUt5/vnn+eyzzzj55JM56aST4prnFJqXwb4ROMR2ewiwLsk1jSLSEegFbLVfYIxZLiK7gGOAOu+Gq5RSym+ZzMC9kq79bLLWs/Pnz+d3v/sdAIsWLWLIkCH079+fHj160KNHD77whS/w9ttvByrYe5mN/wZwuIiMEJHOwHRgoeOahcClkc8vAF4wxpjIczoCiMgw4EhgtYdjVUopVUJ27drFzp07Y5//7W9/45hjjom7prm5GWOsRWN769krr7yBDjdpAAAgAElEQVSS+vp66uvrGTx4MFOmTOGVV15h3759tLa2smTJEkaNGuX715SKZzP7yB78VcBfgTLgQWPMeyJyE1BnjFkIPAA8JiIfYc3op0ee/nlgtojsBQ4AVxhjNns1VqXSeeKDJ1i0clHcfWePPJsLj7iwQCNSSuUjkxa3mbaeHTVqFBMnTuTYY4+lQ4cOzJw5s90fDoXmaVEdY8wiYJHjvhtsn+8G2v22NMY8Bjzm5diUSsUZ3Os2WDtItZW1sdt1G+oOXiMbONv0aP/DrJQKpExa3GbTevbaa6/l2muvdW18btMKekoB28pe5tOy15nxnLVH5wzutZW1cTP5dn8MyB7qZA+LnjuYtKczf6VUUGiwV/4JQBGdZD4te51Ws5aG9dZRvO4cwRcGn8n/fXlWwusvPOLCuED+xENfZNHejbEjeHWyJ37mjwZ/pVThaLBX/ilwEZ1U+vfowuZdQxne9j3AKrKzti3zSoEXjv4PLlz2ZOxw6ROfvs+iXn1ij6/YusK6ToO9UqoANNgrfwWoX719KX7bvtXUDDqShyZaxYOyLrLjOHd/4UPncKEBJlrno2c8F/xCLkqp4qXBXsXxtBZ+wCxauYgVW1dwZN8jObLvkZw98mxP32/F1hVxQV+X9ZVSftFgr+J4Wgs/gI7seyQPTfS+OpnzDwld1ldK+UmDvWrH01r4BeTMoI/O6v3gTOjTZX2VrTA1xmlpaYjra19VOSmu530QbN++nZkzZ/Luu+8iIjz44IOcfPLJWV+zePFipkyZwogRIwA4//zzueGGGwgaDfaqZNiX7QFflu5T0WX97LUsWU9r/aaU13QfM4Dy8YOA+La2hW5pm48wNcapqpxEs+12S0sDzRC4YH/11VczceJEnnzySdra2mhtbc3pGsisDW6habBXJcWvZft0dFk/N631m9i7voVOg8oTPr53fQutEAv29ra2QWhpm6swNcaprr4oLrDbZ/hBsWPHDl5++WUefvhhADp37kznzp2zviZMNNirolXIZXsgZb97XdbPjHMmHw30A2cdm/D6jfe9w971LWy87x3r+i27GFDeNxRtbUvVs88+G1upcEtVVRVnnXVW0sdXrlzJgAEDmDFjBm+//Tbjxo3jV7/6FT169MjqmqhM2uAWmpeNcJQqqOiyfZSvy/ajL4gvGNS8LL6gkEqoZcl6Nt73Tuxj+58/om3Vp7HHOw0qp/uYAUmf333MgLhZv2nbz4GWvZ6OWYXPvn37ePPNN/nWt77FW2+9RY8ePbj11luzvgYOtsF9++23+a//+i+mTp3q15eRlbQzexG5EHjOGLNTRP4HGAv8xBjzpuejU+Fnr5pXgIp5+SzbN6zfEXfefsqYai4ePzSzJ+fQ71738Nsv03ce0StuDz6d8vGD4q6Vn76Madsfm+ln81rKH6lm4F4ZMmQIQ4YMYfz48QBccMEFzJkzhzFjxgBWffypU6e2uyZRsE/WBrd///4+fCWZy2QZ/3pjzBMi8nngy8DPgXuA8Z6OTBUHe9W8AFXMS8d53LBh/Q6AzIN9lnQP/6BUy/TZ6lDeiQMt1ufO/XxVuqqqqjjkkENYsWIFRx55JM8//zxjx47l//7v/+Kuc15TU1PT7rWam5uprKxEROLa4AZNJsF+f+S/5wD3GGMWiMgc74akik6AquZl6uLxQ+MCe9YV9bKke/jeKOvZmbKenRk449jY7F4pgLvuuotLLrmEtrY2Ro4cmTDpMdk1ubTBLbRMgn2TiNwHnA7cJiJd0L1+FUAFT8hTWUuWgOcVe/Ie6LJ+KRszZgx1dXU5XZNrG9xCyiTYfwWYCPzcGLNdRAYBwW3aq0pW0M7Rq/Sce/TpEvDy0X3MAOynpP1e1m9qepzmDU/Fbre0NFBe3n5ZWCkvZBLs7zPGfC16wxizXkR+BvzNu2EplZugnKN3Q6kk7Lm5R5+KM3nP72X95g1PxQX48vIaqion+ToGVboyCfZxBwZFpAwY581wlCpiKc7dOxVrwp7fy/ZBU15ew7ixcws9DFWCkgZ7EfkB8EOgm4jsiN4NtAG/9WFsygel1OWuoJynEJqtEq7Jgn2xJuz5uWyvlDooabA3xtwC3CIitxhjfuDjmJSPSq3LXaFs+7gHO16wHcdp7gesgxe+Hrur4txz6TPtK/4Pzmd+LdtnQhP2VKlIu4xvjPmBiFQDw+zXG2Ne9nJgyj9h7XIX5Oz7bfP/yA5bY4zWN94AoPsJJyS8vvWNN2h9442455RK8C+UQifsKeWnTCro3QpMBxo4eObeABrsVUEFKfs+XXDvfsIJ8cE7unc/49Gkz7cH/+lb3+fFo/YzA2s5v1iT9fxU6IQ9VVjp2tdm2rr2/fffZ8aMGbz55pvcfPPNfO9734s99txzz3H11Vezf/9+Zs6cyezZswFYtWoV06dPZ+vWrYwdO5bHHnvM8yY7mSTonQccaYzZ4+lIVHGwl8cFz0vk+pl9n6p87o6nn2b3++/T9aijgATBPY0+074Sd60z+A/9eCeXfgxr33+P1n2fsWb8Brgp+MHe74Q8e0tbCHdb23TC1N8+iDJpX5tJ69q+ffty55138pe//CXu/v3793PllVfy97//nSFDhnDCCScwefJkampq+P73v88111zD9OnTufzyy3nggQf41re+5erX55RJsF8JdAI02Kv07OVxIVQlclNx5jEMeuU5+jzxNmsiuQ7RQD/ssUddeb9kwf8oYOuyN+m+dLMr7+M1vxPy7C1tgVC3tU0lTP3tg8jN9rUDBw5k4MCBPPNMfJXQ119/ncMOO4yRI0cCMH36dBYsWMCoUaN44YUXmDvXOpVx6aWXMmfOnEAE+1agXkSexxbwjTHf9mxUKtxCWB43HWf53GefuI2qTWth0DEAdD3qKCrOPdez97cH//fPOZGBTa2s+Vo4kvv8TsirqqrKuaVtWBL2wtTfHqwCQva+9p07/Xfs8w8++DE7W5a7+n49y0dxxBHXJ3080/a1+bSubWpq4pBDDondHjJkCEuWLGHLli307t2bjh07xu5vampK9jKuySTYL4x8KFWynMvqVZvW0jxgKMe7NJPPxvJx/YHN9I3c3v3++wCBDfZhoQl73qiqnIS9W31LSwO9e31WsPHAwfa1d911F+PHj+fqq6/m1ltv5cc//nHsmmjr2vLychYtWsTUqVP58MMPM34PY0y7+0Qk6f1eyyQb/xER6QYMNcasSHe9Ul4pZPa9c0++ecBQlh05Hv+bc8Lbn6vk7c9VxnIV7DN8P7zzj+dY/urijK49Zs9JdK/o7e2AXKIJe96orr6I6uqLYrftM3wg5QzcK4la3Drb1yZrXTt//nx+97vfAbBo0SIGDx6c9D0++eST2O3GxkYGDx5M//792b59O/v27aNjx46x+72WSTb+JKy2tp2BESIyBrjJGDPZ68EpZVfo7Hv7nvx1bnTBy6KinpO9lO70re9T3bzXs2V9Z3BvbHgXgCE1x7S7tnLfIQzYfzC/odveHmxb38TrN86O3Tfq3yZw7OkTXRmbUrlI1OLW2b42WevaK6+8kiuvvDLte5xwwgl8+OGHrFq1iurqaubNm8fcuXMREU499VSefPJJpk+fziOPPMKUKVO8+lJjMlnGnwOcCCwGMMbUi8gID8ekVFJ+Zd87l+3ts3pXZFlRz875B86LR+3nVDoRHd1HmyvY+I99dG18M3bNESdWcvQpmRVLShfch9QckzRgb7zvnbiEvF3btrJj1/a412pseDfu9TX4q0JI1L42l9a1zc3N1NbWsmPHDjp06MAvf/lLGhoaqKio4O677+bLX/4y+/fv5z//8z9je/633XYb06dP53/+5384/vjj+cY3vuH515tJsN9njPnU8UW233RQqog4l+1dT8CrnREf2O0z/DRqNnyOju8dGru9ouv7vDUGWvtaY133oRVc+0b28reWVbHuw+188PqGpK9p/2Ng+auL2bR6FQOGW3/TpwruiTgT8kYwIfa58w+JTatXAWiwV75L1L42l9a1VVVVNDY2Jnzs7LPP5uyz268+jhw5ktdffz3LEecnk2D/rohcDJSJyOHAt4F/eTsspQrPzaN0bvrg9Q1sbmyh/5DE59UH9t5L5YY3GLbvAyAy0x/yOSDx3vnad19m1Zvvs/gR69dBW2szvSqHMu1/b014fT6OPX1iXGCff+NsNq1exXxd5lfKU5kE+/8CfoR17G4u8FfgJ14OSikV771XmmIz82igP++7YwGY8dxdAMyeaE98+vLBT7/2dQ7b/RzDvhufGBX14DW/59MNmwHrbLqR/uzcPpQ//yK3bYBsjPq3CXG3gzbTD8tRPKXSySTYH2mM+RFWwFeqKHm+R58le3CHg0vzgw/vTf8h5RxxYmVWr7f7/fdjCXwr9+9mXZ+edBxoFbdp3d7E4CMPi83kne+9ubEFIGmwt1fJy7ZCXpBn+noUTxWTTIL97SIyCHgCmGeMec/jMSkPaUvbxHLZo09VPjdfzqX6wYf3Tjm7tmfnQ3ztfOfXsbplOzv3f0ZlJNgPGD4iboZ99CnVce/z51+8yebGlqQzfXuVvHwr5AVppq9H8bxjMBhjfDlfXiwSnc/PRibn7E8VkSrgK8BvRaQCmG+M0aX8EHK9pa2HtfD9PlefzR6983vWsH4HQM7B/r2No/lg8yiIBFTnUn0qzuz8FVutchjRYP9JvwqWH3bwHO/Ohm30krKM9+SdqwiJZvpuVclLNNNXxceYTWzZsoV+/fppwM+AMYYtW7bQtWvXnF8jk5k9xphm4E4ReRG4DrgB3bcPLVdb2npYC7/Q5+pTcZbPnZbnufsPNo9ic+sA+kduZ7NUf+ERF8Z1wLPP8KF9dn0vKWNoh8zrgKeb6R+5/TO69/SuY1dQlvWVew7sX8TOnWPZtGlT+osVAF27dmXIkCE5Pz+TojqjgGnABcAWYB7w3ZzfURUfD2vhe3mu3r5PX+g9eoD+HVdzXt9fH7yj2wVAbjXenQYMHxGbya/52tfj9vAhuyI8zj9C9u7ZTyttrozTKd2yfil1uSsq8lmsdazyRyYz+4eAx4EzjTHrPB6PUr6x79N73cjGqV0S3J5q+nexNcPIosiOU58Vn9F75R7mL7Fmw/ZZPbTfw8+2tv6wzh0YUH7wV0frZqF9c1B3pFvWL5Uud5nQlrcqlaTBXkQqjDE7jDEnJXhsqDFmrbdDU8p7hTpL70zA6z+sL0ecOApOudS6IIsiO069V+6h69Z9RDvlOBPwnO1zs62t72xb29qxA2tb97Misqzv1TG9KPuy/ibTmZ4VvXLucuelpqbHad7wVOx2S0sD5eU1KZ6RO215q9JJNbNfDIwFEJHnjTFfsj32l+hjSqncZJqAl4vdfTt6UhQnyp6Qt+mVJj6z1QCA5Mf08uVc1m/b7W/3tGzO3TdveCouwJeX11BVOcmTcYWt5a3yX6pgb0+R7JviMaVUGonOriergJctZwnarlv3sbVnW9KjeInks4dvT+BLd0wvX85l/Z/NucGV181ELufuy8trGDd2rudjUyqdVMHeJPk80W2lAq+QhXPaLdvnUBgnGWe2fffBA2kevCv2uPMonlO+e/h2mRzTCys9d++elpaGuFa3VZWT4trgKvelCvYDReS/sWbx0c+J3M69aoZSSXh9rt7z5jZpeLlsb8+2d3IexXPKdw/fLtExPa+17f5Mj+aFSFXlJJptt1taGmgGDfYeSxXsfwf0TPA5wP2ejUiVLD/O1fuVkOflsn3YeLms372iV+Qz6+hf0Grrq/aqqy+KC+z2Gb7yTtJgb4y50c+BKAX+9av3gr187rEf7qFPGwwaZlUqdHPZPlGbWPvROrfZa99DdvXvvV7W79GnLz369GVaJBtfK+4plVhGFfSUUqk5y+e27tkHXTpyhQfL9s49eufROjfYE/Y69D8H6V5F56F9ALKqf1+IZX2lVHsa7JVygbN87s3fW5z/izYviz9vP/qCWJGdVHv0+XLmMRxobaUDzQyc9UVP3i/MtAWuCotURXWuNsb8SkT+zRjzqp+DUsoNfmbfO/foyz87QEu3Drm/oLO/QB4V9SB1VzwnZ8LeJ9fMz+k9k/FyDx/8q6WvLXBVmKSa2c8AfgXchRbQUVEedrlzm5/Z986jdS3dOrCxT1nuL1g7Iy6wv3PTOSxfAqydnfUefbqueH7yeg/fzxa5ehRPhUmqYL9cRFYDA0TE/lMsgDHG5N/PUoWPh13uvOBnOVz70bp5eXbBc1q+DjbthAGDst+jT9cVz8mZkEenfhz49JOci+7Y5buHn67xjbbIVSqxVNn4F0X62P8VmOzfkJSb5i5Zy4L6gw1W7L3sc+ZhlzuV3ICeeFoCN8pZ+76sAsyn62OVtPIpupNINsv62vhGqdykTNCL9LE/TkQ6A0dE7l5hjNnr+ciUKxbUN8UF+JpBFe0yxwvF6yI6Knf22veW02Of5VN0xymXZf2qqqqsGt/Y9/C14I4qVZn0s/8i8CiwGmsJ/xARudQY87LHY1MuqRlUwfxZJxd6GO34UUTHK1o0xx1eH82zb3eUWsEdbXmr7DI5enc7Vi/7FQAicgRWf/txXg5MlQY3i+j4mX3vZa17SFA4Z6e1jK+yY9/DL6X9e215q5wyCfadooEewBjzgYh08nBMSuXE79r3Xta6b1c4pyeMGuze62dzFC+RfLrkKe9py1vllEmwrxORB4DHIrcvAZZ6NySlcudn9n069vK5YFXZsxfeSSeucI69uE6enFslQz7sRdVbho19rUM36crhutklLxGvz+F7yV5kZ291Cx3KdV6kgiGTYP8t4Erg21h79i8Dv/FyUEqFnTMJsmH9DoCsgn079op6tmp62XIexXvlrT9S1dIX+lq305XDdbNLnlOYW+Q6i+wcaNsPLQUbjlJx0gZ7Y8werH37270fjlLBlG1CnrN87rR8z93baxnkWU0vkebyrRw+6/T0F3rM64Q9L6vrOYvsrHw6j6JKSrlMa+MrlQHfE/KcVfLsFfVcXNIvJX5W11PZaWlpiGt1W1U5Sfvbu0yDvVIZ8jUhz4NOdl6yJ+y5naxn38Pf3NpC956dc3odr6vrNTU9TvOGp2K3P+uyim57vGs9XCyqKifRbLvd0tJAM2iwd5kGexVafh6184OXney8ZE/YcztZz7l6snfPflppc+W13da84SlaWhooL68BoNueEfTe+YUCjyr4qqsvigvs9hm+ck8mRXWOAK4FhtmvN8aclsFzJ2I10ykD7jfG3Op4vAtWwZ5xwBZgmjFmtYicAdwKdAbagGuNMS9k+kWpYHK7Yp7fR+1UYvaEPTeT9aD9Hv5HPw52A87y8hrGjZ0LwMY3tDGOCo5MZvZPAPcCvwP2Z/rCIlIG/Bo4A2gE3hCRhcaYBttl3wC2GWMOE5HpwG3ANGAzMMkYs05EjsGqzx/8dNxi5GKXOy8q5nl11K6YK+Q5G91UtfRlZdfGvM7dK6WCLZNgv88Yc08Or30i8JExZiWAiMwDpgD2YD8FmBP5/EngbhERY8xbtmveA7qKSJfIyQDlJ5e73LlZMc9LBU/I85Cz0U3bAOGjioO7poVsgevk7HL32f5P6VbWy7XX9zI7X6kgySTYPyUiVwB/BmLB1hizNc3zqoFPbLcbgfHJrjHG7BORT4F+WDP7qH8H3koU6EXkMuAygKFD8zi/rFIr0S53xZyQZ290MxA4mi/FHkvXAtdPzi533cp60Wlnf1eK7viRnW8vsgPWWXz78Tyl/JJJsL808t9rbfcZYGSa50mC+0w214jI0VhL+2cmegNjzG+B3wLU1tY6X1upQAtrQl46bpfStXe5S7S9ArkV3fE6O99ZZGfv+hZaQYO9KohMiurkurbYCBxiuz0EWJfkmkYR6Qj0ArYCiMgQrNWErxtjPs5xDEoVJ3s1Pcirop6bvC6l63XRHTc5i+zYZ/hK+S2TbPxOWCVzo2dIFgP3ZdDT/g3gcBEZATQB0wHnmYqFWCsHrwEXAC8YY4yI9AaeAX5gjAl2+q3yTdiP2tlr5R+1bgf9y7vk9kLOnAkPKurlystSukqp3GWyjH8P0ImD9fC/FrlvZqonRfbgr8LKpC8DHjTGvCciNwF1xpiFwAPAYyLyEdaMfnrk6VcBhwHXi8j1kfvONMZszPxLU8XGy6N2XmffO2vlt7btZ3NLjvmm9mp6kLainjP7Pl2jG1WctL99acsk2J9gjDnOdvsFEXk7kxc3xiwCFjnuu8H2+W6gXcqvMeYnwE8yeQ8Vb+6StSyob4rdbli/g5pBFQUckbu8Omrndfb9MTsbKGteHLu9um0zrZ3de/1UnNn36RrdlLJizc7X/vYqk2C/X0QOje6bi8hIsjhvr/y1oL4pLsDXDKpoN6tUifmZfd9aXsmWyqM9ea9E7Nn3mci3372f3GqJW8y187W/vcok2F8LvCgiK7Gy54cBhd8cVEnVDKpg/qyTCz0MIL5qXr4V88LOnn2fdxc8DzkLHeV77t7t7Hw7N1viep2dr1QhZZKN/7yIHA4ciRXs39fiNipT9qp5blTMU95z9rvP59y9ZucrFQxJg72InGaMeUFEznc8dKiIYIz5fx6PTRWJsFTNU+7T7Px4WmQnM9ry1n2pZvZfBF4AJiV4zAAa7FWo2TPwi6n2vQomLbKTGW15642kwd4Y87+RT28yxqyyPxY5O69UqNkz8Iup9n2Yj9o5a+HbS+WGnRbZyYy2vPVGJgl6fwKcKcpPYrWlVcozfhTR8SoD3/fa97aKeq1rL2bvvmo6DekLhOuonbMWflVVVbtjY+m4lZ2vVDFJtWd/FHA00Muxb18BdPV6YKpAXGxpm6+w96v3rfa9s6JeWwudOjcxcNYE7987B+my8+218LPlZnY+FO+5e1V6Us3sjwTOBXoTv2+/E/iml4NSBeRyS9t8eVVEJwjs5XPBqrJ38fgcujc6K+rd+AcXRueNMGXnF/O5e1V6Uu3ZLwAWiMjJxpjgHgpW7ivRlrZ+chY6ali/AyC3YO8Dt4rshCk7X8/dq2KSyZ795SKy3BizHUBE+gC/MMb8p7dDU6p4XTx+aFxgL6UiO8Wiqelxmjc8Fbvd0tJAeXlNAUekVHKZBPtjo4EewBizTUSO93BMSnnCy2Y3gcq+3zOQTl3c6xnlZpGdYtK84am4AF9eXkNVZaKTykoVXibBvoOI9DHGbAMQkb4ZPk+VIHt5XAhWiVwvm934nn1v067RTZeNdK94z5f3LnXl5TWMGzu30MPIiXbBKy2ZBO1fAP8SkWiK9oXAzd4NSYWZvTwuELgSuV42u/Et+z6BuEY3D/2gIGPIlT07f8/gQZT16+fq69uP4gXtGJ69op6f1fS0C17pyaQ2/qMishQ4Fas2/vnGmAbPR6ZCK9fyuH6cq1fB4szOP9DamuTK3NhXbvI9huc2e0U9v6vpaRe80pPRcrwx5j0R2UTkfL2IDDXGrPV0ZKrkhP1cvcrex4eOZNlpp8Zub1+5kr5tba69vv0onhtNcuzn7nuNXkn3Xr1zfi17RT2tpqe8ljbYi8hkrKX8wcBGrBa3y7EK7ijlqmI+V+8rW0U9wKqVUBu8xDpnxby+bW2M3NlS4FEl5szBaNv9WWEGolQOMpnZ/xg4CfiHMeZ4ETkV0I4EqqQVMvs+LWcRpOZIrXkXg71b5+4hvmJemM7dP/f/FhduMEplKZNgv9cYs0VEOohIB2PMiyJym+cjUypPXh61K2T2fdpGN86KevYZvgv03L1S4ZNJsN8uIuXAy8AfRGQjsM/bYalszF2ylgX1TYBVia1mUEWBRxQMXh61A/ez7zMtn9vuqJ3PjW703L3ym/a3z18mwX4K8BlwDXAJ0Au4yctBqewsqG+KBfmaQRXtSrGWMi+P2rkp2/K5cUftiky6Rjn50I544aP97d2RMtiLSBmwwBhzOnAAeMSXUams1QyqYP6skws9DJWjMJXP9ZKXjXLc7oin/KH97d2RMtgbY/aLSKuI9DLGfOrXoJRSpcnLRjludsTzgr3ADvhbZEcVv0yW8XcDy0Tk78Cu6J3GmG97NipVErSIjgq7tt2fudLv3l5gB/wvsqOKXybB/pnIhypGdQ9ZPeyj7L3sM5BPLfwwFdEp9FE7ewZ+u+x7VRDRgjrR0/b59Lu3F9gBLbKj3Jc02Eer5BljdJ++mC17Mj7AV41uf047hXxr4YeliE4hj9pBfAa+39n3bqqrq2PZsmWx2/aCOkHnbGlrOjUzYHgNE8+3TmSEvd+9NsYpbqlm9n8BxgKIyJ+MMf/uz5CU76pGw4zcF29yrYXvNi/P1UNhG91Anhn4HlfUsxfZSVVgx1kxr6qqql1TlqAq5pa22hin+KUK9mL7fKTXA1EqX16fqw8tjyvq2VdyMimwY6+Ylwkvj+JlK8wtbVPRxjjFL1WwN0k+VyqwwnKu3lceV9SzF9lxu8COl0fxQM/dq9KRKtgfJyI7sGb43SKfE7ltjDFapk0pD9kr6n1r/T76l3dhYIHH5Dcvj+LpuXtVSpIGe2NMmZ8DUUod5Kyo17pnP5vZU6DRFKegn7tXyk0Z9bNXqtQU+qjdZDpxOj1itz+VfTT59u5KqWKjwV75JkxFdIJ01A6gqTO81aMD/o3AHWE+apeLTatXuVJkB7SinnKXBnvlG7eL6JTSUbv/Cmmt/DAftcuW84/BfIrsaEW91LQLXvY02CtfuVlER4/aBZP9zH3l1kr69uqb1VG7sDr29IlxgT2fIjtaUS857YKXGw32Kiv5lMf1gh61CxZn9cTWfa0H68m6JEjn7pX/tAtebjTYq6zkWx5XBYRHFfXsZ+4BZr/nbglZr8/dq4Ps5XO1dG74abBXWQtKeVw3FTr73t7oBjxuduNxRT0veXnuHrTITpQ9r0JL5xYHDfZKEbzse0+b3XhcUS+stMjOQfbyuVo6tzhosA+hueXT5XgAAB6RSURBVEvWsqD+4KnrhvU7qBmkBQ3zFaTs+0TsFfWmjKnm4vFD/RpaSdAiO6qYabAPoQX1TXEBvmZQRbuKa0nl2b++kLw+ahdk9n/fhvVW5WoN9t7q1O/vdOrzT5a+2RMgruOdUmGjwT6kagZVMH/Wydk/Mc/+9YVUykftLh4/NBbcp4X0zH3YdOrzT8q6rQas/1dyaWnrZpEdpfKhwb4U5dm/Phv2qnluVMwrlqN2vibk+chZMa9za2faurcVcET52f/Z8Jxb2rpZZAe0op7KjwZ75Sl71bx8K+YVE18T8nzkrJjX1r2Nj7t9HNf69uyRZ6fsd5+toJ67d7PIjlbUU/nSYK8852bVPDfZj9v5fdQO0ifkhVVVVVWsYt4THzzBhpUH8yxWbF0B4FqwL5Vz91pRLzV7+VwtnZuYBntVsuzH7fw+ahc4PhXZsc/w3eD1ufu9e/brufuAs5fP1dK5yWmwVyWt0MftAiHERXa81K2iM+xoY2/kdimfuw8ye/lcLZ2bnAZ7lVLQauErD2iRnYR69OpCj15dGDfFSgjVc/cqzDTYq5QKWQu/lM/Vh4k9A7/Y+9WXKnudfNBa+WGkwV6lVaha+MV0rt7to3b2anpQ2Ip69gz8Yu5XX6qc/55aKz+cNNirQCuWc/VuHrVzVksMQkU9ewa+Si6MRXbsdfJBa+WHlQZ75Sp7ER1wp5COWwrd2c6to3b2anqgFfXc0tT0OM0bnorddrs8rhbZUYWkwV65yl5EBwhUIZ1Cd7ZTlhVbVwSyyE7zhqfiAnwu5XFT0SI7qpA02CvXBbWIDuhRu0JzJncGrchOeXlNzuVx/aRFdlS2NNgr5YHQ175PUWTHWf8+mwz8sBfZUSqsNNgr5YFQ175PU2THWf++lDLwNze2aEU9FUoa7ItdiPrXF9u5+tDWvs+gyE4pZt87j31qRb1gstfJB62VH6XBvtiFqH+92+fqC519r4rL0adUxwV2ragXPPY6+aC18u002JeCLPrXF7o8rpvn6ksp+z5IRXaUKhR7nXzQWvl2GuxDYO6StSyob4rdbli/g5pBFZ68VyHL43qhFLLvg1hkR2XGXmQnDAV2orR8bvh4GuxFZCLwK6AMuN8Yc6vj8S7Ao8A4YAswzRizWkT6AU8CJwAPG2Ou8nKcQbegvikuwNcMqmj3C95N2ZTHDXIRHT8VMvtei+ykZj93bz9z73URnXTsq0z5FtgB/4rsaPnccPIs2ItIGfBr4AygEXhDRBYaYxpsl30D2GaMOUxEpgO3AdOA3cD1wDGRj5JXM6iC+bNOLvQw2glyER0/hTr7vsC8LLJj/1l0nrn3uohOOvYiO/kU2AF/i+xo+dxw8nJmfyLwkTFmJYCIzAOmAPZgPwWYE/n8SeBuERFjzC7gnyJymIfjUy4JchEdP4U2+z6Nup0DWLatC9zy3wA0t3Wnqk8PV17b6yI79nP3ic7ch6WITjpaZEel42WwrwY+sd1uBMYnu8YYs09EPgX6AZszeQMRuQy4DGDoUN2fVMoLyziSZnZRFZk7VrGJ0TS68tpeF9lRSlm8DPaS4D6TwzVJGWN+C/wWoLa2NuPnqeKkR+080rOKqp4cPFef4Nx9qdIiOyosvAz2jcAhtttDgHVJrmkUkY5AL2Crh2NSAeJ2EZ1SOmqnCk+L7Kgw8TLYvwEcLiIjgCZgOuA89LgQuBR4DbgAeMEYozP0EuF2ER3w76hd0Gvf53PuPp/a96VEi+yEg1bUs3gW7CN78FcBf8U6evegMeY9EbkJqDPGLAQeAB4TkY+wZvTTo88XkdVABdBZRKYCZzoy+VURcLOIjp+CnH2f77n7Uq59r4qLVtQ7yNNz9saYRcAix3032D7fDSRMuzXGDPdybErlK6jZ926cuy/F2veFZi+wA1pkxw1aUe8graBX4rItj6tFdFSY2Qvs7J64nI79+hV4RBZnbokW2VFu02Bf4rItj6tFdBSQst99vrwqsuP8OT3Q2sq+vF/VHfYCO6BFdpT7NNirrMrjQnCK6Ph51C7oCXm+StPvPh9eFtlp/fx+Nh9+MLzv22DoXKRnf7TIjnLSYF9sAty/PsxH7YKckJevrLPvM+h3nysvi+w4y+N23tqDnivd/TfUc/cqqDTYF5sA968P81E7CG5CXr6KOfve3uwmGuij5XHX3NG+fG4+9Nx9OJTqUTwN9sUoi/71fgvrUbuwS3fuvliz7+2zea8b3ei5++Ar5aN4GuyVKnLOc/efrf+Qt7a8zp4Gq2VysRfNKZZmNyp/pXwUT4O9UgkUU0Ke89z9tT99na57d2LVrCquZXulVGIa7FVo+Jl9X8wJeQC7O/UsymX7bNnP3YN1PC/aErfQ3C6y49e5ewhukZ1SpsE+gOYuWcuC+qbY7Yb1O6gZVFGQsQSpiI7fjW6KNSHPE/Zz9y6euQf/zt3vfv99gEAEe7eL7Ph57l6L7ASTBvsAWlDfFBfgawZVtNt3zVW2FfOCVkTHz+z7YuE8Wtd17052d+rp3hvYT3u4eOYevD1332faV+ICu32GX2huF9nx89y9FtkJJg32AVUzqIL5s052/XWzrZgH+RXRsZ+tz/dcvcqN82jd7k492dHNxRmd/dy9y73uvTx3r1Qp0WBfgrKtmJcP+9l6N87Ve6WYEvISsR+tm3bfa3m1wA0y+7l6IK6ITiFokZ3gK5Vz9xrslefCcLa+2BPy7PJtgRtkzip5Xp+tT0WL7ARfKZ2712CvVESpJOS50QI3yIJyrl6L7ARfKZ2712CvAsvPo3bFJOta90p5TI/iFZ4GexVYfh+1KxYFr3XvYftbr+m5ez2KV6w02BcDe6e7AHW5c4OXR+3sSXnFnJDnKw/b30bpuftwnbvXo3jBoMG+GNg73eXZ5S5IRXS8Zk/KC3NCXqCW7T1sfwvtz933/Owt9q59m6Utfwbyy77Xc/eqmGmwLxYudbrLt4iO2z3rvVYMSXluLNuH5Sie89z93H+MpY+0xG4XMvteqSDTYF/ksq2YB/kV0fGiZ72Kl2wmn+uyfdiP4m0z5ZwVgOz7TOi5++Ar1nP3GuyLXC4V8/KV67l6zb7PjNsJeMV+FC8o9Nx98BXzuXsN9iXAz4p5+fA6+76YquQVLAGvwJwV8vpIC9tMOP4N9dz9QUE9ilfM5+412KtA8TL7vpSq5AVKHkfxnMF9+/YlAPTuPR6wlvBXHRjo3lgd9CieHsUrFhrsVUkJa0KefZ8+VEVy8jyK5yx/27v3+Lg91Lufm8GKrSv4SI/iAXoUTyWnwV6pELDv0/tRJMe17HwXjuKlKn+rLXAtehRPpaPBPgDmLlnLgvqm2G17L3uv5XuuPshH7cK8R+92xn02Cpmdn23XOm2Bq1RmNNgHwIL6prgAXzOoot0vXK/ke64+n6N2Xmffh3mPvpAlbwuZnR+krnV+sB/F02N4wVQsR/E02AdEzaAK5s86uSDvnc+5esj9qJ0fte/DukcPpZtxH5SudV6z/1Fc6sfwgpqdX0xH8TTYF5lciugUkpfZ9yqgQtwox032o3hBPIZnz873KjMfgp2dX0xH8TTYh5G98Q3ENb8pRBGdoCjGPfqi48jOb9r/Ls0b18Cbfwfyq20f5VWjHCido3j27HwvM/NBs/P9osE+jOyNb6Bd85uwFNFxm+7Re8et7PymQV1p7tArdnv79o7AZ/SO3M53j97L7PxSOopnz87XzPzioME+rFxqfFNswrJHX8hs+2y5mZ3f7tz87m5U7epF9Wnu7NF7mZ1fSkfxVHJhTdjTYF9iCnnUTmvfHxT0mbxdPtn5yY7SxRLwXG6BG2baJOcgTdhznwb7ElPIo3ZuZ9+HeY8eSiPbvtSO0uVKm+QcpAl73tBgX4IKddQO3M2+D/MefSlJe5QuxNn5biXsBb1Jjl9180ET9ryiwV6FWtj36MPKtXK6edbOz4Q9O9/NzPwgJ+yBe9n5ftbND6Ow7OFrsC8C9rP1QTtXb9+n1z36cOzRp5MqYS/bcrdu1M5PxZ6d72ZmPgQ7Yc/N7Hytm59cmPbwNdgXAfvZ+qCdq7fv05fSHn2Ysu2z5UzY+9/Hb+awnv9i6Zu92rWgLfQevT07P2x18/NJ2Cu27PygJuyFaQ9fg30BeNH4xquz9W40unFrnz5Me/TFNJNP57Ce/6Jfl7XA6HYtaFVuNGHvoCAn7IWJBvsCyLrxTYqKeenke9Qun+x7L4Rljx6KN9veuVTfr8tatuwZytlu1bPXhD1N2LMJW8JeUPfwNdgXSFaNb9JUzEsl36N2kF/2fT6KYdm+GDiDu3OpfsueobzSOJanIwl7OSfrgecJe16W0tWEPU3YC/Ievgb7sMijYl6+R+2y4WbhHF22D4Z2Ve8cS/Ur9q7l0yZrWyqf6nqApwl7XpbSBU3Y85Pu4WdPg30IBbmznduFc4K6bF9MCXjOmbtTu6p3DvaEvWyq6/nNy1K6ftCEPYvu4edGg30IednZrpAJeWFati+mmbxz5u5U6Iz6MO/huyXoCXu6h5+cfQ+/Z/kojjji+oKMQ4N9SCXLvvc7Ia9Ulu2LeSafbuaeLdcK7kCo9/Ah2Al7xbSHb1/WD8qSPrTfwy8kDfZFxu+EvFJZti/mmbybM3c3O+QBod7DD3LCXjHt4dv/Pwzakr5zD7+QNNgXIT8T8qA4l+2LaSafiJszebt8OuRlzKVlfa/38L1O2AvyHn6hlvWDvqRfSBrsgyiPc/XZcmOPPldBWrZ3Bvc1a9YAMGzYMCBcM/lME+78osv6B7m1rB/kPfxCL+sHNVO/0DTYB5HjXP0TVSNYVH4AIr+E3My+93OPPtlMPgjL9s5l+mHDhoX2l0SQEu50Wf8gN5f1g7yHX8hlfc3UT06DvQ9yKo9rO1e/6LkZVoDHCkTO7Ht7Ul62CXng3x59kGbyiYR1md7rhLt86LL+QUFe1ndzDz8Rv5b1w5ap7ycN9j7IujxuAqlq39uT8nJJyMtWNnv09tl8kGbyEL90H6aKd+kq2hX8qFwauqx/kH1ZP9clfch/Wd/LPXxd1g8GDfY+yao8bg4yTcrLdo8+36N19tl8oWfyqfblw7Qnn66iXZDpsv5B9j/K883U93pZP9clfdBl/aDQYB9AT9DCItkViD36bJftC7kv7wzmTs6kuzDty9tn80Faps9W2Jf13Zzp25f1vSit69ayvttL+lDYZf1SnelrsA8CR/b9or0bWdG5E9HwnmqPHlLv0yebySfbo082k0+2bO8M7m2rPgWg84hegL/78s4kO6cwBXcn+2w+6Mv02QrLsr7XM323MvXB3WX9+TfOdi15Dwq7rF/KM30N9kHg7GrXuQdH9hiQ0R49pC6c4/VM3pl013lEL0/P1NoV01n4fOvTh5Uvy/r2mX4epXa9nOm7XYAn0bJ+rjN9t5P3Ei3r60zfexrsPZBv9j0ZZAFnUzjHy5l8IZfpw3QWPl0wdybZORXbbD4q0bK+ZzN9l5P33JzpJ8rUD8pMP1HyXrHO9NesWcOaNWvifs8US/DXYO+BbLPv0+3Ru7FsH+UM7o0N7wIwpOYYoP1M3s9l+mLec0939j1MSXZecv5/smTVVpas2hr3x3NWwd8+03fO8iE0M/3WN96g9Y034n4PZBP8daafmHOmn2i1MHpd2IkxptBjcEVtba2pq6sr9DCAg0lHSbPvHXv0M/autPboq8bF7rP/koj+lW8P8Pb/0e0Bft2H2wEYfHhvAHZsrGN/2/v06N0ZaB/cK/cdwrAeR9OjT9+EQ31nbQMflTXToevBvwuPPnwU/3bhlzL5VqQN4HbOYJ5IoYJ7upk5/P/27j/IqvK+4/j7s7uyy/5QE1wgiiI00fEHVgUhKTaOUgtMjMbUjCQm06lptK1Ja51OGqd2xmamTf/oTH8mbVrUmCbSqkWlHRJ/xBhi4i8EEgGJqQqIKCCKsiu4sHz7xzkXzt7u3r337t69ew+f18zO7jn3ued59uze+3me55x7DgPCejx/9r2RFM+SPfXymwDMm3Hk/7Xs8C++MuWWx5Pv0y88sm4E4V986+nVO5L3ozlTjvy/Vhv+Q3X4q70sdvGgoPh9A8oP/8JIP/spnZGM9IcbYEBtL7+bVZjmz54HNN4GGJKejYhhG+Swr4Hhwv6eOy5i5YGdMKEDgF/Qx+kdJ3HHp74HDP/CLn6hbl2/iv6+TbRObKH77L20nfwiUvJYBOjVs5n5+pEXXvuxxx8O9+Iw75r8cyZP3kxHUxsAbxx8m9d3Tae55WIA3nvvYbonb6atrfXw9np7z6a3Z/DDBOUEeFb2hTRYwNZq9DvSafbixwcr75H7yDn806YOMwCoVKUn8mYNN1sIjRv+wx06hPqHv8N+DA11jP5w2A8zkv/Vn+7gwo1B98RkOvyX2yawY8ocmrqSwwA7Wn9JT8d+1NQEwARa+ZX+qczsnwTAS007eHXCXpqam5h67n10dr5JX88JSdnON+jpeT+7N31m0LZv73sDOPLPe+xx/0Zn55scPJi8cNratgL/P8T27z9lwOOFZRgY/h2da5g8eTNdXV1A6cAb7mIxwwVupdvPln92zWeGvWZ8tdu22slT+FcS/MUDgnefeQaA9gsuGPI5lXQGCtP8Q12Do9Sov97hX8tRfznhX6zWnYFxEfaSFgF/DzQDSyPir4sebwW+DcwGdgNXR8Tm9LGbgc8D/cAfRsSDpeqqZ9gXTioqHKNf8O5Krmj+KVO6ktHxPW+uZWVHB7Ql/4yL1l3NGU2zaFIzAFum38+BE9dyiEMANNHE3p2nsWf7WQDsbukFYNLBjgHL3YeS5eYT19I9eTOtTccwoX0XLTGdSxYmswSP/ejj9Pb2sHPH4J/jLQ7jPXueY//+Ew6XL34cSk9XlwroSkfHw9VVzvNLPZ5d9jR7PpQT/lklOwLlhH+xCjoD2fAfbNRfrFRnoDj8i1XaGSge6WcNNuVfLNsZKCf8i1XSGciG/2Cj/mKj2Rmo9DyjqVOnsnjx4lGpu6DuYS+pGXgBuBTYBjwDfDoiNmbK/AFwTkT8nqQlwJURcbWkM4FlwFzgROAR4LSI6B+qvrEM++I3lLk/uZX2noOHlyfNn0zbqa8MeM7bOz7Inu3JqHHSSS/QNeUFIpJw73jfawD0vpWOprt2815PN1vXJWdfH/eBjUydupXOY5Jw7z20n33qQ83JXH3x6DsbkMONWEd7urnUCLfS496jUXep7ReX90g8f4pfq1nDdQSKLXh3JfP3/XDAuhM6Ww936svqDAzh8Em6WR3d0JUcKy6nM1BK8exhsXI6AwVbWk7j1ZaZA9a1TJpES3ey7eE6A+/sXE3P7qGvord/72agdGdgKFMOnkx3/8AZh+xhy3I6A5Uq1Xko7gzkNew/AtwaEQvT5ZsBIuJrmTIPpmWekNQCvA50A1/Jls2WG6q+0Q775bdfQctxu6t6bnF4D7W8d89kAEQTE3fP5dxDyYjgpZNuYX/7FrqOT0b2I5m+rnXAmjWqUh2BchR3FgbrDFTrrL4kIDZMSD4a9kh7Hz9p76t6e8+3JuOkM95rHvTx2T/vZ9amIcdSJc3YlmTIy9OSwcc7x86nt6v6aeuepn30922iqX9f1dso6I/kPbxZySHPGZ2nc0rnB0e83YLutuQ9fdf+18oq/86Bt7jo69ePWv0wPsL+KmBRRPxuuvw5YF5EfDFTZn1aZlu6/CIwD7gVeDIivpOuvw34XkTcW1THdcB1AKeccsrswpTJaKg87IUKZ8UBx+y7kEXX/BVQ+fFcH/81G/9G2lkoZTQ7DjDyzkIpI+koDGaknYWs9w5uoa+/Nn8jqLzzkNew/xSwsCjs50bElzJlNqRlsmE/F/gq8ERR2K+MiP8aqr7xdDa+mZnZWCg37Jtq2IZtwMmZ5WnA9qHKpNP4xwFvlvlcMzMzK0Mtw/4Z4EOSZkiaACwBVhSVWQH8dvrzVcCjkUw1rACWSGqVNAP4EPB0DdtqZmaWWzW7XG5EHJT0ReBBko/e3R4RGyR9FVgdESuA24B/l/S/JCP6JelzN0i6G9gIHARuKHUmvpmZmQ3NF9UxMzNrUOPhmL2ZmZmNAw57MzOznHPYm5mZ5ZzD3szMLOcc9mZmZjnnsDczM8s5h72ZmVnOOezNzMxyzmFvZmaWc7m5gp6kXcDo3eN29JwAvFHvRjQg77fqeL9Vx/utOt5v1RnN/TY9IrqHK5SbsB+vJK0u51KGNpD3W3W836rj/VYd77fq1GO/eRrfzMws5xz2ZmZmOeewr71/rXcDGpT3W3W836rj/VYd77fqjPl+8zF7MzOznPPI3szMLOcc9jUiqU3S05J+JmmDpL+od5saiaRmSWsl/U+929IoJG2W9JykdZJW17s9jULS8ZLulbRJ0vOSPlLvNo13kk5P/88KX+9IurHe7WoEkv44zYT1kpZJahuTej2NXxuSBHRERI+kY4DHgT+KiCfr3LSGIOkmYA5wbERcVu/2NAJJm4E5EeHPPVdA0p3AjyNiqaQJQHtE7Kl3uxqFpGbgVWBeRIzHa52MG5JOIsmCMyNin6S7gZUR8a1a1+2RfY1EoiddPCb9cs+qDJKmAR8Dlta7LZZvko4FPgrcBhARfQ76ii0AXnTQl60FmCipBWgHto9FpQ77GkqnotcBO4GHI+KperepQfwd8GXgUL0b0mACeEjSs5Kuq3djGsRMYBdwR3rYaKmkjno3qsEsAZbVuxGNICJeBf4G2Aq8BrwdEQ+NRd0O+xqKiP6IOBeYBsyVdHa92zTeSboM2BkRz9a7LQ1ofkScDywGbpD00Xo3qAG0AOcD/xwR5wG9wFfq26TGkR72uBy4p95taQSS3gdcAcwATgQ6JH12LOp22I+BdFrwMWBRnZvSCOYDl6fHn/8DuETSd+rbpMYQEdvT7zuB+4C59W1RQ9gGbMvMut1LEv5WnsXAmojYUe+GNIjfAF6OiF0RcQBYDvzaWFTssK8RSd2Sjk9/nkjyR95U31aNfxFxc0RMi4hTSaYHH42IMen5NjJJHZK6Cj8Dvwmsr2+rxr+IeB14RdLp6aoFwMY6NqnRfBpP4VdiK/BhSe3pSdwLgOfHouKWsajkKPUB4M70TNUm4O6I8MfIrFamAPcl7x+0AHdFxPfr26SG8SXgu+mU9EvA79S5PQ1BUjtwKXB9vdvSKCLiKUn3AmuAg8Baxuhqev7onZmZWc55Gt/MzCznHPZmZmY557A3MzPLOYe9mZlZzjnszczMcs5hb2ZmlnMOezMzs5xz2JvliKSe4UuNaPsTJf0ovVgUkq6X9FrRvc1npY99U9L89Ocpku6S9FJ6o54nJF05TF2PSVpYtO5GSd+QtCq9a5iZlcFhb2aVuBZYHhH96fI5wC0RcW7m67n0sXnAk+llQe8HVkXEzIiYTXIp5GnD1LUsLZdVuMPaD4CrR+H3MTsqOOzNckjSTZLWp183Ztb/uaRNkh6WtEzSn1S46WuABzLLs4B1g9R/BvBC2im4BOiLiH8pPB4RWyLiH9Oyn5X0dDor8M3CrAHJTWkuk9SaljuV5E5hj5N0Hq6psO1mRy2HvVnOSJpNcn33ecCHgS9IOk/SHOC3gPOATwJzKtzuBGBmRGzOrD6L5F7whSn869L1i4HvZ8qsGWKbZ5CM0Oent4PuJw3xiNgNPM2Ru0UuAf4zkmt8rwcuqKT9ZkczH/Myy58LgfsiohdA0nLg10k69w9ExL50/X8XniBpJvBnwHERcVV657xvAH3AYxHxXeAEYE/mOScDOyPinEHasJAhbigj6etpG/uAO4HZwDPpTXwmAjszxQtT+Q+k368FiIh+SX2SuiJibwX7xuyo5JG9Wf6owvVExEsR8fnMqk8C90bEF4DL03X7gLZMmXMY5LbN6d3Qjo+I7emqDWTuER8RN5Dc2rM7bdOdmeP9p0fErZnN3Q8skHQ+MDEisjMErcD+oX4nMzvCYW+WP6uAT6T3zO4ArgR+THKs++OS2iR1Ah8rsY1pwCvpz/0AEfEW0CypEPizGCTsgYuBH2aWHwXaJP1+Zl17+v0HwFWSJgNIer+k6YVCEdEDPAbcTua+6ZImAbsi4kCJ38HMUp7GN8uZiFgj6Vskx7sBlkbEWgBJK4CfAVuA1cDbQ2xmG0ngr2PgoOAhkin4R0jC/iJJiwtVkxwuWExycl2hPSHpE8DfSvoysAvoBf40IjZKugV4SFITcAC4IW1fwTJgOQPPzL8YWFnWDjEz38/e7GgiqTMietKp9lXAdWnnYBLwl8ClwFLgH4B/Ipkmfzw9Zo+k84CbIuJzJepYA8yr5ag7PQ/h5oj4Ra3qMMsTh73ZUUTSXcCZJMfe74yIr1WxjWvT5/YPW7gG0k8FLImIb9ejfrNG5LA3MzPLOZ+gZ2ZmlnMOezMzs5xz2JuZmeWcw97MzCznHPZmZmY557A3MzPLOYe9mZlZzjnszczMcu7/AGK0aD4/XbsqAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(figsize=(8,6))\n",
    "\n",
    "bin_mids = (bins[1:] + bins[:-1])/2.\n",
    "for ebin in ebins:\n",
    "    ax.step(bin_mids, e_dists[ebin], label=ebin)\n",
    "\n",
    "ax.set_title('True Energy Distributions for Energy Bins')\n",
    "ax.set_xlabel(r'$\\log_{10}(E$/GeV)')\n",
    "ax.set_ylabel('Fraction of Events')\n",
    "ax.legend()\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9d3746a6-4b3d-4d11-b30c-bdd958725af1",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "py3-v4.1.1: icetray/stable",
   "language": "shell",
   "name": "py3-v4.1.1_stable"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
